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:Calcium (Ca2+) is a second messenger assumed to control changes in synaptic strength in the form of both long-term depression and long-term potentiation at Purkinje cell dendritic spine synapses via inositol trisphosphate (IP3)-induced Ca2+ release. These Ca2+ transients happen in response to stimuli from parallel fibers (PFs) from granule cells and climbing fibers (CFs) from the inferior olivary nucleus. These events occur at low numbers of free Ca2+, requiring stochastic single-particle methods when modeling them. We use the stochastic particle simulation program MCell to simulate Ca2+ transients within a three-dimensional Purkinje cell dendritic spine. The model spine includes the endoplasmic reticulum, several Ca2+ transporters, and endogenous buffer molecules. Our simulations successfully reproduce properties of Ca2+ transients in different dynamical situations. We test two different models of the IP3 receptor (IP3R). The model with nonlinear concentration response of binding of activating Ca2+ reproduces experimental results better than the model with linear response because of the filtering of noise. Our results also suggest that Ca2+-dependent inhibition of the IP3R needs to be slow to reproduce experimental results. Simulations suggest the experimentally observed optimal timing window of CF stimuli arises from the relative timing of CF influx of Ca2+ and IP3 production sensitizing IP3R for Ca2+-induced Ca2+ release. We also model ataxia, a loss of fine motor control assumed to be the result of malfunctioning information transmission at the granule to Purkinje cell synapse, resulting in a decrease or loss of Ca2+ transients. Finally, we propose possible ways of recovering Ca2+ transients under ataxia.
| S-EPMC8390834 | biostudies-literature
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:Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared with approximation-based methods (e.g., finite element method [FEM]), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimensions; its formulation produces an efficient algorithm while maintaining coordinate-free geometric descriptions. However, ambiguities in geometric definitions in the discrete context have led to a lack of consensus on which discrete mesh model is theoretically and numerically optimal; a bijective relationship between the terms contributing to both the energy and forces from the discrete and smooth geometric theories remains to be established. We address this and present an extensible framework, Mem3DG, for modeling 3D mechanochemical dynamics of membranes based on discrete differential geometry (DDG) on triangulated meshes. The formalism of DDG resolves the inconsistency and provides a unifying perspective on how to relate the smooth and discrete energy and forces. To demonstrate, Mem3DG is used to model a sequence of examples with increasing mechanochemical complexity: recovering classical shape transformations such as 1) biconcave disk, dumbbell, and unduloid; and 2) spherical bud on spherical, flat-patch membrane; investigating how the coupling of membrane mechanics with protein mobility jointly affects phase and shape transformation. As high-resolution 3D imaging of membrane ultrastructure becomes more readily available, we envision Mem3DG to be applied as an end-to-end tool to simulate realistic cell geometry under user-specified mechanochemical conditions.
| S-EPMC9495267 | biostudies-literature
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-literature
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:MotivationDenoising diffusion probabilistic models (DDPMs) have recently taken the field of generative modeling by storm, pioneering new state-of-the-art results in disciplines such as computer vision and computational biology for diverse tasks ranging from text-guided image generation to structure-guided protein design. Along this latter line of research, methods have recently been proposed for generating 3D molecules using equivariant graph neural networks (GNNs) within a DDPM framework. However, such methods are unable to learn important geometric and physical properties of 3D molecules during molecular graph generation, as they adopt molecule-agnostic and non-geometric GNNs as their 3D graph denoising networks, which negatively impacts their ability to effectively scale to datasets of large 3D molecules.ResultsIn this work, we address these gaps by introducing the Geometry-Complete Diffusion Model (GCDM) for 3D molecule generation, which outperforms existing 3D molecular diffusion models by significant margins across conditional and unconditional settings for the QM9 dataset as well as for the larger GEOM-Drugs dataset. Importantly, we demonstrate that the geometry-complete denoising process GCDM learns for 3D molecule generation allows the model to generate realistic and stable large molecules at the scale of GEOM-Drugs, whereas previous methods fail to do so with the features they learn. Additionally, we show that extensions of GCDM can not only effectively design 3D molecules for specific protein pockets but also that GCDM's geometric features can effectively be repurposed to directly optimize the geometry and chemical composition of existing 3D molecules for specific molecular properties, demonstrating new, real-world versatility of molecular diffusion models.AvailabilityOur source code and data are freely available on GitHub.
| S-EPMC9934735 | 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