Pools versus Queues: The Variable Dynamics of Stochastic "Steady States".
Ontology highlight
ABSTRACT:
SUBMITTER: Lofgren ET
PROVIDER: S-EPMC4474550 | biostudies-other | 2015
REPOSITORIES: biostudies-other
Publications
Pools versus Queues: The Variable Dynamics of Stochastic "Steady States".
PloS one 20150619 6
Mathematical models in ecology and epidemiology often consider populations "at equilibrium", where in-flows, such as births, equal out-flows, such as death. For stochastic models, what is meant by equilibrium is less clear--should the population size be fixed or growing and shrinking with equal probability? Two different mechanisms to implement a stochastic steady state are considered. Under these mechanisms, both a predator-prey model and an epidemic model have vastly different outcomes, includ ...[more]
Similar Datasets
Project description:Many biochemical networks have complex multidimensional dynamics and there is a long history of methods that have been used for dimensionality reduction for such reaction networks. Usually a deterministic mass action approach is used; however, in small volumes, there are significant fluctuations from the mean which the mass action approach cannot capture. In such cases stochastic simulation methods should be used. In this paper, we evaluate the applicability of one such dimensionality reduction method, the quasi-steady state approximation (QSSA) [L. Menten and M. Michaelis, "Die kinetik der invertinwirkung," Biochem. Z 49, 333369 (1913)] for dimensionality reduction in case of stochastic dynamics. First, the applicability of QSSA approach is evaluated for a canonical system of enzyme reactions. Application of QSSA to such a reaction system in a deterministic setting leads to Michaelis-Menten reduced kinetics which can be used to derive the equilibrium concentrations of the reaction species. In the case of stochastic simulations, however, the steady state is characterized by fluctuations around the mean equilibrium concentration. Our analysis shows that a QSSA based approach for dimensionality reduction captures well the mean of the distribution as obtained from a full dimensional simulation but fails to accurately capture the distribution around that mean. Moreover, the QSSA approximation is not unique. We have then extended the analysis to a simple bistable biochemical network model proposed to account for the stability of synaptic efficacies; the substrate of learning and memory [J. E. Lisman, "A mechanism of memory storage insensitive to molecular turnover: A bistable autophosphorylating kinase," Proc. Natl. Acad. Sci. U.S.A. 82, 3055-3057 (1985)]. Our analysis shows that a QSSA based dimensionality reduction method results in errors as big as two orders of magnitude in predicting the residence times in the two stable states.
| S-EPMC3416873 | biostudies-other
Project description:UnlabelledBackgroundIn order to replicate within their cellular host, many viruses have developed self-assembly strategies for their capsids which are sufficiently robust as to be reconstituted in vitro. Mathematical models for virus self-assembly usually assume that the bonds leading to cluster formation have constant reactivity over the time course of assembly (direct assembly). In some cases, however, binding sites between the capsomers have been reported to be activated during the self-assembly process (hierarchical assembly).ResultsIn order to study possible advantages of such hierarchical schemes for icosahedral virus capsid assembly, we use Brownian dynamics simulations of a patchy particle model that allows us to switch binding sites on and off during assembly. For T1 viruses, we implement a hierarchical assembly scheme where inter-capsomer bonds become active only if a complete pentamer has been assembled. We find direct assembly to be favorable for reversible bonds allowing for repeated structural reorganizations, while hierarchical assembly is favorable for strong bonds with small dissociation rate, as this situation is less prone to kinetic trapping. However, at the same time it is more vulnerable to monomer starvation during the final phase. Increasing the number of initial monomers does have only a weak effect on these general features. The differences between the two assembly schemes become more pronounced for more complex virus geometries, as shown here for T3 viruses, which assemble through homogeneous pentamers and heterogeneous hexamers in the hierarchical scheme. In order to complement the simulations for this more complicated case, we introduce a master equation approach that agrees well with the simulation results.ConclusionsOur analysis shows for which molecular parameters hierarchical assembly schemes can outperform direct ones and suggests that viruses with high bond stability might prefer hierarchical assembly schemes. These insights increase our physical understanding of an essential biological process, with many interesting potential applications in medicine and materials science.
| S-EPMC3563543 | biostudies-literature
Project description:Gene expression has been investigated in relation with growth rate in the yeast Saccharomyces cerevisiae, following different experimental strategies. The expression of some specific gene functional categories increases or decreases with growth rate. Our recently published results have unveiled that these changes in mRNA concentration with growth depend on the relative alteration of mRNA synthesis and decay, and that, in addition to this gene-specific transcriptomic signature of growth, global mRNA turnover increases with growth rate. We discuss here these results in relation with other previous and concurrent publications, and we add new evidence which indicates that growth rate controls mRNA turnover even under non-steady-state conditions.
| S-EPMC5207384 | biostudies-literature
Project description:Even simple active systems can show a plethora of intriguing phenomena and often we find complexity where we would have expected simplicity. One striking example is the occurrence of a quiescent or absorbing state with frozen fluctuations that at first sight seems to be impossible for active matter driven by the incessant input of energy. While such states were reported for externally driven systems through macroscopic shear or agitation, the investigation of frozen active states in inherently active systems like cytoskeletal suspensions or active gels is still at large. Using high-density motility assay experiments, we demonstrate that frozen steady states can arise in active systems if active transport is coupled to growth processes. The experiments are complemented by agent-based simulations which identify the coupling between self-organization, growth, and mechanical properties to be responsible for the pattern formation process.
| S-EPMC3228451 | biostudies-literature
Project description:Plasmin (PLS) and urokinase-type plasminogen activator (UPA) are ubiquitous proteases that regulate the extracellular environment. Although they are secreted in inactive forms, they can activate each other through proteolytic cleavage. This mutual interplay creates the potential for complex dynamics, which we investigated using mathematical modeling and in vitro experiments. We constructed ordinary differential equations to model the conversion of precursor plasminogen into active PLS, and precursor urokinase (scUPA) into active urokinase (tcUPA). Although neither PLS nor UPA exhibits allosteric cooperativity, modeling showed that cooperativity occurred at the system level because of substrate competition. Computational simulations and bifurcation analysis predicted that the system would be bistable over a range of parameters for cooperativity and positive feedback. Cell-free experiments with recombinant proteins tested key predictions of the model. PLS activation in response to scUPA stimulus was found to be cooperative in vitro. Finally, bistability was demonstrated in vitro by the presence of two significantly different steady-state levels of PLS activation for the same levels of stimulus. We conclude that ultrasensitive, bistable activation of UPA-PLS is possible in the presence of substrate competition. An ultrasensitive threshold for activation of PLS and UPA would have ramifications for normal and disease processes, including angiogenesis, metastasis, wound healing, and fibrosis.
| S-EPMC3192967 | biostudies-other
Project description:In biochemical networks, reactions often occur on disparate timescales and can be characterized as either fast or slow. The quasi-steady-state approximation (QSSA) utilizes timescale separation to project models of biochemical networks onto lower-dimensional slow manifolds. As a result, fast elementary reactions are not modeled explicitly, and their effect is captured by nonelementary reaction-rate functions (e.g., Hill functions). The accuracy of the QSSA applied to deterministic systems depends on how well timescales are separated. Recently, it has been proposed to use the nonelementary rate functions obtained via the deterministic QSSA to define propensity functions in stochastic simulations of biochemical networks. In this approach, termed the stochastic QSSA, fast reactions that are part of nonelementary reactions are not simulated, greatly reducing computation time. However, it is unclear when the stochastic QSSA provides an accurate approximation of the original stochastic simulation. We show that, unlike the deterministic QSSA, the validity of the stochastic QSSA does not follow from timescale separation alone, but also depends on the sensitivity of the nonelementary reaction rate functions to changes in the slow species. The stochastic QSSA becomes more accurate when this sensitivity is small. Different types of QSSAs result in nonelementary functions with different sensitivities, and the total QSSA results in less sensitive functions than the standard or the prefactor QSSA. We prove that, as a result, the stochastic QSSA becomes more accurate when nonelementary reaction functions are obtained using the total QSSA. Our work provides an apparently novel condition for the validity of the QSSA in stochastic simulations of biochemical reaction networks with disparate timescales.
| S-EPMC4129492 | biostudies-literature
Project description:Neuronal growth is an extremely complex yet reliable process that is directed by a dynamic lamellipodial structure at the tip of every growing neurite, called the growth cone. Lamellipodial edge fluctuations are controlled by the interplay between actin polymerization pushing the edge forward and molecular motor driven retrograde actin flow retracting the actin network. The leading edge switches randomly between extension and retraction processes. We identify switching of "on/off" states in actin polymerization as the main determinant of lamellipodial advancement. Our analysis of motility statistics allows for a prediction of growth direction. This was used in simulations explaining the amazing signal detection capabilities of neuronal growth by the experimentally found biased stochastic processes. Our measurements show that the intensity of stochastic fluctuations depend on changes in the underlying active intracellular processes and we find a power law eta = a*x(alpha) with exponent alpha = 2.63 +/- 0.12 between noise intensity eta and growth cone activity x, defined as the sum of protrusion and retraction velocity. Differences in the lamellipodial dynamics between primary neurons and a neuronal cell line further suggests that active processes tune the observed stochastic fluctuations. This hints at a possible role of noise intensity in determining signal detection sensitivity.
| S-EPMC2712033 | biostudies-literature
Project description:Gait is a large component and indicator of health. Many factors affect gait including age, disease, and even mood disorders. Few studies have looked at the influence of emotional states on gait. This study aimed to investigate the influence of emotional states on walking performance to understand whether an emotional state may be an important factor to consider when evaluating gait. Thirty-six young adults were recruited (23F, 13M) and performed a neutral baseline condition of walking which included six passes of walking across an 8m walkway (a total of 48m of walking). Participants then completed 6 pseudo-randomized emotional state induction conditions while immersive 360-degree videos were used to induce the following emotional conditions: happiness, excitement, sadness, fear, and anger. Participants viewed the emotion elicitation videos using a virtual reality head-mounted display (HMD), then rated their emotional state using self-assessment manikins and walked (without the HMD) over a pressure sensor walkway. One-way repeated measures ANOVA and pairwise comparisons were used to examine differences in gait parameters across the emotional conditions. Participants walked with significantly reduced step length and speed during the sadness condition compared to the other emotional conditions and the neutral condition. Furthermore, participants adjusted the timing of their walking during the sadness condition and walked with significantly increased step, stance, and swing times compared to other emotional conditions, but not the neutral condition. Step time was significantly reduced during the conditions of excitement and fear compared to the neutral condition. Emotions may impact variety of gait parameters involving pace and rhythm, however have little influence on gait variability and postural control. These results indicate that perhaps the emotions of sadness and excitement should be taken into account as potential confounds for future gait analysis.
| S-EPMC10501668 | biostudies-literature
Project description:Living systems use fuel-driven supramolecular polymers such as actin to control important cell functions. Fuel molecules like ATP are used to control when and where such polymers should assemble and disassemble. The cell supplies fresh ATP to the cytosol and removes waste products to sustain steady states. Artificial fuel-driven polymers have been developed recently, but keeping them in sustained non-equilibrium steady states (NESS) has proven challenging. Here we show a supramolecular polymer that can be kept in NESS, inside a membrane reactor where ATP is added and waste removed continuously. Assembly and disassembly of our polymer is regulated by phosphorylation and dephosphorylation, respectively. Waste products lead to inhibition, causing the reaction cycle to stop. Inside the membrane reactor, however, waste can be removed leading to long-lived NESS conditions. We anticipate that our approach to obtain NESS can be applied to other stimuli-responsive materials to achieve more life-like behaviour.
| S-EPMC5481825 | biostudies-literature
Project description:BackgroundThe rapid and accurate identification of a minimal-size core set of representative microbial species plays an important role in the clustering of microbial community data and interpretation of clustering results. However, the huge dimensionality of microbial metagenomics datasets is a major challenge for the existing methods such as Dirichlet multinomial mixture (DMM) models. In the approach of the existing methods, the computational burden of identifying a small number of representative species from a large number of observed species remains a challenge.ResultsWe propose a novel approach to improve the performance of the widely used DMM approach by combining three ideas: (i) we propose an indicator variable to identify representative operational taxonomic units that substantially contribute to the differentiation among clusters; (ii) to address the computational burden of high-dimensional microbiome data, we propose a stochastic variational inference, which approximates the posterior distribution using a controllable distribution called variational distribution, and stochastic optimization algorithms for fast computation; and (iii) we extend the finite DMM model to an infinite case by considering Dirichlet process mixtures and estimating the number of clusters as a variational parameter. Using the proposed method, stochastic variational variable selection (SVVS), we analyzed the root microbiome data collected in our soybean field experiment, the human gut microbiome data from three published datasets of large-scale case-control studies and the healthy human microbiome data from the Human Microbiome Project.ConclusionsSVVS demonstrates a better performance and significantly faster computation than those of the existing methods in all cases of testing datasets. In particular, SVVS is the only method that can analyze massive high-dimensional microbial data with more than 50,000 microbial species and 1000 samples. Furthermore, a core set of representative microbial species is identified using SVVS that can improve the interpretability of Bayesian mixture models for a wide range of microbiome studies. Video Abstract.
| S-EPMC9789572 | biostudies-literature