Complex intracellular calcium oscillations. A theoretical exploration of possible mechanisms

Mathematical modeling of intracellular calcium in presence of receptor: a homeostatic model for endothelial cell

Calcium is a ubiquitous molecule and second messenger that regulates many cellular functions ranging from exocytosis to cell proliferation at different time scales. In the vasculature, a constant adenosine triphosphate (ATP) concentration is maintained because of ATP released by red blood cells (RBCs). These ATP molecules continuously react with purinergic receptors on the surface of endothelial cells (ECs). Consequently, a cascade of chemical reactions are triggered that result in a transient cytoplasmic calcium (Ca[Formula: see text]), followed by return to its basal concentration. The mathematical models proposed in the literature are able to reproduce the transient peak. However, the trailing concentration is always higher than the basal cytoplasmic Ca[Formula: see text] concentrations, and the Ca[Formula: see text] concentration in endoplasmic reticulum (ER) remains lower than its initial concentration. This means that the intracellular homeostasis is not recovered. We propose, herein, a minimal model of calcium kinetics. We find that the desensitization of EC surface receptors due to phosphorylation and recycling plays a vital role in maintaining calcium homeostasis in the presence of a constant stimulus (ATP). The model is able to capture several experimental observations such as refilling of Ca[Formula: see text] in the ER, variation of cytoplasmic Ca[Formula: see text] transient peak in ECs, the resting cytoplasmic Ca[Formula: see text] concentration, the effect of removing ATP from the plasma on Ca[Formula: see text] homeostasis, and the saturation of cytoplasmic Ca[Formula: see text] transient peak with increase in ATP concentration. Direct confrontation with several experimental results is conducted. This work paves the way for systematic studies on coupling between blood flow and chemical signaling, and should contribute to a better understanding of the relation between (patho)physiological conditions and Ca[Formula: see text] kinetics.

Mathematical Analysis of Membrane Transporters Dynamics: A Calcium Fluxes Case Study

A tight control of intracellular [Ca[Formula: see text]] is essential for the survival and normal function of cells. In this study we investigate key mechanistic steps by which calcium is regulated and calcium oscillations could occur using in silico modeling of membrane transporters. To do so we give a deterministic description of intracellular Ca[Formula: see text] dynamics using nonlinear dynamics in order to understand Ca[Formula: see text] signaling. We first present the ordinary differential equations (ODEs) system for cell calcium kinetics and make a preliminary work on Sobol indices. We then describe and analyze complex transporters action. Besides, we analyze the whole system. We finally perform numerical simulations and compare our results to real data.

Principles, mechanisms and functions of entrainment in biological oscillators

Entrainment is a phenomenon in which two oscillators interact with each other, typically through physical or chemical means, to synchronize their oscillations. This phenomenon occurs in biology to coordinate processes from the molecular to organismal scale. Biological oscillators can be entrained within a single cell, between cells or to an external input. Using six illustrative examples of entrainable biological oscillators, we discuss the distinctions between entrainment and synchrony and explore features that contribute to a system's propensity to entrain. Entrainment can either enhance or reduce the heterogeneity of oscillations within a cell population, and we provide examples and mechanisms of each case. Finally, we discuss the known functions of entrainment and discuss potential functions from an evolutionary perspective.

Deregulation of Ca2+-Signaling Systems in White Adipocytes, Manifested as the Loss of Rhythmic Activity, Underlies the Development of Multiple Hormonal Resistance at Obesity and Type 2 Diabetes

Various types of cells demonstrate ubiquitous rhythmicity registered as simple and complex Ca²⁺-oscillations, spikes, waves, and triggering phenomena mediated by G-protein and tyrosine kinase coupled receptors. Phospholipase C/IP₃-receptors (PLC/IP₃R) and endothelial NO-synthase/Ryanodine receptors (NOS/RyR)–dependent Ca²⁺ signaling systems, organized as multivariate positive feedback generators (PLC-G and NOS-G), underlie this rhythmicity. Loss of rhythmicity at obesity may indicate deregulation of these signaling systems. To issue the impact of cell size, receptors’ interplay, and obesity on the regulation of PLC-G and NOS-G, we applied fluorescent microscopy, immunochemical staining, and inhibitory analysis using cultured adipocytes of epididumal white adipose tissue of mice. Acetylcholine, norepinephrine, atrial natriuretic peptide, bradykinin, cholecystokinin, angiotensin II, and insulin evoked complex [Ca²⁺]_i responses in adipocytes, implicating NOS-G or PLC-G. At low sub-threshold concentrations, acetylcholine and norepinephrine or acetylcholine and peptide hormones (in paired combinations) recruited NOS-G, based on G proteins subunits interplay and signaling amplification. Rhythmicity was cell size- dependent and disappeared in hypertrophied cells filled with lipids. Contrary to control cells, adipocytes of obese hyperglycemic and hypertensive mice, growing on glucose, did not accumulate lipids and demonstrated hormonal resistance being non responsive to any hormone applied. Preincubation of preadipocytes with palmitoyl-L-carnitine (100 nM) provided accumulation of lipids, increased expression and clustering of IP₃R and RyR proteins, and partially restored hormonal sensitivity and rhythmicity (5–15% vs. 30–80% in control cells), while adipocytes of diabetic mice were not responsive at all. Here, we presented a detailed kinetic model of NOS-G and discussed its control. Collectively, we may suggest that universal mechanisms underlie loss of rhythmicity, Ca²⁺-signaling systems deregulation, and development of general hormonal resistance to obesity.

Modeling of Ca2+ transients initiated by GPCR agonists in mesenchymal stromal cells

The integrative study that included experimentation and mathematical modeling was carried out to analyze dynamic aspects of transient Ca²⁺ signaling induced by brief pulses of GPCR agonists in mesenchymal stromal cells from the human adipose tissue (AD-MSCs). The experimental findings argued for IP₃/Ca²⁺-regulated Ca²⁺ release via IP₃ receptors (IP₃Rs) as a key mechanism mediating agonist-dependent Ca²⁺ transients. The consistent signaling circuit was proposed to formalize coupling of agonist binding to Ca²⁺ mobilization for mathematical modeling. The model properly simulated the basic phenomenology of agonist transduction in AD-MSCs, which mostly produced single Ca²⁺ spikes upon brief stimulation. The spike-like responses were almost invariantly shaped at different agonist doses above a threshold, while response lag markedly decreased with stimulus strength. In AD-MSCs, agonists and IP₃ uncaging elicited similar Ca²⁺ transients but IP₃ pulses released Ca²⁺ without pronounced delay. This suggested that IP₃ production was rate-limiting in agonist transduction. In a subpopulation of AD-MSCs, brief agonist pulses elicited Ca²⁺ bursts crowned by damped oscillations. With properly adjusted parameters of IP₃R inhibition by cytosolic Ca²⁺, the model reproduced such oscillatory Ca²⁺ responses as well. GEM-GECO1 and R-CEPIA1er, the genetically encoded sensors of cytosolic and reticular Ca²⁺, respectively, were co-expressed in HEK-293 cells that also responded to agonists in an "all-or-nothing" manner. The experimentally observed Ca²⁺ signals triggered by ACh in both compartments were properly simulated with the suggested signaling circuit. Thus, the performed modeling of the transduction process provides sufficient theoretical basis for deeper interpretation of experimental findings on agonist-induced Ca²⁺ signaling in AD-MSCs.

Modeling and simulation of biological systems using SPICE language

The article deals with BB-SPICE (SPICE for Biochemical and Biological Systems), an extension of the famous Simulation Program with Integrated Circuit Emphasis (SPICE). BB-SPICE environment is composed of three modules: a new textual and compact description formalism for biological systems, a converter that handles this description and generates the SPICE netlist of the equivalent electronic circuit and NGSPICE which is an open-source SPICE simulator. In addition, the environment provides back and forth interfaces with SBML (System Biology Markup Language), a very common description language used in systems biology. BB-SPICE has been developed in order to bridge the gap between the simulation of biological systems on the one hand and electronics circuits on the other hand. Thus, it is suitable for applications at the interface between both domains, such as development of design tools for synthetic biology and for the virtual prototyping of biosensors and lab-on-chip. Simulation results obtained with BB-SPICE and COPASI (an open-source software used for the simulation of biochemical systems) have been compared on a benchmark of models commonly used in systems biology. Results are in accordance from a quantitative viewpoint but BB-SPICE outclasses COPASI by 1 to 3 orders of magnitude regarding the computation time. Moreover, as our software is based on NGSPICE, it could take profit of incoming updates such as the GPU implementation, of the coupling with powerful analysis and verification tools or of the integration in design automation tools (synthetic biology).

Discerning mechanistically rewired biological pathways by cumulative interaction heterogeneity statistics

Changes in response of a biological pathway could be a consequence of either pathway rewiring, changed input, or a combination of both. Most pathway analysis methods are not designed for mechanistic rewiring such as regulatory element variations. This limits our understanding of biological pathway evolution. Here we present a Q-method to discern whether changed pathway response is caused by mechanistic rewiring of pathways due to evolution. The main innovation is a cumulative pathway interaction heterogeneity statistic accounting for rewiring-specific effects on the rate of change of each molecular variable across conditions. The Q-method remarkably outperformed differential-correlation based approaches on data from diverse biological processes. Strikingly, it also worked well in differentiating rewired chaotic systems, whose dynamics are notoriously difficult to predict. Applying the Q-method on transcriptome data of four yeasts, we show that pathway interaction heterogeneity for known metabolic and signaling pathways is indeed a predictor of interspecies genetic rewiring due to unbalanced TATA box-containing genes among the yeasts. The demonstrated effectiveness of the Q-method paves the way to understanding network evolution at the resolution of functional biological pathways.

Estimating the biophysical properties of neurons with intracellular calcium dynamics

We investigate the dynamics of a conductance-based neuron model coupled to a model of intracellular calcium uptake and release by the endoplasmic reticulum. The intracellular calcium dynamics occur on a time scale that is orders of magnitude slower than voltage spiking behavior. Coupling these mechanisms sets the stage for the appearance of chaotic dynamics, which we observe within certain ranges of model parameter values. We then explore the question of whether one can, using observed voltage data alone, estimate the states and parameters of the voltage plus calcium (V+Ca) dynamics model. We find the answer is negative. Indeed, we show that voltage plus another observed quantity must be known to allow the estimation to be accurate. We show that observing both the voltage time course V(t) and the intracellular Ca time course will permit accurate estimation, and from the estimated model state, accurate prediction after observations are completed. This sets the stage for how one will be able to use a more detailed model of V+Ca dynamics in neuron activity in the analysis of experimental data on individual neurons as well as functional networks in which the nodes (neurons) have these biophysical properties.

Modelling the transition from simple to complex Ca²⁺ oscillations in pancreatic acinar cells

A mathematical model is proposed which systematically investigates complex calcium oscillations in pancreatic acinar cells. This model is based on calcium-induced calcium release via inositol trisphosphate receptors (IPR) and ryanodine receptors (RyR) and includes calcium modulation of inositol (1,4,5) trisphosphate (IP3) levels through feedback regulation of degradation and production. In our model, the apical and the basal regions are separated by a region containing mitochondria, which is capable of restricting Ca2+ responses to the apical region. We were able to reproduce the observed oscillatory patterns, from baseline spikes to sinusoidal oscillations. The model predicts that calcium-dependent production and degradation of IP3 is a key mechanism for complex calcium oscillations in pancreatic acinar cells. A partial bifurcation analysis is performed which explores the dynamic behaviour of the model in both apical and basal regions.

Cell-type-specific modelling of intracellular calcium signalling: a urothelial cell model

Calcium signalling plays a central role in regulating a wide variety of cell processes. A number of calcium signalling models exist in the literature that are capable of reproducing a variety of experimentally observed calcium transients. These models have been used to examine in more detail the mechanisms underlying calcium transients, but very rarely has a model been directly linked to a particular cell type and experimentally verified. It is important to show that this can be achieved within the general theoretical framework adopted by these models. Here, we develop a framework designed specifically for modelling cytosolic calcium transients in urothelial cells. Where possible, we draw upon existing calcium signalling models, integrating descriptions of components known to be important in this cell type from a number of studies in the literature. We then add descriptions of several additional pathways that play a specific role in urothelial cell signalling, including an explicit ionic influx term and an active pumping mechanism that drives the cytosolic calcium concentration to a target equilibrium. The resulting one-pool model of endoplasmic reticulum (ER)-dependent calcium signalling relates the cytosolic, extracellular and ER calcium concentrations and can generate a wide range of calcium transients, including spikes, bursts, oscillations and sustained elevations in the cytosolic calcium concentration. Using single-variate robustness and multivariate sensitivity analyses, we quantify how varying each of the parameters of the model leads to changes in key features of the calcium transient, such as initial peak amplitude and the frequency of bursting or spiking, and in the transitions between bursting- and plateau-dominated modes. We also show that, novel to our urothelial cell model, the ionic and purinergic P2Y pathways make distinct contributions to the calcium transient. We then validate the model using human bladder epithelial cells grown in monolayer cell culture and show that the model robustly captures the key features of the experimental data in a way that is not possible using more generic calcium models from the literature.

Membrane associated complexes in calcium dynamics modelling

Mitochondria not only govern energy production, but are also involved in crucial cellular signalling processes. They are one of the most important organelles determining the Ca(2+) regulatory pathway in the cell. Several mathematical models explaining these mechanisms were constructed, but only few of them describe interplay between calcium concentrations in endoplasmic reticulum (ER), cytoplasm and mitochondria. Experiments measuring calcium concentrations in mitochondria and ER suggested the existence of cytosolic microdomains with locally elevated calcium concentration in the nearest vicinity of the outer mitochondrial membrane. These intermediate physical connections between ER and mitochondria are called MAM (mitochondria-associated ER membrane) complexes. We propose a model with a direct calcium flow from ER to mitochondria, which may be justified by the existence of MAMs, and perform detailed numerical analysis of the effect of this flow on the type and shape of calcium oscillations. The model is partially based on the Marhl et al model. We have numerically found that the stable oscillations exist for a considerable set of parameter values. However, for some parameter sets the oscillations disappear and the trajectories of the model tend to a steady state with very high calcium level in mitochondria. This can be interpreted as an early step in an apoptotic pathway.

Acetylcholine promotes Ca2+ and NO-oscillations in adipocytes implicating Ca2+→NO→cGMP→cADP-ribose→Ca2+ positive feedback loop--modulatory effects of norepinephrine and atrial natriuretic peptide

PURPOSE

This study investigated possible mechanisms of autoregulation of Ca(2+) signalling pathways in adipocytes responsible for Ca(2+) and NO oscillations and switching phenomena promoted by acetylcholine (ACh), norepinephrine (NE) and atrial natriuretic peptide (ANP).

METHODS

Fluorescent microscopy was used to detect changes in Ca(2+) and NO in cultures of rodent white adipocytes. Agonists and inhibitors were applied to characterize the involvement of various enzymes and Ca(2+)-channels in Ca(2+) signalling pathways.

RESULTS

ACh activating M3-muscarinic receptors and Gβγ protein dependent phosphatidylinositol 3 kinase induces Ca(2+) and NO oscillations in adipocytes. At low concentrations of ACh which are insufficient to induce oscillations, NE or α1, α2-adrenergic agonists act by amplifying the effect of ACh to promote Ca(2+) oscillations or switching phenomena. SNAP, 8-Br-cAMP, NAD and ANP may also produce similar set of dynamic regimes. These regimes arise from activation of the ryanodine receptor (RyR) with the implication of a long positive feedback loop (PFL): Ca(2+)→NO→cGMP→cADPR→Ca(2+), which determines periodic or steady operation of a short PFL based on Ca(2+)-induced Ca(2+) release via RyR by generating cADPR, a coagonist of Ca(2+) at the RyR. Interplay between these two loops may be responsible for the observed effects. Several other PFLs, based on activation of endothelial nitric oxide synthase or of protein kinase B by Ca(2+)-dependent kinases, may reinforce functioning of main PFL and enhance reliability. All observed regimes are independent of operation of the phospholipase C/Ca(2+)-signalling axis, which may be switched off due to negative feedback arising from phosphorylation of the inositol-3-phosphate receptor by protein kinase G.

CONCLUSIONS

This study presents a kinetic model of Ca(2+)-signalling system operating in adipocytes and integrating signals from various agonists, which describes it as multivariable multi feedback network with a family of nested positive feedback.

Using chemical kinetics to model biochemical pathways

Chemical kinetics is the study of the rate of reactions transforming some chemical entities into other chemical entities. Over the twentieth century it has become one of the cornerstones of biochemistry. When in the second half of the century basic knowledge of cellular processes became sufficient to understand quantitatively metabolic networks, chemical kinetics associated with systems theory led to the development of what would become an important branch of systems biology. In this chapter we introduce basic concepts of chemical and enzyme kinetics, and show how the temporal evolution of a reaction system can be described by ordinary differential equations. Finally we present a method to apply this type of approach to model any regulatory network.

P2X7 receptor-mediated calcium dynamics in HEK293 cells: experimental characterization and modelling approach

The P2X7 receptor (P2X7R) induces ionotropic Ca²⁺ signalling in different cell types. It plays an important role in the immune response and in the nervous system. Here, the mechanisms underlying intracellular Ca²⁺ variations evoked by 3'-O-(4-benzoyl)benzoyl-ATP (BzATP), a potent agonist of the P2X7R, in transfected HEK293 cells, are investigated both experimentally and theoretically. We propose a minimal model of P2X7R that is capable of reproducing, qualitatively and quantitatively, the experimental data. This approach was also adopted for the P2X7R variant, which lacks the entire C-terminus tail (trP2X7R). Then we introduce a biophysical model describing the Ca²⁺ dynamics in HEK293. Our model gives an account of the ionotropic Ca²⁺ influx evoked by BzATP on the basis of the kinetics model of P2X7R. To explain the complex Ca²⁺ responses evoked by BzATP, the model predicted that an impairment in Ca²⁺ extrusion flux through the plasma membrane is a key factor for Ca²⁺ homeostasis in HEK293 cells.

Intercellular synchronization of diffusively coupled Ca(2+) oscillators

We examine the synchrony in the dynamics of localized [Ca(2 + )](i) oscillations among a group of cells exhibiting such complex Ca(2 + ) oscillations, connected in the form of long chain, via diffusing coupling where cytosolic Ca(2 + ) and inositol 1,4,5-triphosphate are coupling molecules. Based on our numerical results, we could able to identify three regimes, namely desynchronized, transition and synchronized regimes in the (T - k(e)) (time period-coupling constant) and (A - k(e)) (amplitude-coupling constant) spaces which are supported by phase plots (Δϕ verses time) and recurrence plots, respectively. We further show the increase of synchronization among the cells as the number of coupling molecules increases in the (T - k(e)) and (A - k(e)) spaces.

A graphical method for reducing and relating models in systems biology

Motivation: In Systems Biology, an increasing collection of models of various biological processes is currently developed and made available in publicly accessible repositories, such as biomodels.net for instance, through common exchange formats such as SBML. To date, however, there is no general method to relate different models to each other by abstraction or reduction relationships, and this task is left to the modeler for re-using and coupling models. In mathematical biology, model reduction techniques have been studied for a long time, mainly in the case where a model exhibits different time scales, or different spatial phases, which can be analyzed separately. These techniques are however far too restrictive to be applied on a large scale in systems biology, and do not take into account abstractions other than time or phase decompositions. Our purpose here is to propose a general computational method for relating models together, by considering primarily the structure of the interactions and abstracting from their dynamics in a first step.

Results: We present a graph-theoretic formalism with node merge and delete operations, in which model reductions can be studied as graph matching problems. From this setting, we derive an algorithm for deciding whether there exists a reduction from one model to another, and evaluate it on the computation of the reduction relations between all SBML models of the biomodels.net repository. In particular, in the case of the numerous models of MAPK signalling, and of the circadian clock, biologically meaningful mappings between models of each class are automatically inferred from the structure of the interactions. We conclude on the generality of our graphical method, on its limits with respect to the representation of the structure of the interactions in SBML, and on some perspectives for dealing with the dynamics.

Availability: The algorithms described in this article are implemented in the open-source software modeling platform BIOCHAM available at http://contraintes.inria.fr/biocham The models used in the experiments are available from http://www.biomodels.net/

Contact:francois.fages@inria.fr

A scalable and integrative system for pathway bioinformatics and systems biology

MOTIVATION

Progress in systems biology depends on developing scalable informatics tools to predictively model, visualize, and flexibly store information about complex biological systems. Scalability of these tools, as well as their ability to integrate within larger frameworks of evolving tools, is critical to address the multi-scale and size complexity of biological systems.

RESULTS

Using current software technology, such as self-generation of database and object code from UML schemas, facilitates rapid updating of a scalable expert assistance system for modeling biological pathways. Distribution of key components along with connectivity to external data sources and analysis tools is achieved via a web service interface.

AVAILABILITY

All sigmoid modeling software components and supplementary information are available through: http://www.igb.uci.edu/servers/sb.html.

Cross-talk between signaling pathways can generate robust oscillations in calcium and cAMP

Background

To control and manipulate cellular signaling, we need to understand cellular strategies for information transfer, integration, and decision-making. A key feature of signal transduction is the generation of only a few intracellular messengers by many extracellular stimuli.

Methodology/Principal Findings

Here we model molecular cross-talk between two classic second messengers, cyclic AMP (cAMP) and calcium, and show that the dynamical complexity of the response of both messengers increases substantially through their interaction. In our model of a non-excitable cell, both cAMP and calcium concentrations can oscillate. If mutually inhibitory, cross-talk between the two second messengers can increase the range of agonist concentrations for which oscillations occur. If mutually activating, cross-talk decreases the oscillation range, but can generate ‘bursting’ oscillations of calcium and may enable better filtering of noise.

Conclusion

We postulate that this increased dynamical complexity allows the cell to encode more information, particularly if both second messengers encode signals. In their native environments, it is unlikely that cells are exposed to one stimulus at a time, and cross-talk may help generate sufficiently complex responses to allow the cell to discriminate between different combinations and concentrations of extracellular agonists.

Nonlinear time series analysis of nodulation factor induced calcium oscillations: evidence for deterministic chaos?

Legume plants form beneficial symbiotic interactions with nitrogen fixing bacteria (called rhizobia), with the rhizobia being accommodated in unique structures on the roots of the host plant. The legume/rhizobial symbiosis is responsible for a significant proportion of the global biologically available nitrogen. The initiation of this symbiosis is governed by a characteristic calcium oscillation within the plant root hair cells and this signal is activated by the rhizobia. Recent analyses on calcium time series data have suggested that stochastic effects have a large role to play in defining the nature of the oscillations. The use of multiple nonlinear time series techniques, however, suggests an alternative interpretation, namely deterministic chaos. We provide an extensive, nonlinear time series analysis on the nature of this calcium oscillation response. We build up evidence through a series of techniques that test for determinism, quantify linear and nonlinear components, and measure the local divergence of the system. Chaos is common in nature and it seems plausible that properties of chaotic dynamics might be exploited by biological systems to control processes within the cell. Systems possessing chaotic control mechanisms are more robust in the sense that the enhanced flexibility allows more rapid response to environmental changes with less energetic costs. The desired behaviour could be most efficiently targeted in this manner, supporting some intriguing speculations about nonlinear mechanisms in biological signaling.

Using chemical organization theory for model checking

Motivation: The increasing number and complexity of biomodels makes automatic procedures for checking the models' properties and quality necessary. Approaches like elementary mode analysis, flux balance analysis, deficiency analysis and chemical organization theory (OT) require only the stoichiometric structure of the reaction network for derivation of valuable information. In formalisms like Systems Biology Markup Language (SBML), however, information about the stoichiometric coefficients required for an analysis of chemical organizations can be hidden in kinetic laws.

Results: First, we introduce an algorithm that uncovers stoichiometric information that might be hidden in the kinetic laws of a reaction network. This allows us to apply OT to SBML models using modifiers. Second, using the new algorithm, we performed a large-scale analysis of the 185 models contained in the manually curated BioModels Database. We found that for 41 models (22%) the set of organizations changes when modifiers are considered correctly. We discuss one of these models in detail (BIOMD149, a combined model of the ERK- and Wnt-signaling pathways), whose set of organizations drastically changes when modifiers are considered. Third, we found inconsistencies in 5 models (3%) and identified their characteristics. Compared with flux-based methods, OT is able to identify those species and reactions more accurately [in 26 cases (14%)] that can be present in a long-term simulation of the model. We conclude that our approach is a valuable tool that helps to improve the consistency of biomodels and their repositories.

Availability: All data and a JAVA applet to check SBML-models is available from http://www.minet.uni-jena.de/csb/prj/ot/tools

Contact:dittrich@minet.uni-jena.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Computational modeling of paroxysmal depolarization shifts in neurons induced by the glutamate release from astrocytes

Recent experimental studies have shown that astrocytes respond to external stimuli with a transient increase of the intracellular calcium concentration or can exhibit self-sustained spontaneous activity. Both evoked and spontaneous astrocytic calcium oscillations are accompanied by exocytosis of glutamate caged in astrocytes leading to paroxysmal depolarization shifts (PDS) in neighboring neurons. Here, we present a simple mathematical model of the interaction between astrocytes and neurons that is able to numerically reproduce the experimental results concerning the initiation of the PDS. The timing of glutamate release from the astrocyte is studied by means of a combined modeling of a vesicle cycle and the dynamics of SNARE-proteins. The neuronal slow inward currents (SICs), induced by the astrocytic glutamate and leading to PDS, are modeled via the activation of presynaptic glutamate receptors. The dependence of the bidirectional communication between neurons and astrocytes on the concentration of glutamate transporters is analyzed, as well. Our numerical results are in line with experimental findings showing that astrocyte can induce synchronous PDSs in neighboring neurons, resulting in a transient synchronous spiking activity.

A mesoscopic stochastic mechanism of cytosolic calcium oscillations

Based on a model of intracellular calcium (Ca(2+)) oscillation with self-modulation of inositol 1,4,5-trisphosphate signal, the mesoscopic stochastic differential equations for the intracellular Ca(2+) oscillations are theoretically derived by using the chemical Langevin equation method. The effects of the finite biochemical reaction molecule number on both simple and complex cytosolic Ca(2+) oscillations are numerically studied. In the case of simple intracellular Ca(2+) oscillation, it is found that, with the increase of molecule number, the coherence resonance or autonomous resonance phenomena can occur for some external stimulation parameter values. In the cases of complex cytosolic Ca(2+) oscillations, each extremum of concentration of cytosolic Ca(2+) oscillations corresponds to a peak in the histogram of Ca(2+) concentration, and the most probability appeared during the bursting plateau level for bursting, but at the largest minimum of Ca(2+) concentration for chaos. For quasi-periodicity, however, there are only two peaks in the histogram of Ca(2+) concentration, and the most probability is located at low concentration state.

A quantitative kinetic model for ATP-induced intracellular Ca2+ oscillations

A quantitative kinetic model is proposed to simulate the ATP-induced intracellular Ca(2+) oscillations. The quantitative effect of ATP concentration upon the oscillations was successfully simulated. Our simulation results support previous experimental explanations that the Ca(2+) oscillations are mainly due to interaction of Ca(2+) release from the endoplasmic reticulum (ER) and the ATP-dependent Ca(2+) pump back into the ER, and the oscillations are prolonged by extracellular Ca(2+) entry that maintains the constant Ca(2+) supplies to its intracellular stores. The model is also able to simulate the sudden disappearance phenomenon of the Ca(2+) oscillations observed in some cell types by taking into account of the biphasic characteristic of the Ca(2+) release from the endoplasmic reticulum (ER). Moreover, the model simulation results for the Ca(2+) oscillations characteristics such as duration, peak [Ca(2+)](cyt), and average interval, etc., lead to prediction of some possible factors responsible for the variations of Ca(2+) oscillations in different types of cells.

Differential regulation of proteins by bursting calcium oscillations--a theoretical study

Calcium in ionic form is a second messenger connecting several input signals to several target processes in the cell. The question arises how one second messenger can transmit more than one signal simultaneously (bow-tie structure of signalling). Experimental data on calcium dynamics often show patterns of successive low-peak and high-peak oscillatory phases, known as bursting. Here, we propose that bursting calcium oscillations can perform the function of simultaneous transmission of two signals at physiological calcium concentrations, for example, by selective activation of two calcium-binding proteins. This differential regulation by periodic bursting is investigated in a theoretical model. The two proteins are assumed to be activated by calcium, and one of them is assumed to be subject to biphasic regulation due to additional inhibitory binding sites. To explore which characteristics of the complex signal could be responsible for independent regulation of low-peak activated and spike activated targets, different bursting patterns of simplified square pulses are applied. Depending on the change in the bursting pattern, one protein can be gradually activated at a constant level of the other protein's activity, or the two proteins can be activated simultaneously, or one protein can be activated while the other one is deactivated simultaneously. Thus, the two proteins can be regulated virtually independently.

A model of intracellular Ca2+ oscillations based on the activity of the intermediate-conductance Ca2+-activated K+ channels

Intracellular Ca2+ oscillations are observed in a large number of non-excitable cells. While most appear to reflect an intermittent Ca2+ release from intracellular stores, in some instances intracellular Ca2+ oscillations strongly depend on Ca2+ influx, and are coupled to oscillations of the membrane potential, suggesting that a plasma membrane-based mechanism may be involved. We have developed a theoretical model for the latter type of intracellular Ca2+ oscillations based on the Ca2+-dependent modulation of the intermediate-conductance, Ca2+-activated K+ (IKCa) channel. The functioning of this model relies on the Ca2+-dependent activation, and the much slower Ca2+-dependent rundown of this channel. We have shown that Ca2+-dependent activation of the IKCa channels, the consequent membrane hyperpolarization and the resulting increase in Ca2+ influx may confer the positive feedback mechanism required for the ascending phase of the oscillation. The much slower Ca2+-dependent rundown process will conversely halt this positive loop, and establish the descending phase of the intracellular Ca2+ oscillation. We found that this simple model gives rise to intracellular Ca2+ oscillations when using physiologically reasonable parameters, suggesting that IKCa channels could participate in the generation of intracellular Ca2+ oscillations.

Release currents of IP(3) receptor channel clusters and concentration profiles

We simulate currents and concentration profiles generated by Ca(2+) release from the endoplasmic reticulum (ER) to the cytosol through IP(3) receptor channel clusters. Clusters are described as conducting pores in the lumenal membrane with a diameter from 6 nm to 36 nm. The endoplasmic reticulum is modeled as a disc with a radius of 1-12 microm and an inner height of 28 nm. We adapt the dependence of the currents on the trans Ca(2+) concentration (intralumenal) measured in lipid bilayer experiments to the cellular geometry. Simulated currents are compared with signal mass measurements in Xenopus oocytes. We find that release currents depend linearly on the concentration of free Ca(2+) in the lumen. The release current is approximately proportional to the square root of the number of open channels in a cluster. Cytosolic concentrations at the location of the cluster range from 25 microM to 170 microM. Concentration increase due to puffs in a distance of a few micrometers from the puff site is found to be in the nanomolar range. Release currents decay biexponentially with timescales of <1 s and a few seconds. Concentration profiles decay with timescales of 0.125-0.250 s upon termination of release.

Local dissipation and coupling properties of cellular oscillators: a case study on calcium oscillations

Synchronised signal transduction between cells is crucial, since it assures fast and immutable information processing, which is vital for flawless functioning of living organisms. The question arises how to recognise the ability of a cell to be easily coupled with other cells. In the present paper, we investigate the system properties that determine best coupling abilities and assure the most efficient signal transduction between cells. A case study is done for intercellular calcium oscillations. For a particular diffusion-like coupled system of cellular oscillators, we determined the minimal gap-junctional permeability that is necessary for synchronisation of initially asynchronous oscillators. Our results show that dissipation is a crucial system property that determines the coupling ability of cellular oscillators. We found that low dissipation assures synchronisation of coupled cells already at very low gap-junctional permeability, whereas highly dissipative oscillators require much higher gap-junctional permeability in order to synchronise. The results are discussed in the sense of their biological importance for systems where the synchronous responses of cells were recognised to be indispensable for appropriate physiological functioning of the tissue.

On the encoding and decoding of calcium signals in hepatocytes

Many different agonists use calcium as a second messenger. Despite intensive research in intracellular calcium signalling it is an unsolved riddle how the different types of information represented by the different agonists, is encoded using the universal carrier calcium. It is also still not clear how the information encoded is decoded again into the intracellular specific information at the site of enzymes and genes. After the discovery of calcium oscillations, one likely mechanism is that information is encoded in the frequency, amplitude and waveform of the oscillations. This hypothesis has received some experimental support. However, the mechanism of decoding of oscillatory signals is still not known. Here, we study a mechanistic model of calcium oscillations, which is able to reproduce both spiking and bursting calcium oscillations. We use the model to study the decoding of calcium signals on the basis of co-operativity of calcium binding to various proteins. We show that this co-operativity offers a simple way to decode different calcium dynamics into different enzyme activities.

A theoretical study of effects of cytosolic Ca2+ oscillations on activation of glycogen phosphorylase

Taking into account the Ca(2+)-stimulated degradation of inositol 1,4,5-trisphosphate (IP(3)) by a 3-kinase, we have theoretically explored the effects of both simple and complex Ca(2+) oscillations on the regulation of a phosphorylation-dephosphorylation cycle process involved in glycogen degradation by glycogen phosphorylase a-form, respectively. For the case of simple Ca(2+) oscillations, the roles of cytosolic Ca(2+) oscillations in the regulation of active phosphorylase depend upon the maximum rate of IP(3) degradation by the 3-kinase, V(M5). In particular, the smaller the values of V(M5) are, the lower the effective Ca(2+) threshold for the activation of glycogen phosphorylase will be. For the case of complex Ca(2+) oscillations, the average level of fraction of active phosphorylase is nearly independent from the level of stimulation increasing in the bursting oscillatory domain. Both simple and complex Ca(2+) oscillations can contribute to increase the efficiency and specificity of cellular signalling, and some theoretical results of activation of glycogen phosphorylase regulated by Ca(2+) oscillations are close to the experimental results for gene expression in lymphocytes.

Internal stochastic resonance under two-parameter modulation in intercellular calcium ion oscillations

Internal stochastic resonance (ISR) in a model of intercellular calcium ion oscillations is investigated under the modulation of two parameters, viz., degree of extracellular stimulation (beta) and leak rate (k(f)). ISR can occur when either beta or k(f) is subjected to a noise. Internal stochastic biresonance (ISBR) can occur when noise is added to the two parameters simultaneously. The distance to the bifurcation point is found to be able to enhance or suppress the ISBR, and to affect the number of peaks of ISR.

Under what conditions signal transduction pathways are highly flexible in response to external forcing? A case study on calcium oscillations

Sensitivity and flexibility are typical properties of biological systems. These properties are here investigated in a model for simple and complex intracellular calcium oscillations. In particular, the influence of external periodic forcing is studied. The main point of the study is to compare responses of the system in a chaotic regime with those obtained in a regular periodic regime. We show that the response to external signals in terms of the range of synchronization is not significantly different in regular and chaotic Ca2+ oscillations. However, both types of oscillation are highly flexible in regimes with weak dissipation. Therefore, we conclude that dissipation of free energy is a suitable index characterizing flexibility. For biological systems this appears to be of special importance since for thermodynamic reasons, notably in view of low free energy consumption, dissipation should be minimized.

Sensitivity and flexibility of regular and chaotic calcium oscillations

Sensitivity and flexibility are important properties of biological systems. These properties are here investigated for intracellular calcium oscillations. For a particular model, we comparatively investigate sensitivity and flexibility of regular and chaotic Ca(2+) oscillations. For this model, we obtain two main results. First, sensitivity of the model system to parameter shifting does not depend on the complexity of Ca(2+) oscillations. We observe, however, that both regular and chaotic Ca(2+) oscillations are highly sensitive in regions close to bifurcation points. Second, also flexibility of Ca(2+) oscillations does not significantly depend on the type of Ca(2+) oscillations. Our results show that regular as well as chaotic Ca(2+) oscillations in the studied model are highly flexible in regimes with weak dissipation. Both results are discussed in the sense of possible biological importance.

Modelling of simple and complex calcium oscillations. From single-cell responses to intercellular signalling

This review provides a comparative overview of recent developments in the modelling of cellular calcium oscillations. A large variety of mathematical models have been developed for this wide-spread phenomenon in intra- and intercellular signalling. From these, a general model is extracted that involves six types of concentration variables: inositol 1,4,5-trisphosphate (IP3), cytoplasmic, endoplasmic reticulum and mitochondrial calcium, the occupied binding sites of calcium buffers, and the fraction of active IP3 receptor calcium release channels. Using this framework, the models of calcium oscillations can be classified into 'minimal' models containing two variables and 'extended' models of three and more variables. Three types of minimal models are identified that are all based on calcium-induced calcium release (CICR), but differ with respect to the mechanisms limiting CICR. Extended models include IP3--calcium cross-coupling, calcium sequestration by mitochondria, the detailed gating kinetics of the IP3 receptor, and the dynamics of G-protein activation. In addition to generating regular oscillations, such models can describe bursting and chaotic calcium dynamics. The earlier hypothesis that information in calcium oscillations is encoded mainly by their frequency is nowadays modified in that some effect is attributed to amplitude encoding or temporal encoding. This point is discussed with reference to the analysis of the local and global bifurcations by which calcium oscillations can arise. Moreover, the question of how calcium binding proteins can sense and transform oscillatory signals is addressed. Recently, potential mechanisms leading to the coordination of oscillations in coupled cells have been investigated by mathematical modelling. For this, the general modelling framework is extended to include cytoplasmic and gap-junctional diffusion of IP3 and calcium, and specific models are compared. Various suggestions concerning the physiological significance of oscillatory behaviour in intra- and intercellular signalling are discussed. The article is concluded with a discussion of obstacles and prospects.

Fold–Hopf Bursting in a Model for Calcium Signal Transduction

We study a recent model for calcium signal transduction. This model displays spiking, bursting and chaotic oscillations in accordance with experimental results. We calculate bifurcation diagrams and study the bursting behaviour in detail. This behaviour is classified according to the dynamics of separated slow and fast subsystems. It is shown to be of the Fold–Hopf type, a type which was previously only described in the context of neuronal systems, but not in the context of signal transduction in the cell.

Effects of noise on the off rate of Ca(2+) binding proteins in a coupled biochemical cell system

We investigated numerically the kinetic properties of calcium binding proteins by using a three-calcium store model and discussed the response of a two-way coupled biochemical cell system, whose subsystems were coupled via diffusion-like cytosolic calcium transfer through gap junctions, to the external stimulation. When we used noise to modulate the off rate of the Ca(2+) binding proteins in simulation, an SR-like phenomenon of synchronous oscillations is observed. In addition, the interaction involving noise, coupling and the dynamics of Ca(2+) binding proteins is discussed.

Birhythmicity, trirhythmicity and chaos in bursting calcium oscillations

We have analyzed various types of complex calcium oscillations. The oscillations are explained with a model based on calcium-induced calcium release (CICR). In addition to the endoplasmic reticulum as the main intracellular Ca2+ store, mitochondrial and cytosolic Ca2+ binding proteins are also taken into account. This model was previously proposed for the study of the physiological role of mitochondria and the cytosolic proteins in gene rating complex Ca2+ oscillations [1]. Here, we investigated the occurrence of different types of Ca2+ oscillations obtained by the model, i.e. simple oscillations, bursting, and chaos. In a bifurcation diagram, we have shown that all these various modes of oscillatory behavior are obtained by a change of only one model parameter, which corresponds to the physiological variability of an agonist. Bursting oscillations were studied in more detail because they express birhythmicity, trirhythmicity and chaotic behavior. Two different routes to chaos are observed in the model: in addition to the usual period doubling cascade, we also show intermittency. For the characterization of the chaotic behavior, we made use of return maps and Lyapunov exponents. The potential biological role of chaos in intracellular signaling is discussed.

From simple to complex oscillatory behavior in metabolic and genetic control networks

We present an overview of mechanisms responsible for simple or complex oscillatory behavior in metabolic and genetic control networks. Besides simple periodic behavior corresponding to the evolution toward a limit cycle we consider complex modes of oscillatory behavior such as complex periodic oscillations of the bursting type and chaos. Multiple attractors are also discussed, e.g., the coexistence between a stable steady state and a stable limit cycle (hard excitation), or the coexistence between two simultaneously stable limit cycles (birhythmicity). We discuss mechanisms responsible for the transition from simple to complex oscillatory behavior by means of a number of models serving as selected examples. The models were originally proposed to account for simple periodic oscillations observed experimentally at the cellular level in a variety of biological systems. In a second stage, these models were modified to allow for complex oscillatory phenomena such as bursting, birhythmicity, or chaos. We consider successively (1) models based on enzyme regulation, proposed for glycolytic oscillations and for the control of successive phases of the cell cycle, respectively; (2) a model for intracellular Ca(2+) oscillations based on transport regulation; (3) a model for oscillations of cyclic AMP based on receptor desensitization in Dictyostelium cells; and (4) a model based on genetic regulation for circadian rhythms in Drosophila. Two main classes of mechanism leading from simple to complex oscillatory behavior are identified, namely (i) the interplay between two endogenous oscillatory mechanisms, which can take multiple forms, overt or more subtle, depending on whether the two oscillators each involve their own regulatory feedback loop or share a common feedback loop while differing by some related process, and (ii) self-modulation of the oscillator through feedback from the system's output on one of the parameters controlling oscillatory behavior. However, the latter mechanism may also be viewed as involving the interplay between two feedback processes, each of which might be capable of producing oscillations. Although our discussion primarily focuses on the case of autonomous oscillatory behavior, we also consider the case of nonautonomous complex oscillations in a model for circadian oscillations subjected to periodic forcing by a light-dark cycle and show that the occurrence of entrainment versus chaos in these conditions markedly depends on the wave form of periodic forcing. (c) 2001 American Institute of Physics.

Excitable calcium wave propagation in the presence of localized stores

We study the propagation of calcium waves in the presence of a discrete distribution of calcium stores. Calcium-induced calcium release coupled to diffusion can be used to produce a criterion for wave propagation across connected clusters of stores. The velocity of the resulting wave and its relationship to the frequency of the excitatory stimulus can then be described using percolation theory. Simulations show a homogeneous and a fractal regime and are in agreement with both experiments and theory.

Switching from simple to complex oscillations in calcium signaling

We present a new model for calcium oscillations based on experiments in hepatocytes. The model considers feedback inhibition on the initial agonist receptor complex by calcium and activated phospholipase C, as well as receptor type-dependent self-enhanced behavior of the activated G(alpha) subunit. It is able to show simple periodic oscillations and periodic bursting, and it is the first model to display chaotic bursting in response to agonist stimulations. Moreover, our model offers a possible explanation for the differences in dynamic behavior observed in response to different agonists in hepatocytes.

Complex calcium oscillations and the role of mitochondria and cytosolic proteins

Intracellular calcium oscillations, which are oscillatory changes of cytosolic calcium concentration in response to agonist stimulation, are experimentally well observed in various living cells. Simple calcium oscillations represent the most common pattern and many mathematical models have been published to describe this type of oscillation. On the other hand, relatively few theoretical studies have been proposed to give an explanation of complex intracellular calcium oscillations, such as bursting and chaos. In this paper, we develop a new possible mechanism for complex calcium oscillations based on the interplay between three calcium stores in the cell: the endoplasmic reticulum (ER), mitochondria and cytosolic proteins. The majority ( approximately 80%) of calcium released from the ER is first very quickly sequestered by mitochondria. Afterwards, a much slower release of calcium from the mitochondria serves as the calcium supply for the intermediate calcium exchanges between the ER and the cytosolic proteins causing bursting calcium oscillations. Depending on the permeability of the ER channels and on the kinetic properties of calcium binding to the cytosolic proteins, different patterns of complex calcium oscillations appear. With our model, we are able to explain simple calcium oscillations, bursting and chaos. Chaos is also observed for calcium oscillations in the bursting mode.