Mathematical Model of Stem Cell Differentiation and Tissue Regeneration with Stochastic Noise

Differentiation and self-renewal of stem cells is an essential process for the maintenance of tissue composition. The promise of novel medical therapies combined with the complexity of this process encourage us to employ numerical and mathematical methods. This will allow us to understand better the mechanisms which regulate stem cell behaviour. Perturbations to the cellular environment may have an influence on the death rate, proliferation rate and on the fraction of self-renewal at every stage of differentiation. In this paper, we present mathematical study of the effect of stochastic noise on the process of tissue regeneration. Here, a system of Itô stochastic differential equations with linear diffusion coefficients that is based on a deterministic model of multistage cell lineages is investigated. Numerical simulations of the stochastic model are shown for a different number of stages of differentiation. Interactions between the noise, added to the different stages, are characterised using numerical simulations. The long-time behaviour of the two-dimensional version of the model is fully characterised; asymptotic stability of the related Markov semigroup is proved using the theory of the Markov semigroups and the method of the Khasminskií function.     

Deterministic Versus Stochastic Cell Polarisation Through Wave-Pinning

Cell polarization is an important part of the response of eukaryotic cells to stimuli, and forms a primary step in cell motility, differentiation, and many cellular functions. Among the important biochemical players implicated in the onset of intracellular asymmetries that constitute the early phases of polarization are the Rho GTPases, such as Cdc42, Rac, and Rho, which present high active concentration levels in a spatially localized manner. Rho GTPases exhibit positive feedback-driven interconversion between distinct active and inactive forms, the former residing on the cell membrane, and the latter predominantly in the cytosol. A?deterministic model of the dynamics of a single Rho GTPase described earlier by Mori et al.?exhibits sustained polarization by a wave-pinning mechanism. It remained, however, unclear how such polarization behaves at typically low cellular concentrations, as stochasticity could significantly affect the dynamics. We therefore study the low copy number dynamics of this model, using a stochastic kinetics framework based on the Gillespie algorithm, and propose statistical and analytic techniques which help us analyse the equilibrium behaviour of our stochastic system. We use local perturbation analysis to predict parameter regimes for initiation of polarity and wave-pinning in our deterministic system, and compare these predictions with deterministic and stochastic spatial simulations. Comparing the behaviour of the stochastic with the deterministic system, we determine the threshold number of molecules required for robust polarization in a given effective reaction volume. We show that when the molecule number is lowered wave-pinning behaviour is lost due to an increasingly large transition zone as well as increasing fluctuations in the pinning position, due to which a broadness can be reached that is unsustainable, causing the collapse of the wave, while the variations in the high and low equilibrium levels are much less affected.     

Coexistence of multiple pathogen strains in stochastic epidemic models with density-dependent mortality

Stochastic differential equations that model an SIS epidemic with multiple pathogen strains are derived from a system of ordinary differential equations. The stochastic model assumes there is demographic variability. The dynamics of the deterministic model are summarized. Then the dynamics of the stochastic model are compared to the deterministic model. In the deterministic model, there can be either disease extinction, competitive exclusion, where only one strain persists, or coexistence, where more than one strain persists. In the stochastic model, all strains are eventually eliminated because the disease-free state is an absorbing state. However, if the population size and the initial number of infected individuals are sufficiently large, it may take a long time until all strains are eliminated. Numerical simulations of the stochastic model show that coexistence cases predicted by the deterministic model are an unlikely occurrence in the stochastic model even for short time periods. In the stochastic model, either disease extinction or competitive exclusion occur. The initial number of infected individuals, the basic reproduction numbers, and other epidemiological parameters are important determinants of the dominant strain in the stochastic epidemic model.     

The probability of extinction in a bovine respiratory syncytial virus epidemic model

Backward bifurcation is a relatively recent yet well-studied phenomenon associated with deterministic epidemic models. It allows for the presence of multiple subcritical endemic equilibria, and is generally found only in models possessing a reasonable degree of complexity. One particular aspect of backward bifurcation that appears to have been virtually overlooked in the literature is the potential influence its presence might have on the behaviour of any analogous stochastic model. Indeed, the primary aim of this paper is to investigate this possibility. Our approach is to compare the theoretical probabilities of extinction, calculated via a particular stochastic formulation of a deterministic model exhibiting backward bifurcation, with those obtained from a series of stochastic simulations. We have found some interesting links in the behaviour between the deterministic and stochastic models, and are able to offer plausible explanations for our observations.     

Langevin equation,Fokker-Planck equation and cell migration

Cell migration can be characterized by two independent variables: the speed,v, and the migration angle, ϕ. Each variable can be described by a stochastic differential equation—a Langevin equation. The migration behaviour of an ensemble of cells can be predicted due to the stochastic processes involved in the signal transduction/response system of each cell. Distribution functions, correlation functions, etc. are determined by using the corresponding Fokker-Planck equation. The model assumptions are verified by experimental results. The theoretical predictions are mainly compared with the galvanotactic response of human granulocytes. The coefficient characterizing the mean effect of the signal transduction/response system of the cell is experimentally determined to 0.08 mm/V sec (galvanotaxis) or 0.7 mm/sec (chemotaxis) and the characteristic time characterizing stochastic effects in the signal transduction/response system is experimentally determined as 30 sec. The temporal directed response induced by electric field pulses is investigated: the experimental cells react slower but are more sensitive than predicted by theory.     

A spatially structured metapopulation model with patch dynamics

Metapopulation models that incorporate both spatial and temporal structure are studied in this paper. The existence and stability of equilibria are provided, and an extinction threshold condition is derived which depends on patch dynamics (patch destruction and creation) and metapopulation dynamics (patch colonization and extinction). These results refine threshold conditions given by previous metapopulation models. By comparing landscapes with different spatial heterogeneities with respect to weighted long-term patch occupancies, we conclude that the pattern of a landscape is of overwhelming importance in determining metapopulation persistence and patch occupancy. We show that the same conclusion holds when a rescue effect is considered. We also derive a stochastic differential equations (SDE) model of the It? type based on our deterministic model. Our simulations reveal good agreement between the deterministic model and the SDE model.     

A stochastic model for adhesion-mediated cell random motility and haptotaxis

The active migration of blood and tissue cells is important in a number of physiological processes including inflammation, wound healing, embryogenesis, and tumor cell metastasis. These cells move by transmitting cytoplasmic force through membrane receptors which are bound specifically to adhesion ligands in the surrounding substratum. Recently, much research has focused on the influence of the composition of extracellular matrix and the distribution of its components on the speed and direction of cell migration. It is commonly believed that the magnitude of the adhesion influences cell speed and/or random turning behavior, whereas a gradient of adhesion may bias the net direction of the cell movement, a phenomenon known as haptotaxis. The mechanisms underlying these responses are presently not understood.A stochastic model is presented to provide a mechanistic understanding of how the magnitude and distribution of adhesion ligands in the substratum influence cell movement. The receptor-mediated cell migration is modeled as an interrelation of random processes on distinct time scales. Adhesion receptors undergo rapid binding and transport, resulting in a stochastic spatial distribution of bound receptors fluctuating about some mean distribution. This results in a fluctuating spatio-temporal pattern of forces on the cell, which in turn affects the speed and turning behavior on a longer time scale. The model equations are a system of nonlinear stochastic differential equations (SDE's) which govern the time evolution of the spatial distribution of bound and free receptors, and the orientation and position of the cell. These SDE's are integrated numerically to simulate the behavior of the model cell on both a uniform substratum, and on a gradient of adhesion ligand concentration.Furthermore, analysis of the governing SDE system and corresponding Fokker-Planck equation (FPE) yields analytical expressions for indices which characterize cell movement on multiple time scales in terms of cell cytomechanical, morphological, and receptor binding and transport parameters. For a uniform adhesion ligand concentration, this analysis provides expressions for traditional cell movement indices such as mean speed, directional persistence time, and random motility coefficient. In a small gradient of adhesion, a perturbation analysis of the FPE yields a constitutive cell flux expression which includes a drift term for haptotactic directional cell migration. The haptotactic drift contains terms identified as contributions from directional orientation bias (taxis).     

Mathematical modelling of whole chromosome replication

All chromosomes must be completely replicated prior to cell division, a requirement that demands the activation of a sufficient number of appropriately distributed DNA replication origins. Here we investigate how the activity of multiple origins on each chromosome is coordinated to ensure successful replication. We present a stochastic model for whole chromosome replication where the dynamics are based upon the parameters of individual origins. Using this model we demonstrate that mean replication time at any given chromosome position is determined collectively by the parameters of all origins. Combining parameter estimation with extensive simulations we show that there is a range of model parameters consistent with mean replication data, emphasising the need for caution in interpreting such data. In contrast, the replicated-fraction at time points through S phase contains more information than mean replication time data and allowed us to use our model to uniquely estimate many origin parameters. These estimated parameters enable us to make a number of predictions that showed agreement with independent experimental data, confirming that our model has predictive power. In summary, we demonstrate that a stochastic model can recapitulate experimental observations, including those that might be interpreted as deterministic such as ordered origin activation times.     

Stochastic simulations on a model of circadian rhythm generation

Biological phenomena are often modeled by differential equations, where states of a model system are described by continuous real values. When we consider concentrations of molecules as dynamical variables for a set of biochemical reactions, we implicitly assume that numbers of the molecules are large enough so that their changes can be regarded as continuous and they are described deterministically. However, for a system with small numbers of molecules, changes in their numbers are apparently discrete and molecular noises become significant. In such cases, models with deterministic differential equations may be inappropriate, and the reactions must be described by stochastic equations. In this study, we focus a clock gene expression for a circadian rhythm generation, which is known as a system involving small numbers of molecules. Thus it is appropriate for the system to be modeled by stochastic equations and analyzed by methodologies of stochastic simulations. The interlocked feedback model proposed by Ueda et al. as a set of deterministic ordinary differential equations provides a basis of our analyses. We apply two stochastic simulation methods, namely Gillespie's direct method and the stochastic differential equation method also by Gillespie, to the interlocked feedback model. To this end, we first reformulated the original differential equations back to elementary chemical reactions. With those reactions, we simulate and analyze the dynamics of the model using two methods in order to compare them with the dynamics obtained from the original deterministic model and to characterize dynamics how they depend on the simulation methodologies.     

Geographical patterns of adaptation within a species' range: interactions between drift and gene flow

We use individual-based stochastic simulations and analytical deterministic predictions to investigate the interaction between drift, natural selection and gene flow on the patterns of local adaptation across a fragmented species' range under clinally varying selection. Migration between populations follows a stepping-stone pattern and density decreases from the centre to the periphery of the range. Increased migration worsens gene swamping in small marginal populations but mitigates the effect of drift by replenishing genetic variance and helping purge deleterious mutations. Contrary to the deterministic prediction that increased connectivity within the range always inhibits local adaptation, simulations show that low intermediate migration rates improve fitness in marginal populations and attenuate fitness heterogeneity across the range. Such migration rates are optimal in that they maximize the total mean fitness at the scale of the range. Optimal migration rates increase with shallower environmental gradients, smaller marginal populations and higher mutation rates affecting fitness.     

Stochastic differential equation model for cerebellar granule cell excitability

Neurons in the brain express intrinsic dynamic behavior which is known to be stochastic in nature. A crucial question in building models of neuronal excitability is how to be able to mimic the dynamic behavior of the biological counterpart accurately and how to perform simulations in the fastest possible way. The well-established Hodgkin-Huxley formalism has formed to a large extent the basis for building biophysically and anatomically detailed models of neurons. However, the deterministic Hodgkin-Huxley formalism does not take into account the stochastic behavior of voltage-dependent ion channels. Ion channel stochasticity is shown to be important in adjusting the transmembrane voltage dynamics at or close to the threshold of action potential firing, at the very least in small neurons. In order to achieve a better understanding of the dynamic behavior of a neuron, a new modeling and simulation approach based on stochastic differential equations and Brownian motion is developed. The basis of the work is a deterministic one-compartmental multi-conductance model of the cerebellar granule cell. This model includes six different types of voltage-dependent conductances described by Hodgkin-Huxley formalism and simple calcium dynamics. A new model for the granule cell is developed by incorporating stochasticity inherently present in the ion channel function into the gating variables of conductances. With the new stochastic model, the irregular electrophysiological activity of an in vitro granule cell is reproduced accurately, with the same parameter values for which the membrane potential of the original deterministic model exhibits regular behavior. The irregular electrophysiological activity includes experimentally observed random subthreshold oscillations, occasional spontaneous spikes, and clusters of action potentials. As a conclusion, the new stochastic differential equation model of the cerebellar granule cell excitability is found to expand the range of dynamics in comparison to the original deterministic model. Inclusion of stochastic elements in the operation of voltage-dependent conductances should thus be emphasized more in modeling the dynamic behavior of small neurons. Furthermore, the presented approach is valuable in providing faster computation times compared to the Markov chain type of modeling approaches and more sophisticated theoretical analysis tools compared to previously presented stochastic modeling approaches.     

Individualism in plant populations: using stochastic differential equations to model individual neighbourhood-dependent plant growth

We study individual plant growth and size hierarchy formation in an experimental population of Arabidopsis thaliana, within an integrated analysis that explicitly accounts for size-dependent growth, size- and space-dependent competition, and environmental stochasticity. It is shown that a Gompertz-type stochastic differential equation (SDE) model, involving asymmetric competition kernels and a stochastic term which decreases with the logarithm of plant weight, efficiently describes individual plant growth, competition, and variability in the studied population. The model is evaluated within a Bayesian framework and compared to its deterministic counterpart, and to several simplified stochastic models, using distributional validation. We show that stochasticity is an important determinant of size hierarchy and that SDE models outperform the deterministic model if and only if structural components of competition (asymmetry; size- and space-dependence) are accounted for. Implications of these results are discussed in the context of plant ecology and in more general modelling situations.     

Integrating Biosystem Models Using Waveform Relaxation

Modelling in systems biology often involves the integration of component models into larger composite models. How to do this systematically and efficiently is a significant challenge: coupling of components can be unidirectional or bidirectional, and of variable strengths. We adapt the waveform relaxation (WR) method for parallel computation of ODEs as a general methodology for computing systems of linked submodels. Four test cases are presented: (i) a cascade of unidirectionally and bidirectionally coupled harmonic oscillators, (ii) deterministic and stochastic simulations of calcium oscillations, (iii) single cell calcium oscillations showing complex behaviour such as periodic and chaotic bursting, and (iv) a multicellular calcium model for a cell plate of hepatocytes. We conclude that WR provides a flexible means to deal with multitime-scale computation and model heterogeneity. Global solutions over time can be captured independently of the solution techniques for the individual components, which may be distributed in different computing environments.     

A Minimal Model of Signaling Network Elucidates Cell-to-Cell Stochastic Variability in Apoptosis

Background

Signaling networks are designed to sense an environmental stimulus and adapt to it. We propose and study a minimal model of signaling network that can sense and respond to external stimuli of varying strength in an adaptive manner. The structure of this minimal network is derived based on some simple assumptions on its differential response to external stimuli.

Methodology

We employ stochastic differential equations and probability distributions obtained from stochastic simulations to characterize differential signaling response in our minimal network model. Gillespie''s stochastic simulation algorithm (SSA) is used in this study.

Conclusions/Significance

We show that the proposed minimal signaling network displays two distinct types of response as the strength of the stimulus is decreased. The signaling network has a deterministic part that undergoes rapid activation by a strong stimulus in which case cell-to-cell fluctuations can be ignored. As the strength of the stimulus decreases, the stochastic part of the network begins dominating the signaling response where slow activation is observed with characteristic large cell-to-cell stochastic variability. Interestingly, this proposed stochastic signaling network can capture some of the essential signaling behaviors of a complex apoptotic cell death signaling network that has been studied through experiments and large-scale computer simulations. Thus we claim that the proposed signaling network is an appropriate minimal model of apoptosis signaling. Elucidating the fundamental design principles of complex cellular signaling pathways such as apoptosis signaling remains a challenging task. We demonstrate how our proposed minimal model can help elucidate the effect of a specific apoptotic inhibitor Bcl-2 on apoptotic signaling in a cell-type independent manner. We also discuss the implications of our study in elucidating the adaptive strategy of cell death signaling pathways.     

Screening stochastic Life Cycle assessment inventory models

A screening methodology is presented that utilizes the linear structure within the deterministic life cycle inventory (LCI) model. The methodology ranks each input data element based upon the amount it contributes toward the final output. The identified data elements along with their position in the deterministic model are then sorted into descending order based upon their individual contributions. This enables practitioners and model users to identify those input data elements that contribute the most in the inventory stage. Percentages of the top ranked data elements are then selected, and their corresponding data quality index (DQI) value is upgraded in the stochastic LCI model. Monte Carlo computer simulations are obtained and used to compare the output variance of the original stochastic model with modified stochastic model. The methodology is applied to four real-world beverage delivery system LCA inventory models for verification. This research assists LCA practitioners by streamlining the conversion process when converting a deterministic LCI model to a stochastic model form. Model users and decision-makers can benefit from the reduction in output variance and the increase in ability to discriminate between product system alternatives.     

Stochastic models for circadian rhythms: effect of molecular noise on periodic and chaotic behaviour

Circadian rhythms are endogenous oscillations that occur with a period close to 24 h in nearly all living organisms. These rhythms originate from the negative autoregulation of gene expression. Deterministic models based on such genetic regulatory processes account for the occurrence of circadian rhythms in constant environmental conditions (e.g., constant darkness), for entrainment of these rhythms by light-dark cycles, and for their phase-shifting by light pulses. When the numbers of protein and mRNA molecules involved in the oscillations are small, as may occur in cellular conditions, it becomes necessary to resort to stochastic simulations to assess the influence of molecular noise on circadian oscillations. We address the effect of molecular noise by considering the stochastic version of a deterministic model previously proposed for circadian oscillations of the PER and TIM proteins and their mRNAs in Drosophila. The model is based on repression of the per and tim genes by a complex between the PER and TIM proteins. Numerical simulations of the stochastic version of the model are performed by means of the Gillespie method. The predictions of the stochastic approach compare well with those of the deterministic model with respect both to sustained oscillations of the limit cycle type and to the influence of the proximity from a bifurcation point beyond which the system evolves to stable steady state. Stochastic simulations indicate that robust circadian oscillations can emerge at the cellular level even when the maximum numbers of mRNA and protein molecules involved in the oscillations are of the order of only a few tens or hundreds. The stochastic model also reproduces the evolution to a strange attractor in conditions where the deterministic PER-TIM model admits chaotic behaviour. The difference between periodic oscillations of the limit cycle type and aperiodic oscillations (i.e. chaos) persists in the presence of molecular noise, as shown by means of Poincaré sections. The progressive obliteration of periodicity observed as the number of molecules decreases can thus be distinguished from the aperiodicity originating from chaotic dynamics. As long as the numbers of molecules involved in the oscillations remain sufficiently large (of the order of a few tens or hundreds, or more), stochastic models therefore provide good agreement with the predictions of the deterministic model for circadian rhythms.     

A reduced 1D stochastic model of bleb-driven cell migration

Blebs are pressure-driven protrusions that have been observed in cells undergoing apoptosis, cytokinesis, or migration, including tumor cells that use blebs to escape their organs of origin. Here, we present a minimal 1D model of bleb-driven cell motion that combines a simple mechanical model with turnover kinetics of the actin cortex and adhesions between the membrane and the cortex. The deterministic version of this model is used to study the properties of individual blebbing events. We further introduce stochastic turnover of the adhesions, which allows for spontaneous initiation of repeated blebbing events, thus leading to sustained cell travel. We explore how the main parameters of the system control the properties of the blebbing events and the speed of cell travel. Finally, we derive a further simplification by deriving a Langevin approximation to this stochastic model.     

A computational model of cell migration coupling the growth of focal adhesions with oscillatory cell protrusions

Cell migration is a highly integrated process where actin turnover, actomyosin contractility, and adhesion dynamics are all closely linked. In this paper, we propose a computational model investigating the coupling of these fundamental processes within the context of spontaneous (i.e. unstimulated) cell migration. In the unstimulated cell, membrane oscillations originating from the interaction between passive hydrostatic pressure and contractility are sufficient to lead to the formation of adhesion spots. Cell contractility then leads to the maturation of these adhesion spots into focal adhesions. Due to active actin polymerization, which reinforces protrusion at the leading edge, the traction force required for cell translocation can be generated. Computational simulations first show that the model hypotheses allow one to reproduce the main features of fibroblast cell migration and established results on the biphasic aspect of the cell speed as a function of adhesion strength. The model also demonstrates that certain temporal parameters, such as the adhesion proteins recycling time and adhesion lifetimes, influence cell motion patterns, particularly cell speed and persistence of the direction of migration. This study provides some elements, which allow a better understanding of spontaneous cell migration and enables a first glance at how an individual cell would potentially react once exposed to a stimulus.