Stochastic modeling of calcium in 3D geometry.
Ontology highlight
ABSTRACT:
SUBMITTER: Mazel T
PROVIDER: S-EPMC2996128 | biostudies-other | 2009 Mar
REPOSITORIES: biostudies-other
Publications
Stochastic modeling of calcium in 3D geometry.
Mazel Tomás T Raymond Rebecca R Raymond-Stintz Mary M Jett Stephen S Wilson Bridget S BS
Biophysical journal 20090301 5
Release of inflammatory mediators by mast cells in type 1 immediate-hypersensitivity allergic reactions relies on antigen-dependent increases in cytosolic calcium. Here, we used a series of electron microscopy images to build a 3D reconstruction representing a slice through a rat tumor mast cell, which then served as a basis for stochastic modeling of inositol-trisphosphate-mediated calcium responses. The stochastic approach was verified by reaction-diffusion modeling within the same geometry. L ...[more]
Similar Datasets
Project description:Deterministic models of biochemical processes at the subcellular level might become inadequate when a cascade of chemical reactions is induced by a few molecules. Inherent randomness of such phenomena calls for the use of stochastic simulations. However, being computationally intensive, such simulations become infeasible for large and complex reaction networks. To improve their computational efficiency in handling these networks, we present a hybrid approach, in which slow reactions and fluxes are handled through exact stochastic simulation and their fast counterparts are treated partially deterministically through chemical Langevin equation. The classification of reactions as fast or slow is accompanied by the assumption that in the time-scale of fast reactions, slow reactions do not occur and hence do not affect the probability of the state. Our new approach also handles reactions with complex rate expressions such as Michaelis-Menten kinetics. Fluxes which cannot be modeled explicitly through reactions, such as flux of Ca(2+) from endoplasmic reticulum to the cytosol through inositol 1,4,5-trisphosphate receptor channels, are handled deterministically. The proposed hybrid algorithm is used to model the regulation of the dynamics of cytosolic calcium ions in mouse macrophage RAW 264.7 cells. At relatively large number of molecules, the response characteristics obtained with the stochastic and deterministic simulations coincide, which validates our approach in the limit of large numbers. At low doses, the response characteristics of some key chemical species, such as levels of cytosolic calcium, predicted with stochastic simulations, differ quantitatively from their deterministic counterparts. These observations are ubiquitous throughout dose response, sensitivity, and gene-knockdown response analyses. While the relative differences between the peak-heights of the cytosolic [Ca(2+)] time-courses obtained from stochastic (mean of 16 realizations) and deterministic simulations are merely 1%-4% for most perturbations, it is specially sensitive to levels of G(βγ) (relative difference as large as 90% at very low G(βγ)).
| S-EPMC2998048 | biostudies-other
Project description:Gaussian random fields pervade all areas of science. However, it is often the departures from Gaussianity that carry the crucial signature of the nonlinear mechanisms at the heart of diverse phenomena, ranging from structure formation in condensed matter and cosmology to biomedical imaging. The standard test of non-Gaussianity is to measure higher-order correlation functions. In the present work, we take a different route. We show how geometric and topological properties of Gaussian fields, such as the statistics of extrema, are modified by the presence of a non-Gaussian perturbation. The resulting discrepancies give an independent way to detect and quantify non-Gaussianities. In our treatment, we consider both local and nonlocal mechanisms that generate non-Gaussian fields, both statically and dynamically through nonlinear diffusion.
| S-EPMC3523821 | biostudies-other
Project description:Phosphorylation reactions, driven by competing kinases and phosphatases, are central elements of cellular signal transduction. We reconstituted a native eukaryotic lipid kinase-phosphatase reaction that drives the interconversion of phosphatidylinositol-4-phosphate [PI(4)P] and phosphatidylinositol-4,5-phosphate [PI(4,5)P2] on membrane surfaces. This system exhibited bistability and formed spatial composition patterns on supported membranes. In smaller confined regions of membrane, rapid diffusion ensures the system remains spatially homogeneous, but the final outcome-a predominantly PI(4)P or PI(4,5)P2 membrane composition-was governed by the size of the reaction environment. In larger confined regions, interplay between the reactions, diffusion, and confinement created a variety of differentially patterned states, including polarization. Experiments and kinetic modeling reveal how these geometric confinement effects arise from a mechanism based on stochastic fluctuations in the copy number of membrane-bound kinases and phosphatases. The underlying requirements for such behavior are unexpectedly simple and likely to occur in natural biological signaling systems.
| S-EPMC6660746 | biostudies-literature
Project description:Because multivariate autoregressive models have failed to adequately account for the complexity of neural signals, researchers have predominantly relied on non-parametric methods when studying the relations between brain and behavior. Using medial temporal lobe (MTL) recordings from 96 neurosurgical patients, we show that time series models with volatility described by a multivariate stochastic latent-variable process and lagged interactions between signals in different brain regions provide new insights into the dynamics of brain function. The implied volatility inferred from our process positively correlates with high-frequency spectral activity, a signal that correlates with neuronal activity. We show that volatility features derived from our model can reliably decode memory states, and that this classifier performs as well as those using spectral features. Using the directional connections between brain regions during complex cognitive process provided by the model, we uncovered perirhinal-hippocampal desynchronization in the MTL regions that is associated with successful memory encoding.
| S-EPMC6697415 | biostudies-literature
Project description:Pose estimation of objects in real scenes is critically important for biological and machine visual systems, but little is known of how humans infer 3D poses from 2D retinal images. We show unexpectedly remarkable agreement in the 3D poses different observers estimate from pictures. We further show that all observers apply the same inferential rule from all viewpoints, utilizing the geometrically derived back-transform from retinal images to actual 3D scenes. Pose estimations are altered by a fronto-parallel bias, and by image distortions that appear to tilt the ground plane. We used pictures of single sticks or pairs of joined sticks taken from different camera angles. Observers viewed these from five directions, and matched the perceived pose of each stick by rotating an arrow on a horizontal touchscreen. The projection of each 3D stick to the 2D picture, and then onto the retina, is described by an invertible trigonometric expression. The inverted expression yields the back-projection for each object pose, camera elevation, and observer viewpoint. We show that a model that uses the back-projection, modulated by just two free parameters, explains 560 pose estimates per observer. By considering changes in retinal image orientations due to position and elevation of limbs, the model also explains perceived limb poses in a complex scene of two bodies lying on the ground. The inferential rules simply explain both perceptual invariance and dramatic distortions in poses of real and pictured objects, and show the benefits of incorporating projective geometry of light into mental inferences about 3D scenes.
| S-EPMC6065005 | biostudies-literature
Project description:It is now well-established that the host's adaptive immune system plays an important role in identifying and eliminating cancer cells in much the same way that intracellular pathogens are cleared during an adaptive immune response to infection. From a therapeutic standpoint, the adaptive immune system is unique in that it can co-evolve alongside a developing tumor. Tumor acquisition of immune evasive phenotypes, such as class-I MHC down-regulation, remains a major limitation of successful T-cell immunotherapy. Here, we consider a population dynamical model coupling tumor and adaptive immune compartments in order to study the dynamics and survival of an evolving threat when faced with adaptive immune pressure. We demonstrate that predicted optimal growth strategies depend on whether or not the threat may acquire an immune-evasive phenotype as well as the mode of immune detection. We parameterize adaptive immune functioning by T-cell turnover and repertoire diversity and predict that decreases in the latter quantity which occur in advanced age may substantially affect the ability to recognize, and therefore control, an immune evasive threat like cancer. This framework recapitulates general features of age-dependent AML incidence, thereby providing a probable association between cancer frequency and adaptive immune functioning. Lastly, we quantify therapeutic efficacy of adjuvant immunotherapeutic strategies, and predict their benefits and limitations with regard to handling immune evasion. Our model generates survival behavior consistent with known growth-dependent characteristics, and serves as a first attempt at modeling stochastic cancer evolution alongside an adaptive immune compartment.
| S-EPMC6225903 | biostudies-literature
Project description:We present results of a study of a simple, stochastic, agent-based model of influenza A infection, simulating its dynamics over the course of one flu season. Building on an early work of Bartlett, we define a model with a limited number of parameters and rates that have clear epidemiological interpretation and can be constrained by data. We demonstrate the occurrence of recurrent behavior in the infected number [more than one peak in a season], which is observed in data, in our simulations for populations consisting of cohorts with strong intra- and weak inter-cohort transmissibility. We examine the dependence of the results on epidemiological and population characteristics by investigating their dependence on a range of parameter values. Finally, we study infection with two strains of influenza, inspired by observations, and show a counter-intuitive result for the effect of inoculation against the strain that leads to the first wave of infection.
| S-EPMC7173783 | biostudies-literature
Project description:High risk forms of human papillomaviruses (HPVs) promote cancerous lesions and are implicated in almost all cervical cancer. Of particular relevance to cancer progression is regulation of the early promoter that controls gene expression in the initial phases of infection and can eventually lead to pre-cancer progression. Our goal was to develop a stochastic model to investigate the control mechanisms that regulate gene expression from the HPV early promoter. Our model integrates modules that account for transcriptional, post-transcriptional, translational and post-translational regulation of E1 and E2 early genes to form a functioning gene regulatory network. Each module consists of a set of biochemical steps whose stochastic evolution is governed by a chemical Master Equation and can be simulated using the Gillespie algorithm. To investigate the role of noise in gene expression, we compared our stochastic simulations with solutions to ordinary differential equations for the mean behavior of the system that are valid under the conditions of large molecular abundances and quasi-equilibrium for fast reactions. The model produced results consistent with known HPV biology. Our simulation results suggest that stochasticity plays a pivotal role in determining the dynamics of HPV gene expression. In particular, the combination of positive and negative feedback regulation generates stochastic bursts of gene expression. Analysis of the model reveals that regulation at the promoter affects burst amplitude and frequency, whereas splicing is more specialized to regulate burst frequency. Our results also suggest that splicing enhancers are a significant source of stochasticity in pre-mRNA abundance and that the number of viruses infecting the host cell represents a third important source of stochasticity in gene expression.
| S-EPMC6937396 | biostudies-literature
Project description:A methodology was developed for the prediction of the electrical properties of carbon nanotube-polymer nanocomposites via Monte Carlo computational simulations. A two-dimensional microstructure that takes into account waviness, fiber length and diameter distributions is used as a representative volume element. Fiber interactions in the microstructure are identified and then modeled as an equivalent electrical circuit, assuming one-third metallic and two-thirds semiconductor nanotubes. Tunneling paths in the microstructure are also modeled as electrical resistors, and crossing fibers are accounted for by assuming a contact resistance associated with them. The equivalent resistor network is then converted into a set of linear equations using nodal voltage analysis, which is then solved by means of the Gauss-Jordan elimination method. Nodal voltages are obtained for the microstructure, from which the percolation probability, equivalent resistance and conductivity are calculated. Percolation probability curves and electrical conductivity values are compared to those found in the literature.
| S-EPMC5455403 | biostudies-other
Project description:In cardiac ventricular myocytes, events crucial to excitation-contraction coupling take place in spatially restricted microdomains known as dyads. The movement and dynamics of calcium (Ca2+) ions in the dyad have often been described by assigning continuously valued Ca2+ concentrations to one or more dyadic compartments. However, even at its peak, the estimated number of free Ca2+ ions present in a single dyad is small (approximately 10-100 ions). This in turn suggests that modeling dyadic calcium dynamics using laws of mass action may be inappropriate. In this study, we develop a model of stochastic molecular signaling between L-type Ca2+ channels (LCCs) and ryanodine receptors (RyR2s) that describes: a), known features of dyad geometry, including the space-filling properties of key dyadic proteins; and b), movement of individual Ca2+ ions within the dyad, as driven by electrodiffusion. The model enables investigation of how local Ca2+ signaling is influenced by dyad structure, including the configuration of key proteins within the dyad, the location of Ca2+ binding sites, and membrane surface charges. Using this model, we demonstrate that LCC-RyR2 signaling is influenced by both the stochastic dynamics of Ca2+ ions in the dyad as well as the shape and relative positioning of dyad proteins. Results suggest the hypothesis that the relative placement and shape of the RyR2 proteins helps to "funnel" Ca2+ ions to RyR2 binding sites, thus increasing excitation-contraction coupling gain.
| S-EPMC1853149 | biostudies-literature