Stochastic hierarchical systems: excitable dynamics.
Ontology highlight
ABSTRACT:
SUBMITTER: Leonhardt H
PROVIDER: S-EPMC2652551 | biostudies-other | 2008 Oct
REPOSITORIES: biostudies-other
Similar Datasets
Project description:This study investigates the contributions of network topology features to the dynamic behavior of hierarchically organized excitable networks. Representatives of different types of hierarchical networks as well as two biological neural networks are explored with a three-state model of node activation for systematically varying levels of random background network stimulation. The results demonstrate that two principal topological aspects of hierarchical networks, node centrality and network modularity, correlate with the network activity patterns at different levels of spontaneous network activation. The approach also shows that the dynamic behavior of the cerebral cortical systems network in the cat is dominated by the network's modular organization, while the activation behavior of the cellular neuronal network of Caenorhabditis elegans is strongly influenced by hub nodes. These findings indicate the interaction of multiple topological features and dynamic states in the function of complex biological networks.
| S-EPMC2542420 | biostudies-literature
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:Spatiotemporal pattern formation governs dynamics and functions in various biological systems. In the heart, excitable waves can form complex oscillatory and chaotic patterns even at an abnormally higher frequency than normal heart beats, which increase the risk of fatal heart conditions by inhibiting normal blood circulation. Previous studies suggested that line defects (nodal lines) play a critical role in stabilizing those undesirable patterns. However, it remains unknown if the line defects are static or dynamically changing structures in heart tissue. Through in vitro experiments of heart tissue observation, we reveal the spatiotemporal dynamics of line defects in rotating spiral waves. We combined a novel signaling over-sampling technique with a multi-dimensional Fourier analysis, showing that line defects can translate, merge, collapse and form stable singularities with even and odd parity while maintaining a stable oscillation of the spiral wave in the tissue. These findings provide insights into a broad class of complex periodic systems, with particular impact to the control and understanding of heart diseases.
| S-EPMC5552747 | biostudies-other
Project description:Complex feedback systems are ubiquitous in biology. Modeling such systems with mass action laws or master equations requires information rarely measured directly. Thus rates and reaction topologies are often treated as adjustable parameters. Here we present a general stochastic modeling method for small chemical and biochemical systems with emphasis on feedback systems. The method, Maximum Caliber (MaxCal), is more parsimonious than others in constructing dynamical models requiring fewer model assumptions and parameters to capture the effects of feedback. MaxCal is the dynamical analogue of Maximum Entropy. It uses average rate quantities and correlations obtained from short experimental trajectories to construct dynamical models. We illustrate the method on the bistable genetic toggle switch. To test our method, we generate synthetic data from an underlying stochastic model. MaxCal reliably infers the statistics of the stochastic bistability and other full dynamical distributions of the simulated data, without having to invoke complex reaction schemes. The method should be broadly applicable to other systems.
| S-EPMC3098004 | biostudies-literature
Project description:Tumors are often heterogeneous in which tumor cells of different phenotypes have distinct properties. For scientific and clinical interests, it is of fundamental importance to understand their properties and the dynamic variations among different phenotypes, specifically under radio- and/or chemo-therapy. Currently there are two controversial models describing tumor heterogeneity, the cancer stem cell (CSC) model and the stochastic model. To clarify the controversy, we measured probabilities of different division types and transitions of cells via in situ immunofluorescence. Based on the experiment data, we constructed a model that combines the CSC with the stochastic concepts, showing the existence of both distinctive CSC subpopulations and the stochastic transitions from NSCCs to CSCs. The results showed that the dynamic variations between CSCs and non-stem cancer cells (NSCCs) can be simulated with the model. Further studies also showed that the model can be used to describe the dynamics of the two subpopulations after radiation treatment. More importantly, analysis demonstrated that the experimental detectable equilibrium CSC proportion can be achieved only when the stochastic transitions from NSCCs to CSCs occur, indicating that tumor heterogeneity may exist in a model coordinating with both the CSC and the stochastic concepts. The mathematic model based on experimental parameters may contribute to a better understanding of the tumor heterogeneity, and provide references on the dynamics of CSC subpopulation during radiotherapy.
| S-EPMC3886990 | biostudies-literature
Project description:Current estimates of the HIV epidemic indicate a decrease in the incidence of the disease in the undiagnosed subpopulation over the past 10 years. However, a lack of access to care has not been considered when modeling the population. Populations at high risk for contracting HIV are twice as likely to lack access to reliable medical care. In this paper, we consider three contributors to the HIV population dynamics: at-risk population exhaustion, lack of access to care, and usage of anti-retroviral therapy (ART) by diagnosed individuals. An extant problem in the mathematical study of this system is deriving parameter estimates due to a portion of the population being unobserved. We approach this problem by looking at the proportional change in the infected subpopulations. We obtain conservative estimates for the proportional change of the infected subpopulations using hierarchical Bayesian statistics. The estimated proportional change is used to derive epidemic parameter estimates for a system of stochastic differential equations (SDEs). Model fit is quantified to determine the best parametric explanation for the observed dynamics in the infected subpopulations. Parameter estimates derived using these methods produce simulations that closely follow the dynamics observed in the data, as well as values that are generally in agreement with prior understanding of transmission and diagnosis rates. Simulations suggest that the undiagnosed population may be larger than currently estimated without significantly affecting the population dynamics.
| S-EPMC6059410 | biostudies-literature
Project description:We study the effects of network topology on the response of networks of coupled discrete excitable systems to an external stochastic stimulus. We extend recent results that characterize the response in terms of spectral properties of the adjacency matrix by allowing distributions in the transmission delays and in the number of refractory states and by developing a nonperturbative approximation to the steady state network response. We confirm our theoretical results with numerical simulations. We find that the steady state response amplitude is inversely proportional to the duration of refractoriness, which reduces the maximum attainable dynamic range. We also find that transmission delays alter the time required to reach steady state. Importantly, neither delays nor refractoriness impact the general prediction that criticality and maximum dynamic range occur when the largest eigenvalue of the adjacency matrix is unity.
| S-EPMC3183795 | biostudies-other
Project description:Amoeboid cell migration is characterized by frequent changes of the direction of motion and resembles a persistent random walk on long time scales. Although it is well known that cell migration is typically driven by the actin cytoskeleton, the cause of this migratory behavior remains poorly understood. We analyze the spontaneous dynamics of actin assembly due to nucleation promoting factors, where actin filaments lead to an inactivation of these factors. We show that this system exhibits excitable dynamics and can spontaneously generate waves, which we analyze in detail. By using a phase-field approach, we show that these waves can generate cellular random walks. We explore how the characteristics of these persistent random walks depend on the parameters governing the actin-nucleator dynamics. In particular, we find that the effective diffusion constant and the persistence time depend strongly on the speed of filament assembly and the rate of nucleator inactivation. Our findings point to a deterministic origin of the random walk behavior and suggest that cells could adapt their migration pattern by modifying the pool of available actin.
| S-EPMC7850500 | biostudies-literature
Project description:Viewing of ambiguous stimuli can lead to bistable perception alternating between the possible percepts. During continuous presentation of ambiguous stimuli, percept changes occur as single events, whereas during intermittent presentation of ambiguous stimuli, percept changes occur at more or less regular intervals either as single events or bursts. Response patterns can be highly variable and have been reported to show systematic differences between patients with schizophrenia and healthy controls. Existing models of bistable perception often use detailed assumptions and large parameter sets which make parameter estimation challenging. Here we propose a parsimonious stochastic model that provides a link between empirical data analysis of the observed response patterns and detailed models of underlying neuronal processes. Firstly, we use a Hidden Markov Model (HMM) for the times between percept changes, which assumes one single state in continuous presentation and a stable and an unstable state in intermittent presentation. The HMM captures the observed differences between patients with schizophrenia and healthy controls, but remains descriptive. Therefore, we secondly propose a hierarchical Brownian model (HBM), which produces similar response patterns but also provides a relation to potential underlying mechanisms. The main idea is that neuronal activity is described as an activity difference between two competing neuronal populations reflected in Brownian motions with drift. This differential activity generates switching between the two conflicting percepts and between stable and unstable states with similar mechanisms on different neuronal levels. With only a small number of parameters, the HBM can be fitted closely to a high variety of response patterns and captures group differences between healthy controls and patients with schizophrenia. At the same time, it provides a link to mechanistic models of bistable perception, linking the group differences to potential underlying mechanisms.
| S-EPMC5714404 | biostudies-literature
Project description:Many animal cells initiate crawling by protruding lamellipodia, consisting of a dense network of actin filaments, at their leading edge. We imaged XTC cells that exhibit flat lamellipodia on poly-L-lysine-coated coverslips. Using active contours, we tracked the leading edge and measured the total amount of F-actin by summing the pixel intensities within a 5-?m band. We observed protrusion and retraction with period 130-200 s and local wavelike features. Positive (negative) velocities correlated with minimum (maximum) integrated actin concentration. Approximately constant retrograde flow indicated that protrusions and retractions were driven by fluctuations of the actin polymerization rate. We present a model of these actin dynamics as an excitable system in which a diffusive, autocatalytic activator causes actin polymerization; F-actin accumulation in turn inhibits further activator accumulation. Simulations of the model reproduced the pattern of actin polymerization seen in experiments. To explore the model's assumption of an autocatalytic activation mechanism, we imaged cells expressing markers for both F-actin and the p21 subunit of the Arp2/3 complex. We found that integrated Arp2/3-complex concentrations spike several seconds before spikes of F-actin concentration. This suggests that the Arp2/3 complex participates in an activation mechanism that includes additional diffuse components. Response of cells to stimulation by fetal calf serum could be reproduced by the model, further supporting the proposed dynamical picture.
| S-EPMC3318140 | biostudies-literature