A comparison of deterministic and stochastic simulations of neuronal vesicle release models.
Ontology highlight
ABSTRACT:
SUBMITTER: Modchang C
PROVIDER: S-EPMC2892017 | biostudies-literature | 2010 May
REPOSITORIES: biostudies-literature
Publications
A comparison of deterministic and stochastic simulations of neuronal vesicle release models.
Modchang Charin C Nadkarni Suhita S Bartol Thomas M TM Triampo Wannapong W Sejnowski Terrence J TJ Levine Herbert H Rappel Wouter-Jan WJ
Physical biology 20100526 2
We study the calcium-induced vesicle release into the synaptic cleft using a deterministic algorithm and MCell, a Monte Carlo algorithm that tracks individual molecules. We compare the average vesicle release probability obtained using both algorithms and investigate the effect of the three main sources of noise: diffusion, sensor kinetics and fluctuations from the voltage-dependent calcium channels (VDCCs). We find that the stochastic opening kinetics of the VDCCs are the main contributors to d ...[more]
Similar Datasets
Project description:The lac operon has been a paradigm for genetic regulation with positive feedback, and several modeling studies have described its dynamics at various levels of detail. However, it has not yet been analyzed how stochasticity can enrich the system's behavior, creating effects that are not observed in the deterministic case. To address this problem we use a comparative approach. We develop a reaction network for the dynamics of the lac operon genetic switch and derive corresponding deterministic and stochastic models that incorporate biological details. We then analyze the effects of key biomolecular mechanisms, such as promoter strength and binding affinities, on the behavior of the models. No assumptions or approximations are made when building the models other than those utilized in the reaction network. Thus, we are able to carry out a meaningful comparison between the predictions of the two models to demonstrate genuine effects of stochasticity. Such a comparison reveals that in the presence of stochasticity, certain biomolecular mechanisms can profoundly influence the region where the system exhibits bistability, a key characteristic of the lac operon dynamics. For these cases, the temporal asymptotic behavior of the deterministic model remains unchanged, indicating a role of stochasticity in modulating the behavior of the system.
| S-EPMC2996131 | biostudies-other
Project description:Experimental studies of the flowering of Arabidopsis thaliana have shown that a large complex gene regulatory network (GRN) is responsible for its regulation. This process has been mathematically modelled with deterministic differential equations by considering the interactions between gene activators and inhibitors (Valentim et al. in PLoS ONE 10(2):e0116973, 2015; van Mourik et al. in BMC Syst Biol 4(1):1, 2010). However, due to complexity of the model, the properties of the network and the roles of the individual genes cannot be deducted from the numerical solution the published work offers. Here, we propose simplifications of the model, based on decoupling of the original GRN to motifs, described with three and two differential equations. A stable solution of the original model is sought by linearisation of the original model which contributes to further investigation of the role of the individual genes to the flowering. Furthermore, we study the role of noise by introducing and investigating two types of stochastic elements into the model. The deterministic and stochastic nonlinear dynamic models of Arabidopsis flowering time are considered by following the deterministic delayed model introduced in Valentim et al. (2015). Steady-state regimes and stability of the deterministic original model are investigated analytically and numerically. By decoupling some concentrations, the system was reduced to emphasise the role played by the transcription factor Suppressor of Overexpression of Constants1 ([Formula: see text]) and the important floral meristem identity genes, Leafy ([Formula: see text]) and Apetala1 ([Formula: see text]). Two-dimensional motifs, based on the dynamics of [Formula: see text] and [Formula: see text], are obtained from the reduced network and parameter ranges ensuring flowering are determined. Their stability analysis shows that [Formula: see text] and [Formula: see text] are regulating each other for flowering, matching experimental findings. New sufficient conditions of mean square stability in the stochastic model are obtained using a stochastic Lyapunov approach. Our numerical simulations demonstrate that the reduced models of Arabidopsis flowering time, describing specific motifs of the GRN, can capture the essential behaviour of the full system and also introduce the conditions of flowering initiation. Additionally, they show that stochastic effects can change the behaviour of the stability region through a stability switch. This study thus contributes to a better understanding of the role of [Formula: see text] and [Formula: see text] in Arabidopsis flowering.
| S-EPMC6320361 | biostudies-other
Project description:BackgroundGillespie's stochastic simulation algorithm (SSA) for chemical reactions admits three kinds of elementary processes, namely, mass action reactions of 0th, 1st or 2nd order. All other types of reaction processes, for instance those containing non-integer kinetic orders or following other types of kinetic laws, are assumed to be convertible to one of the three elementary kinds, so that SSA can validly be applied. However, the conversion to elementary reactions is often difficult, if not impossible. Within deterministic contexts, a strategy of model reduction is often used. Such a reduction simplifies the actual system of reactions by merging or approximating intermediate steps and omitting reactants such as transient complexes. It would be valuable to adopt a similar reduction strategy to stochastic modelling. Indeed, efforts have been devoted to manipulating the chemical master equation (CME) in order to achieve a proper propensity function for a reduced stochastic system. However, manipulations of CME are almost always complicated, and successes have been limited to relative simple cases.ResultsWe propose a rather general strategy for converting a deterministic process model into a corresponding stochastic model and characterize the mathematical connections between the two. The deterministic framework is assumed to be a generalized mass action system and the stochastic analogue is in the format of the chemical master equation. The analysis identifies situations: where a direct conversion is valid; where internal noise affecting the system needs to be taken into account; and where the propensity function must be mathematically adjusted. The conversion from deterministic to stochastic models is illustrated with several representative examples, including reversible reactions with feedback controls, Michaelis-Menten enzyme kinetics, a genetic regulatory motif, and stochastic focusing.ConclusionsThe construction of a stochastic model for a biochemical network requires the utilization of information associated with an equation-based model. The conversion strategy proposed here guides a model design process that ensures a valid transition between deterministic and stochastic models.
| S-EPMC3236013 | biostudies-literature
Project description:BackgroundThe lac operon genetic switch is considered as a paradigm of genetic regulation. This system has a positive feedback loop due to the LacY permease boosting its own production by the facilitated transport of inducer into the cell and the subsequent de-repression of the lac operon genes. Previously, we have investigated the effect of stochasticity in an artificial lac operon network at the single cell level by comparing corresponding deterministic and stochastic kinetic models.ResultsThis work focuses on the dynamics of cell populations by incorporating the above kinetic scheme into two Monte Carlo (MC) simulation frameworks. The first MC framework assumes stochastic reaction occurrence, accounts for stochastic DNA duplication, division and partitioning and tracks all daughter cells to obtain the statistics of the entire cell population. In order to better understand how stochastic effects shape cell population distributions, we develop a second framework that assumes deterministic reaction dynamics. By comparing the predictions of the two frameworks, we conclude that stochasticity can create or destroy bimodality, and may enhance phenotypic heterogeneity.ConclusionsOur results show how various sources of stochasticity act in synergy with the positive feedback architecture, thereby shaping the behavior at the cell population level. Further, the insights obtained from the present study allow us to construct simpler and less computationally intensive models that can closely approximate the dynamics of heterogeneous cell populations.
| S-EPMC3181209 | biostudies-literature
Project description:Mechanical stimulation has been shown to reduce apnea of prematurity (AOP), a major concern in preterm infants. Previous work suggested that the underlying mechanism is stochastic resonance, amplification of a subthreshold signal by stochastic stimulation. We hypothesized that the mechanism behind the reduction of apnea length may not be a solely stochastic phenomenon, and suggest that a purely deterministic, non-random mechanical stimulation could be equally as effective. Mice and rats were anesthetized, tracheostomized, and mechanically ventilated to halt spontaneous breathing. Two miniature motors controlled by a microcontroller were attached around the abdomen. Ventilation was paused, stimulations were applied, and the time to the rodent's first spontaneous breath (T) was measured. Six spectrally different signals were compared to one another and the no-stimulation control in mice. The most successful deterministic stimulation (D) at reducing apnea was then compared to a pseudo-random noise (PRN) signal of comparable amplitude and frequency. CO2%, CO2 stabilization time (Ts), O2 saturation (SpO2%), and T were also measured. D significantly reduced T compared to no stimulation for medium and high amplitudes. PRN also reduced T, without a difference between D and PRN. Furthermore, both stimulations significantly reduced Ts with no significant differences between the respective stimulations. However, there was no effect of D or PRN on SpO2%. The lack of differences between D and PRN led to an additional series of experiment comparing the same D to a band-limited white noise (WN) signal in young rats. Both D and WN were shown to significantly reduce T, with D showing statistical superiority in reduction of apnea. We further speculate that both deterministic and stochastic mechanical stimulations induce some form of mechanotransduction which is responsible for their efficacy, and our findings suggest that mechanical stimulation may be effective in treating AOP.Supplementary informationThe online version contains supplementary material available at 10.1007/s13534-021-00203-x.
| S-EPMC8409477 | biostudies-literature
Project description:Hybrid deterministic-stochastic methods provide an efficient alternative to a fully stochastic treatment of models which include components with disparate levels of stochasticity. However, general-purpose hybrid solvers for spatially resolved simulations of reaction-diffusion systems are not widely available. Here we describe fundamentals of a general-purpose spatial hybrid method. The method generates realizations of a spatially inhomogeneous hybrid system by appropriately integrating capabilities of a deterministic partial differential equation solver with a popular particle-based stochastic simulator, Smoldyn. Rigorous validation of the algorithm is detailed, using a simple model of calcium 'sparks' as a testbed. The solver is then applied to a deterministic-stochastic model of spontaneous emergence of cell polarity. The approach is general enough to be implemented within biologist-friendly software frameworks such as Virtual Cell.
| S-EPMC5154471 | biostudies-other
Project description:BackgroundStochastic and asymptotic methods are powerful tools in developing multiscale systems biology models; however, little has been done in this context to compare the efficacy of these methods. The majority of current systems biology modelling research, including that of auxin transport, uses numerical simulations to study the behaviour of large systems of deterministic ordinary differential equations, with little consideration of alternative modelling frameworks.ResultsIn this case study, we solve an auxin-transport model using analytical methods, deterministic numerical simulations and stochastic numerical simulations. Although the three approaches in general predict the same behaviour, the approaches provide different information that we use to gain distinct insights into the modelled biological system. We show in particular that the analytical approach readily provides straightforward mathematical expressions for the concentrations and transport speeds, while the stochastic simulations naturally provide information on the variability of the system.ConclusionsOur study provides a constructive comparison which highlights the advantages and disadvantages of each of the considered modelling approaches. This will prove helpful to researchers when weighing up which modelling approach to select. In addition, the paper goes some way to bridging the gap between these approaches, which in the future we hope will lead to integrative hybrid models.
| S-EPMC2873313 | biostudies-literature
Project description:We present deterministic and stochastic programming models for the workover rig problem, one of the most challenging problems in the oil industry. In the deterministic approach, an integer linear programming model is used to determine the rig fleet size and schedule needed to service wells while maximizing oil production and minimizing rig usage cost. The stochastic approach is an extension of the deterministic method and relies on a two-stage stochastic programming model to define the optimal rig fleet size considering uncertainty in the intervention time. In this approach, different scenario-generation methods are compared. Several experiments were performed using instances based on real-world problems. The results suggest that the proposed methodology can be used to solve large instances and produces quality solutions in computationally reasonable times.
| S-EPMC6156096 | biostudies-literature
Project description:Detection of different extracellular stimuli leading to functionally distinct outcomes is ubiquitous in cell biology, and is often mediated by differential regulation of positive and negative feedback loops that are a part of the signaling network. In some instances, these cellular responses are stimulated by small numbers of molecules, and so stochastic effects could be important. Therefore, we studied the influence of stochastic fluctuations on a simple signaling model with dueling positive and negative feedback loops. The class of models we have studied is characterized by single deterministic steady states for all parameter values, but the stochastic response is bimodal; a behavior that is distinctly different from models studied in the context of gene regulation. For example, when positive and negative regulation is roughly balanced, a unique deterministic steady state with an intermediate value for the amount of a downstream signaling product is found. However, for small numbers of signaling molecules, stochastic effects result in a bimodal distribution for this quantity, with neither mode corresponding to the deterministic solution; i.e., cells are in "on" or "off" states, not in some intermediate state. For a large number of molecules, the stochastic solution converges to the mean-field result. When fluctuations are important, we find that signal output scales with control parameters "anomalously" compared with mean-field predictions. The necessary and sufficient conditions for the phenomenon we report are quite common. So, our findings are expected to be of broad relevance, and suggest that stochastic effects can enable binary cellular decisions.
| S-EPMC2141890 | biostudies-other
Project description:Estimation of epidemiological and population parameters from molecular sequence data has become central to the understanding of infectious disease dynamics. Various models have been proposed to infer details of the dynamics that describe epidemic progression. These include inference approaches derived from Kingman's coalescent theory. Here, we use recently described coalescent theory for epidemic dynamics to develop stochastic and deterministic coalescent susceptible-infected-removed (SIR) tree priors. We implement these in a Bayesian phylogenetic inference framework to permit joint estimation of SIR epidemic parameters and the sample genealogy. We assess the performance of the two coalescent models and also juxtapose results obtained with a recently published birth-death-sampling model for epidemic inference. Comparisons are made by analyzing sets of genealogies simulated under precisely known epidemiological parameters. Additionally, we analyze influenza A (H1N1) sequence data sampled in the Canterbury region of New Zealand and HIV-1 sequence data obtained from known United Kingdom infection clusters. We show that both coalescent SIR models are effective at estimating epidemiological parameters from data with large fundamental reproductive number [Formula: see text] and large population size [Formula: see text]. Furthermore, we find that the stochastic variant generally outperforms its deterministic counterpart in terms of error, bias, and highest posterior density coverage, particularly for smaller [Formula: see text] and [Formula: see text]. However, each of these inference models is shown to have undesirable properties in certain circumstances, especially for epidemic outbreaks with [Formula: see text] close to one or with small effective susceptible populations.
| S-EPMC4317665 | biostudies-literature