Survival probability for a diffusive process on a growing domain

Smoothing in linear multicompartment biological processes subject to stochastic input

Many physical and biological systems rely on the progression of material through multiple independent stages. In viral replication, for example, virions enter a cell to undergo a complex process comprising several disparate stages before the eventual accumulation and release of replicated virions. While such systems may have some control over the internal dynamics that make up this progression, a challenge for many is to regulate behavior under what are often highly variable external environments acting as system inputs. In this work, we study a simple analog of this problem through a linear multicompartment model subject to a stochastic input in the form of a mean-reverting Ornstein-Uhlenbeck process, a type of Gaussian process. By expressing the system as a multidimensional Gaussian process, we derive several closed-form analytical results relating to the covariances and autocorrelations of the system, quantifying the smoothing effect discrete compartments afford multicompartment systems. Semianalytical results demonstrate that feedback and feedforward loops can enhance system robustness, and simulation results probe the intractable problem of the first passage time distribution, which has specific relevance to eventual cell lysis in the viral replication cycle. Finally, we demonstrate that the smoothing seen in the process is a consequence of the discreteness of the system, and does not manifest in systems with continuous transport. While we make progress through analysis of a simple linear problem, many of our insights are applicable more generally, and our work enables future analysis into multicompartment processes subject to stochastic inputs.

The mechanism of pattern transitions between formation and dispersion

On the surface of living organisms, a wide variety of patterns can be observed, some of which change during their growth process. For instance, Pelodiscus sinensis exhibits distinct black patterns on its vivid orange plastron during the embryonic and juvenile stages, but as it matures, the black patterns gradually disappear, resulting in a whitened plastron. This pattern transition is a mysterious phenomenon that forms and vanishes on the plastron, a ventral part with low visibility to both predators and peers. Our research aims to focus on understanding the mechanisms behind such pattern transitions and proposes a model capable of representing pattern formation and dispersion. To understand the changing patterns, we propose a hypothesis based on a reaction-diffusion system with a time-dependent growing spatial domain. This mathematical framework suggests the occurrence of the dispersion phenomenon. Specifically, we focus on the dilution term within the system under the growing-domain condition. While previous studies have investigated the effects of growth domains, this study specifically addresses the role of the time-dependently growing domain effects - change of diffusion coefficient and dilution - in reaction-diffusion systems. Our research sheds light on the intricate phenomenon of pattern formation and dispersion on the surface of living organisms, proposing a natural system based on the effects of growing domain, namely, a model of reaction-dilution-diffusion systems.

Influence of electric field, blood velocity, and pharmacokinetics on electrochemotherapy efficiency

The convective delivery of chemotherapeutic drugs in cancerous tissues is directly proportional to the blood perfusion rate, which in turns can be transiently reduced by the application of high-voltage and short-duration electric pulses due to vessel vasoconstriction. However, electric pulses can also increase vessel wall and cell membrane permeabilities, boosting the extravasation and cell internalization of drug. These opposite effects, as well as possible adverse impacts on the viability of tissues and endothelial cells, suggest the importance of conducting in silico studies about the influence of physical parameters involved in electric-mediated drug transport. In the present work, the global method of approximate particular solutions for axisymmetric domains, together with two solution schemes (Gauss-Seidel iterative and linearization+successive over-relaxation), is applied for the simulation of drug transport in electroporated cancer tissues, using a continuum tumor cord approach and considering both the electropermeabilization and vasoconstriction phenomena. The developed global method of approximate particular solutions algorithm is validated with numerical and experimental results previously published, obtaining a satisfactory accuracy and convergence. Then, a parametric study about the influence of electric field magnitude and inlet blood velocity on the internalization efficacy, drug distribution uniformity, and cell-kill capacity of the treatment, as expressed by the number of internalized moles into viable cells, homogeneity of exposure to bound intracellular drug, and cell survival fraction, respectively, is analyzed for three pharmacokinetic profiles, namely one-short tri-exponential, mono-exponential, and uniform. According to numerical results, the trade-off between vasoconstriction and electropermeabilization effects and, consequently, the influence of electric field magnitude and inlet blood velocity on the assessment parameters considered here (efficacy, uniformity, and cell-kill capacity) is different for each pharmacokinetic profile deemed.

Lévy Walk Dynamics in an External Constant Force Field in Non-Static Media

Based on the recognition of the huge change of the transport properties for diffusion particles in non-static media, we consider a Lévy walk model subjected to an external constant force in non-static media. Since the physical and comoving coordinates of non-static media are related by scale factor, we equivalently transfer the process from physical coordinate into comoving coordinate and derive the master equation governing the probability density function of the position of the particles in comoving coordinate. Utilizing the Hermite orthogonal polynomial expansions, some statistical properties are obtained, including the asymptotic behaviors of the first two moments in both coordinates and kurtosis. For some representative types of non-static media and Lévy walks, the striking and interesting phenomena originating from the interplay between non-static media, external force, and intrinsic stochastic motion are observed. The stationary distribution are also analyzed for some cases through numerical simulations.

The role of mechanical interactions in EMT

The detachment of cells from the boundary of an epithelial tissue and the subsequent invasion of these cells into surrounding tissues is important for cancer development and wound healing, and is strongly associated with the epithelial-mesenchymal transition (EMT). Chemical signals, such as TGF-β, produced by surrounding tissue can be uptaken by cells and induce EMT. In this work, we present a novel cell-based discrete mathematical model of mechanical cellular relaxation, cell proliferation, and cell detachment driven by chemically-dependent EMT in an epithelial tissue. A continuum description of the model is then derived in the form of a novel nonlinear free boundary problem. Using the discrete and continuum models we explore how the coupling of chemical transport and mechanical interactions influences EMT, and postulate how this could be used to help control EMT in pathological situations.

Diffusion in heterogeneous discs and spheres: New closed-form expressions for exit times and homogenization formulas

Mathematical models of diffusive transport underpin our understanding of chemical, biochemical, and biological transport phenomena. Analysis of such models often focuses on relatively simple geometries and deals with diffusion through highly idealized homogeneous media. In contrast, practical applications of diffusive transport theory inevitably involve dealing with more complicated geometries as well as dealing with heterogeneous media. One of the most fundamental properties of diffusive transport is the concept of mean particle lifetime or mean exit time, which are particular applications of the concept of first passage time and provide the mean time required for a diffusing particle to reach an absorbing boundary. Most formal analysis of mean particle lifetime applies to relatively simple geometries, often with homogeneous (spatially invariant) material properties. In this work, we present a general framework that provides exact mathematical insight into the mean particle lifetime, and higher moments of particle lifetime, for point particles diffusing in heterogeneous discs and spheres with radial symmetry. Our analysis applies to geometries with an arbitrary number and arrangement of distinct layers, where transport in each layer is characterized by a distinct diffusivity. We obtain exact closed-form expressions for the mean particle lifetime for a diffusing particle released at an arbitrary location, and we generalize these results to give exact, closed-form expressions for any higher-order moment of particle lifetime for a range of different boundary conditions. Finally, using these results, we construct new homogenization formulas that provide an accurate simplified description of diffusion through heterogeneous discs and spheres.

A free boundary mechanobiological model of epithelial tissues

In this study, we couple intracellular signalling and cell-based mechanical properties to develop a novel free boundary mechanobiological model of epithelial tissue dynamics. Mechanobiological coupling is introduced at the cell level in a discrete modelling framework, and new reaction-diffusion equations are derived to describe tissue-level outcomes. The free boundary evolves as a result of the underlying biological mechanisms included in the discrete model. To demonstrate the accuracy of the continuum model, we compare numerical solutions of the discrete and continuum models for two different signalling pathways. First, we study the Rac-Rho pathway where cell- and tissue-level mechanics are directly related to intracellular signalling. Second, we study an activator-inhibitor system which gives rise to spatial and temporal patterning related to Turing patterns. In all cases, the continuum model and free boundary condition accurately reflect the cell-level processes included in the discrete model.

Diffusion in an expanding medium: Fokker-Planck equation, Green's function, and first-passage properties

We present a classical, mesoscopic derivation of the Fokker-Planck equation for diffusion in an expanding medium. To this end, we take a conveniently generalized Chapman-Kolmogorov equation as the starting point. We obtain an analytical expression for the Green's function (propagator) and investigate both analytically and numerically how this function and the associated moments behave. We also study first-passage properties in expanding hyperspherical geometries. We show that in all cases the behavior is determined to a great extent by the so-called Brownian conformal time τ(t), which we define via the relation τ[over ̇]=1/a^{2}, where a(t) is the expansion scale factor. If the medium expansion is driven by a power law [a(t)∝t^{γ} with γ>0], then we find interesting crossover effects in the mixing effectiveness of the diffusion process when the characteristic exponent γ is varied. Crossover effects are also found at the level of the survival probability and of the moments of the first passage-time distribution with two different regimes separated by the critical value γ=1/2. The case of an exponential scale factor is analyzed separately both for expanding and contracting media. In the latter situation, a stationary probability distribution arises in the long-time limit.

Exact Solutions of Coupled Multispecies Linear Reaction-Diffusion Equations on a Uniformly Growing Domain

Embryonic development involves diffusion and proliferation of cells, as well as diffusion and reaction of molecules, within growing tissues. Mathematical models of these processes often involve reaction–diffusion equations on growing domains that have been primarily studied using approximate numerical solutions. Recently, we have shown how to obtain an exact solution to a single, uncoupled, linear reaction–diffusion equation on a growing domain, 0 < x < L(t), where L(t) is the domain length. The present work is an extension of our previous study, and we illustrate how to solve a system of coupled reaction–diffusion equations on a growing domain. This system of equations can be used to study the spatial and temporal distributions of different generations of cells within a population that diffuses and proliferates within a growing tissue. The exact solution is obtained by applying an uncoupling transformation, and the uncoupled equations are solved separately before applying the inverse uncoupling transformation to give the coupled solution. We present several example calculations to illustrate different types of behaviour. The first example calculation corresponds to a situation where the initially–confined population diffuses sufficiently slowly that it is unable to reach the moving boundary at x = L(t). In contrast, the second example calculation corresponds to a situation where the initially–confined population is able to overcome the domain growth and reach the moving boundary at x = L(t). In its basic format, the uncoupling transformation at first appears to be restricted to deal only with the case where each generation of cells has a distinct proliferation rate. However, we also demonstrate how the uncoupling transformation can be used when each generation has the same proliferation rate by evaluating the exact solutions as an appropriate limit.