1
|
Friedhoff VN, Lindner B, Falcke M. Modeling IP 3-induced Ca 2+ signaling based on its interspike interval statistics. Biophys J 2023; 122:2818-2831. [PMID: 37312455 PMCID: PMC10398346 DOI: 10.1016/j.bpj.2023.06.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Revised: 04/24/2023] [Accepted: 06/06/2023] [Indexed: 06/15/2023] Open
Abstract
Inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ signaling is a second messenger system used by almost all eukaryotic cells. Recent research demonstrated randomness of Ca2+ signaling on all structural levels. We compile eight general properties of Ca2+ spiking common to all cell types investigated and suggest a theory of Ca2+ spiking starting from the random behavior of IP3 receptor channel clusters mediating the release of Ca2+ from the endoplasmic reticulum capturing all general properties and pathway-specific behavior. Spike generation begins after the absolute refractory period of the previous spike. According to its hierarchical spreading from initiating channel openings to cell level, we describe it as a first passage process from none to all clusters open while the cell recovers from the inhibition which terminated the previous spike. Our theory reproduces the exponential stimulation response relation of the average interspike interval Tav and its robustness properties, random spike timing with a linear moment relation between Tav and the interspike interval SD and its robustness properties, sensitive dependency of Tav on diffusion properties, and nonoscillatory local dynamics. We explain large cell variability of Tav observed in experiments by variability of channel cluster coupling by Ca2+-induced Ca2+ release, the number of clusters, and IP3 pathway component expression levels. We predict the relation between puff probability and agonist concentration and [IP3] and agonist concentration. Differences of spike behavior between cell types and stimulating agonists are explained by the different types of negative feedback terminating spikes. In summary, the hierarchical random character of spike generation explains all of the identified general properties.
Collapse
Affiliation(s)
- Victor Nicolai Friedhoff
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany; Department of Physics, Humboldt University, Berlin, Germany
| | - Benjamin Lindner
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany; Department of Physics, Humboldt University, Berlin, Germany
| | - Martin Falcke
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany; Department of Physics, Humboldt University, Berlin, Germany.
| |
Collapse
|
2
|
Ramlow L, Falcke M, Lindner B. An integrate-and-fire approach to Ca 2+ signaling. Part I: Renewal model. Biophys J 2023; 122:713-736. [PMID: 36635961 PMCID: PMC9989887 DOI: 10.1016/j.bpj.2023.01.007] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Revised: 12/13/2022] [Accepted: 01/06/2023] [Indexed: 01/13/2023] Open
Abstract
In computational neuroscience integrate-and-fire models capture the spike generation by a subthreshold dynamics supplemented by a simple fire-and-reset rule; they allow for a numerically efficient and analytically tractable description of stochastic single cell as well as network dynamics. Stochastic spiking is also a prominent feature of Ca2+ signaling which suggests to adopt the integrate-and-fire approach for this fundamental biophysical process. The model introduced here consists of two components describing 1) activity of clusters of inositol-trisphosphate receptor channels and 2) dynamics of the global Ca2+ concentrations in the cytosol. The cluster dynamics is given in terms of a cyclic Markov chain, capturing the puff, i.e., the punctuated release of Ca2+ from intracellular stores. The cytosolic Ca2+ concentration is described by an integrate-and-fire dynamics driven by the puff current. For the cyclic Markov chain we derive expressions for the statistics of the interpuff interval, the single-puff strength and the puff current assuming constant cytosolic Ca2+. The latter condition is often well approximated because cytosolic Ca2+ varies much slower than the cluster activity does. Furthermore, because the detailed two-component model is numerically expensive to simulate and difficult to treat analytically, we develop an analytical framework to approximate the driving puff current of the stochastic cytosolic Ca2+ dynamics by a temporally uncorrelated Gaussian noise. This approximation reduces our two-component system to an integrate-and-fire model with a nonlinear drift function and a multiplicative Gaussian white noise, a model that is known to generate a renewal spike train, i.e., a point process with statistically independent interspike intervals. The model allows for fast numerical simulations, permits to derive analytical expressions for the rate of Ca2+ spiking and the coefficient of variation of the interspike interval, and to approximate the interspike interval density and the spike train power spectrum. Comparison of these statistics to experimental data is discussed.
Collapse
Affiliation(s)
- Lukas Ramlow
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany; Physics Department of Humboldt University Berlin, Berlin, Germany.
| | - Martin Falcke
- Physics Department of Humboldt University Berlin, Berlin, Germany; Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Benjamin Lindner
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany; Physics Department of Humboldt University Berlin, Berlin, Germany
| |
Collapse
|
3
|
Stochastic reaction-diffusion modeling of calcium dynamics in 3D dendritic spines of Purkinje cells. Biophys J 2021; 120:2112-2123. [PMID: 33887224 PMCID: PMC8390834 DOI: 10.1016/j.bpj.2021.03.027] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Revised: 02/22/2021] [Accepted: 03/18/2021] [Indexed: 02/07/2023] Open
Abstract
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.
Collapse
|
4
|
Zhou A, Liu X, Zhang S, Huo B. Effects of store-operated and receptor-operated calcium channels on synchronization of calcium oscillations in astrocytes. Biosystems 2020; 198:104233. [PMID: 32858094 DOI: 10.1016/j.biosystems.2020.104233] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2019] [Revised: 08/04/2020] [Accepted: 08/20/2020] [Indexed: 10/23/2022]
Abstract
Intercellular calcium signaling allows cells to communicate with each other and to interact with adjacent cells. Gap junction is the most common and important way for cellular communication. Recently, mathematical models have been widely used to gain a precise and quantitative understanding of the dynamics of intracellular calcium ions (Ca2+). In this paper, we establish a mathematical model considering the gap junction permeable to Ca2+ and to IP3 for describing the calcium oscillations in coupled astrocytes. Store-operated calcium entry (SOCE) is viewed as the main process which controls the non-excitable cells, hence, we focus on the effect of store-operated calcium channel (SOCC) and receptor-operated calcium channel (ROCC) on the intercellular synchronization, respectively. By employing bifurcation analysis on this model, the dynamic behaviors of the coupled system with different physiological state cells is obtained with changes in the maximum capacity of the SOCC and the ROCC. The synchronization boundaries for different conditions are gained in the two parameters space of the channel parameters and the coupling strength. The results suggest that the variation of the maximum flow for different calcium channels determines the stable oscillations of the coupled system, as well as for the frequency and amplitude of oscillations. The SOCC has an expected effect on the change of the oscillatory interval while the ROCC demonstrated the influence on the amplitude modulation. Furthermore, the coupling strength and channel parameters could induce 1:1 locking of intercellular Ca2+ oscillations and the synchronization region like Arnol'd tongue is found.
Collapse
Affiliation(s)
- Anqi Zhou
- Department of Mechanics, Tianjin University, Tianjin, 300354, PR China
| | - Xijun Liu
- Department of Mechanics, Tianjin University, Tianjin, 300354, PR China
| | - Suxia Zhang
- Department of Mechanics, Tianjin University, Tianjin, 300354, PR China.
| | - Bing Huo
- College of Mechanical Engineering, Tianjin University of Science & Technology, Tianjin, 300222, PR China
| |
Collapse
|
5
|
Falcke M, Friedhoff VN. The stretch to stray on time: Resonant length of random walks in a transient. CHAOS (WOODBURY, N.Y.) 2018; 28:053117. [PMID: 29857685 DOI: 10.1063/1.5023164] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
First-passage times in random walks have a vast number of diverse applications in physics, chemistry, biology, and finance. In general, environmental conditions for a stochastic process are not constant on the time scale of the average first-passage time or control might be applied to reduce noise. We investigate moments of the first-passage time distribution under an exponential transient describing relaxation of environmental conditions. We solve the Laplace-transformed (generalized) master equation analytically using a novel method that is applicable to general state schemes. The first-passage time from one end to the other of a linear chain of states is our application for the solutions. The dependence of its average on the relaxation rate obeys a power law for slow transients. The exponent ν depends on the chain length N like ν=-N/(N+1) to leading order. Slow transients substantially reduce the noise of first-passage times expressed as the coefficient of variation (CV), even if the average first-passage time is much longer than the transient. The CV has a pronounced minimum for some lengths, which we call resonant lengths. These results also suggest a simple and efficient noise control strategy and are closely related to the timing of repetitive excitations, coherence resonance, and information transmission by noisy excitable systems. A resonant number of steps from the inhibited state to the excitation threshold and slow recovery from negative feedback provide optimal timing noise reduction and information transmission.
Collapse
Affiliation(s)
- Martin Falcke
- Max Delbrück Center for Molecular Medicine, Robert Rössle Str. 10, 13125 Berlin, Germany
| | | |
Collapse
|
6
|
Modelling mechanism of calcium oscillations in pancreatic acinar cells. J Bioenerg Biomembr 2014; 46:403-20. [DOI: 10.1007/s10863-014-9561-0] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2014] [Accepted: 06/24/2014] [Indexed: 01/18/2023]
|
7
|
Cao P, Donovan G, Falcke M, Sneyd J. A stochastic model of calcium puffs based on single-channel data. Biophys J 2013; 105:1133-42. [PMID: 24010656 PMCID: PMC3852038 DOI: 10.1016/j.bpj.2013.07.034] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2013] [Revised: 07/03/2013] [Accepted: 07/24/2013] [Indexed: 11/17/2022] Open
Abstract
Calcium puffs are local transient Ca(2+) releases from internal Ca(2+) stores such as the endoplasmic reticulum or the sarcoplasmic reticulum. Such release occurs through a cluster of inositol 1,4,5-trisphosphate receptors (IP3Rs). Based on the IP3R model (which is determined by fitting to stationary single-channel data) and nonstationary single-channel data, we construct a new IP3R model that includes time-dependent rates of mode switches. A point-source model of Ca(2+) puffs is then constructed based on the new IP3R model and is solved by a hybrid Gillespie method with adaptive timing. Model results show that a relatively slow recovery of an IP3R from Ca(2+) inhibition is necessary to reproduce most of the experimental outcomes, especially the nonexponential interpuff interval distributions. The number of receptors in a cluster could be severely underestimated when the recovery is sufficiently slow. Furthermore, we find that, as the number of IP3Rs increases, the average duration of puffs initially increases but then becomes saturated, whereas the average decay time keeps increasing linearly. This gives rise to the observed asymmetric puff shape.
Collapse
Affiliation(s)
- Pengxing Cao
- Department of Mathematics, The University of Auckland, Auckland, New Zealand
| | - Graham Donovan
- Department of Mathematics, The University of Auckland, Auckland, New Zealand
| | - Martin Falcke
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - James Sneyd
- Department of Mathematics, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
8
|
Comparison of models for IP3 receptor kinetics using stochastic simulations. PLoS One 2013; 8:e59618. [PMID: 23630568 PMCID: PMC3629942 DOI: 10.1371/journal.pone.0059618] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2012] [Accepted: 02/15/2013] [Indexed: 12/07/2022] Open
Abstract
Inositol 1,4,5-trisphosphate receptor (IP3R) is a ubiquitous intracellular calcium (Ca2+) channel which has a major role in controlling Ca2+ levels in neurons. A variety of computational models have been developed to describe the kinetic function of IP3R under different conditions. In the field of computational neuroscience, it is of great interest to apply the existing models of IP3R when modeling local Ca2+ transients in dendrites or overall Ca2+ dynamics in large neuronal models. The goal of this study was to evaluate existing IP3R models, based on electrophysiological data. This was done in order to be able to suggest suitable models for neuronal modeling. Altogether four models (Othmer and Tang, 1993; Dawson etal., 2003; Fraiman and Dawson, 2004; Doi etal., 2005) were selected for a more detailed comparison. The selection was based on the computational efficiency of the models and the type of experimental data that was used in developing the model. The kinetics of all four models were simulated by stochastic means, using the simulation software STEPS, which implements the Gillespie stochastic simulation algorithm. The results show major differences in the statistical properties of model functionality. Of the four compared models, the one by Fraiman and Dawson (2004) proved most satisfactory in producing the specific features of experimental findings reported in literature. To our knowledge, the present study is the first detailed evaluation of IP3R models using stochastic simulation methods, thus providing an important setting for constructing a new, realistic model of IP3R channel kinetics for compartmental modeling of neuronal functions. We conclude that the kinetics of IP3R with different concentrations of Ca2+ and IP3 should be more carefully addressed when new models for IP3R are developed.
Collapse
|
9
|
Actin-based confinement of calcium responses during Shigella invasion. Nat Commun 2013; 4:1567. [DOI: 10.1038/ncomms2561] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2012] [Accepted: 01/28/2013] [Indexed: 12/18/2022] Open
|
10
|
Moenke G, Falcke M, Thurley K. Hierarchic stochastic modelling applied to intracellular Ca(2+) signals. PLoS One 2012; 7:e51178. [PMID: 23300536 PMCID: PMC3531454 DOI: 10.1371/journal.pone.0051178] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2012] [Accepted: 10/30/2012] [Indexed: 11/19/2022] Open
Abstract
Important biological processes like cell signalling and gene expression have noisy components and are very complex at the same time. Mathematical analysis of such systems has often been limited to the study of isolated subsystems, or approximations are used that are difficult to justify. Here we extend a recently published method (Thurley and Falcke, PNAS 2011) which is formulated in observable system configurations instead of molecular transitions. This reduces the number of system states by several orders of magnitude and avoids fitting of kinetic parameters. The method is applied to Ca(2+) signalling. Ca(2+) is a ubiquitous second messenger transmitting information by stochastic sequences of concentration spikes, which arise by coupling of subcellular Ca(2+) release events (puffs). We derive analytical expressions for a mechanistic Ca(2+) model, based on recent data from live cell imaging, and calculate Ca(2+) spike statistics in dependence on cellular parameters like stimulus strength or number of Ca(2+) channels. The new approach substantiates a generic Ca(2+) model, which is a very convenient way to simulate Ca(2+) spike sequences with correct spiking statistics.
Collapse
Affiliation(s)
- Gregor Moenke
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Martin Falcke
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Keven Thurley
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
- Institute for Theoretical Biology, Charité Universitätsmedizin, Berlin, Germany
| |
Collapse
|
11
|
Siekmann I, Wagner LE, Yule D, Crampin EJ, Sneyd J. A kinetic model for type I and II IP3R accounting for mode changes. Biophys J 2012; 103:658-68. [PMID: 22947927 PMCID: PMC3443778 DOI: 10.1016/j.bpj.2012.07.016] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2012] [Revised: 07/11/2012] [Accepted: 07/13/2012] [Indexed: 10/28/2022] Open
Abstract
Based upon an extensive single-channel data set, a Markov model for types I and II inositol trisphosphate receptors (IP(3)R) is developed. The model aims to represent accurately the kinetics of both receptor types of IP(3)R depending on the concentrations of inositol trisphosphate (IP(3)), adenosine trisphosphate (ATP), and intracellular calcium (Ca(2+)). In particular, the model takes into account that for some combinations of ligands the IP(3)R switches between extended periods of inactivity alternating with intervals of bursting activity (mode changes). In a first step, the inactive and active modes are modeled separately. It is found that, within modes, both receptor types are ligand-independent. In a second step, the submodels are connected by transition rates. Ligand-dependent regulation of the channel activity is achieved by modulating these transitions between active and inactive modes. As a result, a compact representation of the IP(3)R is obtained that accurately captures stochastic single-channel dynamics including mode changes in a model with six states and 10 rate constants, only two of which are ligand-dependent.
Collapse
Affiliation(s)
- Ivo Siekmann
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand.
| | | | | | | | | |
Collapse
|
12
|
Thurley K, Smith IF, Tovey SC, Taylor CW, Parker I, Falcke M. Timescales of IP(3)-evoked Ca(2+) spikes emerge from Ca(2+) puffs only at the cellular level. Biophys J 2012; 101:2638-44. [PMID: 22261051 DOI: 10.1016/j.bpj.2011.10.030] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2011] [Revised: 09/29/2011] [Accepted: 10/03/2011] [Indexed: 10/14/2022] Open
Abstract
The behavior of biological systems is determined by the properties of their component molecules, but the interactions are usually too complex to understand fully how molecular behavior generates cellular behavior. Ca(2+) signaling by inositol trisphosphate receptors (IP(3)R) offers an opportunity to understand this relationship because the cellular behavior is defined largely by Ca(2+)-mediated interactions between IP(3)R. Ca(2+) released by a cluster of IP(3)R (giving a local Ca(2+) puff) diffuses and ignites the behavior of neighboring clusters (to give repetitive global Ca(2+) spikes). We use total internal reflection fluorescence microscopy of two mammalian cell lines to define the temporal relationships between Ca(2+) puffs (interpuff intervals, IPI) and Ca(2+) spikes (interspike intervals) evoked by flash photolysis of caged IP(3). We find that IPI are much shorter than interspike intervals, that puff activity is stochastic with a recovery time that is much shorter than the refractory period of the cell, and that IPI are not periodic. We conclude that Ca(2+) spikes do not arise from oscillatory dynamics of IP(3)R clusters, but that repetitive Ca(2+) spiking with its longer timescales is an emergent property of the dynamics of the whole cluster array.
Collapse
Affiliation(s)
- Kevin Thurley
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Berlin, Germany.
| | | | | | | | | | | |
Collapse
|
13
|
Fundamental properties of Ca2+ signals. Biochim Biophys Acta Gen Subj 2011; 1820:1185-94. [PMID: 22040723 DOI: 10.1016/j.bbagen.2011.10.007] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2011] [Revised: 10/16/2011] [Accepted: 10/17/2011] [Indexed: 11/21/2022]
Abstract
BACKGROUND Ca2+ is a ubiquitous and versatile second messenger that transmits information through changes of the cytosolic Ca2+ concentration. Recent investigations changed basic ideas on the dynamic character of Ca2+ signals and challenge traditional ideas on information transmission. SCOPE OF REVIEW We present recent findings on key characteristics of the cytosolic Ca2+ dynamics and theoretical concepts that explain the wide range of experimentally observed Ca2+ signals. Further, we relate properties of the dynamical regulation of the cytosolic Ca2+ concentration to ideas about information transmission by stochastic signals. MAJOR CONCLUSIONS We demonstrate the importance of the hierarchal arrangement of Ca2+ release sites on the emergence of cellular Ca2+ spikes. Stochastic Ca2+ signals are functionally robust and adaptive to changing environmental conditions. Fluctuations of interspike intervals (ISIs) and the moment relation derived from ISI distributions contain information on the channel cluster open probability and on pathway properties. GENERAL SIGNIFICANCE Robust and reliable signal transduction pathways that entail Ca2+ dynamics are essential for eukaryotic organisms. Moreover, we expect that the design of a stochastic mechanism which provides robustness and adaptivity will be found also in other biological systems. Ca2+ dynamics demonstrate that the fluctuations of cellular signals contain information on molecular behavior. This article is part of a Special Issue entitled Biochemical, biophysical and genetic approaches to intracellular calcium signaling.
Collapse
|
14
|
LaMar MD, Kemper P, Smith GD. Reduction of calcium release site models via moment fitting of phase-type distributions. Phys Biol 2011; 8:026015. [PMID: 21471635 DOI: 10.1088/1478-3975/8/2/026015] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Models of calcium (Ca(2 +)) release sites derived from continuous-time Markov chain (CTMC) models of intracellular Ca(2 +) channels exhibit collective gating reminiscent of the experimentally observed phenomenon of Ca(2 +) puffs and sparks. In order to overcome the state-space explosion that occurs in compositionally defined Ca(2 +) release site models, we have implemented an automated procedure for model reduction that replaces aggregated states of the full release site model with much simpler CTMCs that have similar within-group phase-type sojourn times and inter-group transitions. Error analysis based on comparison of full and reduced models validates the method when applied to release site models composed of 20 three-state channels that are both activated and inactivated by Ca(2 +). Although inspired by existing techniques for fitting moments of phase-type distributions, the automated reduction method for compositional Ca(2 +) release site models is unique in several respects and novel in this biophysical context.
Collapse
Affiliation(s)
- M Drew LaMar
- Department of Applied Science, The College of William and Mary, Williamsburg, VA 23187, USA.
| | | | | |
Collapse
|
15
|
|
16
|
Skupin A, Kettenmann H, Falcke M. Calcium signals driven by single channel noise. PLoS Comput Biol 2010; 6:e1000870. [PMID: 20700497 PMCID: PMC2917103 DOI: 10.1371/journal.pcbi.1000870] [Citation(s) in RCA: 78] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2009] [Accepted: 06/29/2010] [Indexed: 12/21/2022] Open
Abstract
Usually, the occurrence of random cell behavior is appointed to small copy numbers of molecules involved in the stochastic process. Recently, we demonstrated for a variety of cell types that intracellular Ca2+ oscillations are sequences of random spikes despite the involvement of many molecules in spike generation. This randomness arises from the stochastic state transitions of individual Ca2+ release channels and does not average out due to the existence of steep concentration gradients. The system is hierarchical due to the structural levels channel--channel cluster--cell and a corresponding strength of coupling. Concentration gradients introduce microdomains which couple channels of a cluster strongly. But they couple clusters only weakly; too weak to establish deterministic behavior on cell level. Here, we present a multi-scale modelling concept for stochastic hierarchical systems. It simulates active molecules individually as Markov chains and their coupling by deterministic diffusion. Thus, we are able to follow the consequences of random single molecule state changes up to the signal on cell level. To demonstrate the potential of the method, we simulate a variety of experiments. Comparisons of simulated and experimental data of spontaneous oscillations in astrocytes emphasize the role of spatial concentration gradients in Ca2+ signalling. Analysis of extensive simulations indicates that frequency encoding described by the relation between average and standard deviation of interspike intervals is surprisingly robust. This robustness is a property of the random spiking mechanism and not a result of control.
Collapse
Affiliation(s)
- Alexander Skupin
- Max-Planck-Institute of Molecular Plant Physiology, Potsdam, Germany.
| | | | | |
Collapse
|
17
|
Song SO, Varner J. Modeling and analysis of the molecular basis of pain in sensory neurons. PLoS One 2009; 4:e6758. [PMID: 19750220 PMCID: PMC2735677 DOI: 10.1371/journal.pone.0006758] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2009] [Accepted: 07/23/2009] [Indexed: 12/04/2022] Open
Abstract
Intracellular calcium dynamics are critical to cellular functions like pain transmission. Extracellular ATP plays an important role in modulating intracellular calcium levels by interacting with the P2 family of surface receptors. In this study, we developed a mechanistic mathematical model of ATP-induced P2 mediated calcium signaling in archetype sensory neurons. The model architecture, which described 90 species connected by 162 interactions, was formulated by aggregating disparate molecular modules from literature. Unlike previous models, only mass action kinetics were used to describe the rate of molecular interactions. Thus, the majority of the 252 unknown model parameters were either association, dissociation or catalytic rate constants. Model parameters were estimated from nine independent data sets taken from multiple laboratories. The training data consisted of both dynamic and steady-state measurements. However, because of the complexity of the calcium network, we were unable to estimate unique model parameters. Instead, we estimated a family or ensemble of probable parameter sets using a multi-objective thermal ensemble method. Each member of the ensemble met an error criterion and was located along or near the optimal trade-off surface between the individual training data sets. The model quantitatively reproduced experimental measurements from dorsal root ganglion neurons as a function of extracellular ATP forcing. Hypothesized architecture linking phosphoinositide regulation with P2X receptor activity explained the inhibition of P2X-mediated current flow by activated metabotropic P2Y receptors. Sensitivity analysis using individual and the whole system outputs suggested which molecular subsystems were most important following P2 activation. Taken together, modeling and analysis of ATP-induced P2 mediated calcium signaling generated qualitative insight into the critical interactions controlling ATP induced calcium dynamics. Understanding these critical interactions may prove useful for the design of the next generation of molecular pain management strategies.
Collapse
Affiliation(s)
- Sang Ok Song
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York, United States of America
| | - Jeffrey Varner
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York, United States of America
- * E-mail:
| |
Collapse
|
18
|
Gin E, Wagner LE, Yule DI, Sneyd J. Inositol trisphosphate receptor and ion channel models based on single-channel data. CHAOS (WOODBURY, N.Y.) 2009; 19:037104. [PMID: 19792029 PMCID: PMC5848693 DOI: 10.1063/1.3184540] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/11/2009] [Accepted: 07/01/2009] [Indexed: 05/28/2023]
Abstract
The inositol trisphosphate receptor (IPR) plays an important role in controlling the dynamics of intracellular Ca(2+). Single-channel patch-clamp recordings are a typical way to study these receptors as well as other ion channels. Methods for analyzing and using this type of data have been developed to fit Markov models of the receptor. The usual method of parameter fitting is based on maximum-likelihood techniques. However, Bayesian inference and Markov chain Monte Carlo techniques are becoming more popular. We describe the application of the Bayesian methods to real experimental single-channel data in three ion channels: the ryanodine receptor, the K(+) channel, and the IPR. One of the main aims of all three studies was that of model selection with different approaches taken. We also discuss the modeling implications for single-channel data that display different levels of channel activity within one recording.
Collapse
Affiliation(s)
- Elan Gin
- Department of Mathematics, The University of Auckland, Auckland, New Zealand
| | | | | | | |
Collapse
|
19
|
Thul R, Thurley K, Falcke M. Toward a predictive model of Ca2+ puffs. CHAOS (WOODBURY, N.Y.) 2009; 19:037108. [PMID: 19792033 DOI: 10.1063/1.3183809] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
We investigate the key characteristics of Ca(2+) puffs in deterministic and stochastic frameworks that all incorporate the cellular morphology of IP(3) receptor channel clusters. In the first step, we numerically study the Ca(2+) liberation in a three-dimensional representation of a cluster environment with reaction-diffusion dynamics in both the cytosol and the lumen. These simulations reveal that Ca(2+) concentrations at a releasing cluster range from 80 to 170 microM and equilibrate almost instantaneously on the time scale of the release duration. These highly elevated Ca(2+) concentrations eliminate Ca(2+) oscillations in a deterministic model of an IP(3)R channel cluster at physiological parameter values as revealed by a linear stability analysis. The reason lies in the saturation of all feedback processes in the IP(3)R gating dynamics, so that only fluctuations can restore experimentally observed Ca(2+) oscillations. In this spirit, we derive master equations that allow us to analytically quantify the onset of Ca(2+) puffs and hence the stochastic time scale of intracellular Ca(2+) dynamics. Moving up the spatial scale, we suggest to formulate cellular dynamics in terms of waiting time distribution functions. This approach prevents the state space explosion that is typical for the description of cellular dynamics based on channel states and still contains information on molecular fluctuations. We illustrate this method by studying global Ca(2+) oscillations.
Collapse
Affiliation(s)
- R Thul
- School of Mathematical Sciences, University of Nottingham, Nottingham, United Kingdom
| | | | | |
Collapse
|
20
|
Shuai JW, Yang DP, Pearson JE, Rüdiger S. An investigation of models of the IP3R channel in Xenopus oocyte. CHAOS (WOODBURY, N.Y.) 2009; 19:037105. [PMID: 19792030 PMCID: PMC2771705 DOI: 10.1063/1.3156402] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2009] [Accepted: 05/29/2009] [Indexed: 05/28/2023]
Abstract
We consider different models of inositol 1,4,5-trisphosphate (IP(3)) receptor (IP(3)R) channels in order to fit nuclear membrane patch clamp data of the stationary open probability, mean open time, and mean close time of channels in the Xenopus oocyte. Our results indicate that rather than to treat the tetrameric IP(3)R as four independent and identical subunits, one should assume sequential binding-unbinding processes of Ca(2+) ions and IP(3) messengers. Our simulations also favor the assumption that a channel opens through a conformational transition from a close state to an active state.
Collapse
Affiliation(s)
- J W Shuai
- Department of Physics, Xiamen University, Xiamen, China.
| | | | | | | |
Collapse
|
21
|
Warren NJ, Tawhai MH, Crampin EJ. Mathematical modelling of calcium wave propagation in mammalian airway epithelium: evidence for regenerative ATP release. Exp Physiol 2009; 95:232-49. [PMID: 19700517 DOI: 10.1113/expphysiol.2009.049585] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Airway epithelium has been shown to exhibit intracellular calcium waves after mechanical stimulation. Two classes of mechanism have been proposed to explain calcium wave propagation: diffusion through gap junctions of the intracellular messenger inositol 1,4,5-trisphosphate (IP3), and diffusion of paracrine extracellular messengers such as ATP. We have used single cell recordings of airway epithelium to parameterize a model of an airway epithelial cell. This was then incorporated into a spatial model of a cell culture where both mechanisms for calcium wave propagation are possible. It is shown that a decreasing return on the radius of Ca2+ wave propagation is achieved as the amount of ATP released from the stimulated cell increases. It is therefore shown that for a Ca2+ wave to propagate large distances, a significant fraction of the intracellular ATP pool would be required to be released. Further to this, the radial distribution of maximal calcium response from the stimulated cell does not produce the same flat profile of maximal calcium response seen in experiential studies. This suggests that an additional mechanism is important in Ca2+ wave propagation, such as regenerative release of ATP from cells downstream of the stimulated cell.
Collapse
Affiliation(s)
- N J Warren
- Auckland Bioengineering Institute, Level 6, 70 Symonds Street, CBD, Auckland, New Zealand.
| | | | | |
Collapse
|
22
|
Modeling Ca2+ feedback on a single inositol 1,4,5-trisphosphate receptor and its modulation by Ca2+ buffers. Biophys J 2008; 95:3738-52. [PMID: 18641077 DOI: 10.1529/biophysj.108.137182] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The inositol 1,4,5-trisphosphate receptor/channel (IP(3)R) is a major regulator of intracellular Ca(2+) signaling, and liberates Ca(2+) ions from the endoplasmic reticulum in response to binding at cytosolic sites for both IP(3) and Ca(2+). Although the steady-state gating properties of the IP(3)R have been extensively studied and modeled under conditions of fixed [IP(3)] and [Ca(2+)], little is known about how Ca(2+) flux through a channel may modulate the gating of that same channel by feedback onto activating and inhibitory Ca(2+) binding sites. We thus simulated the dynamics of Ca(2+) self-feedback on monomeric and tetrameric IP(3)R models. A major conclusion is that self-activation depends crucially on stationary cytosolic Ca(2+) buffers that slow the collapse of the local [Ca(2+)] microdomain after closure. This promotes burst-like reopenings by the rebinding of Ca(2+) to the activating site; whereas inhibitory actions are substantially independent of stationary buffers but are strongly dependent on the location of the inhibitory Ca(2+) binding site on the IP(3)R in relation to the channel pore.
Collapse
|
23
|
Thul R, Bellamy TC, Roderick HL, Bootman MD, Coombes S. Calcium oscillations. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2008; 641:1-27. [PMID: 18783168 DOI: 10.1007/978-0-387-09794-7_1] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Changes in cellular Ca2+ concentration control a wide range of physiological processes, from the subsecond release of synaptic neurotransmitters, to the regulation of gene expression over months or years. Ca2+ can also trigger cell death through both apoptosis and necrosis, and so the regulation of cellular Ca2+ concentration must be tightly controlled through the concerted action of pumps, channels and buffers that transport Ca2+ into and out of the cell cytoplasm. A hallmark of cellular Ca2+ signalling is its spatiotemporal complexity: stimulation of cells by a hormone or neurotransmitter leads to oscillations in cytoplasmic Ca2+ concentration that can vary markedly in time course, amplitude, frequency, and spatial range. In this chapter we review some of the biological roles of Ca2+, the experimental characterisation of complex dynamic changes in Ca2+ concentration, and attempts to explain this complexity using computational models. We consider the 'toolkit' of cellular proteins which influence Ca2+ concentrarion, describe mechanistic models of key elements of the toolkit, and fit these into the framework of whole cell models of Ca2+ oscillations and waves. Finally, we will touch on recent efforts to use stochastic modelling to elucidate elementary Ca2+ signal events, and how these may evolve into global signals.
Collapse
Affiliation(s)
- Ruediger Thul
- School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | | | | | | | | |
Collapse
|
24
|
Abstract
The inositol 1,4,5-trisphosphate (InsP3) receptors (InsP3Rs) are a family of Ca2+ release channels localized predominately in the endoplasmic reticulum of all cell types. They function to release Ca2+ into the cytoplasm in response to InsP3 produced by diverse stimuli, generating complex local and global Ca2+ signals that regulate numerous cell physiological processes ranging from gene transcription to secretion to learning and memory. The InsP3R is a calcium-selective cation channel whose gating is regulated not only by InsP3, but by other ligands as well, in particular cytoplasmic Ca2+. Over the last decade, detailed quantitative studies of InsP3R channel function and its regulation by ligands and interacting proteins have provided new insights into a remarkable richness of channel regulation and of the structural aspects that underlie signal transduction and permeation. Here, we focus on these developments and review and synthesize the literature regarding the structure and single-channel properties of the InsP3R.
Collapse
Affiliation(s)
- J Kevin Foskett
- Department of Physiology, University of Pennsylvania, Philadelphia 19104-6085, USA.
| | | | | | | |
Collapse
|
25
|
Ionescu L, Cheung KH, Vais H, Mak DOD, White C, Foskett JK. Graded recruitment and inactivation of single InsP3 receptor Ca2+-release channels: implications for quantal [corrected] Ca2+release. J Physiol 2006; 573:645-62. [PMID: 16644799 PMCID: PMC1779751 DOI: 10.1113/jphysiol.2006.109504] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Modulation of cytoplasmic free Ca2+ concentration ([Ca2+]i) by receptor-mediated generation of inositol 1,4,5-trisphosphate (InsP3) and activation of its receptor (InsP3R), a Ca2+-release channel in the endoplasmic reticulum, is a ubiquitous signalling mechanism. A fundamental aspect of InsP3-mediated signalling is the graded release of Ca2+ in response to incremental levels of stimuli. Ca2+ release has a transient fast phase, whose rate is proportional to [InsP3], followed by a much slower one even in constant [InsP3]. Many schemes have been proposed to account for quantal Ca2+ release, including the presence of heterogeneous channels and Ca2+ stores with various mechanisms of release termination. Here, we demonstrate that mechanisms intrinsic to the single InsP3R channel can account for quantal Ca2+ release. Patch-clamp electrophysiology of isolated insect Sf9 cell nuclei revealed a consistent and high probability of detecting functional endogenous InsP3R channels, enabling InsP3-induced channel inactivation to be identified as an inevitable consequence of activation, and allowing the average number of activated channels in the membrane patch (N(A)) to be accurately quantified. InsP3-activated channels invariably inactivated, with average duration of channel activity reduced by high [Ca2+]i and suboptimal [InsP3]. Unexpectedly, N(A) was found to be a graded function of both [Ca2+]i and [InsP3]. A qualitative model involving Ca2+-induced InsP3R sequestration and inactivation can account for these observations. These results suggest that apparent heterogeneous ligand sensitivity can be generated in a homogeneous population of InsP3R channels, providing a mechanism for graded Ca2+ release that is intrinsic to the InsP3R Ca2+ release channel itself.
Collapse
Affiliation(s)
- Lucian Ionescu
- Department of Physiology, B39 Anatomy-Chemistry Building, 414 Guardian Drive, University of Pennsylvania, Philadelphia, PA 19104-6085, USA.
| | | | | | | | | | | |
Collapse
|
26
|
Sneyd J, Falcke M. Models of the inositol trisphosphate receptor. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2005; 89:207-45. [PMID: 15950055 DOI: 10.1016/j.pbiomolbio.2004.11.001] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
The inositol (1,4,5)-trisphosphate receptor (IPR) plays a crucial role in calcium dynamics in a wide range of cell types, and is often a central feature in quantitative models of calcium oscillations and waves. We review deterministic and stochastic mathematical models of the IPR, from the earliest ones of the 1970s and 1980s, to the most recent. The effects of IPR stochasticity on Ca2+ dynamics are briefly discussed.
Collapse
Affiliation(s)
- J Sneyd
- Department of Mathematics, University of Auckland, Auckland, New Zealand.
| | | |
Collapse
|
27
|
Roux E, Noble PJ, Noble D, Marhl M. Modelling of calcium handling in airway myocytes. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2005; 90:64-87. [PMID: 15982722 DOI: 10.1016/j.pbiomolbio.2005.05.004] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Airway myocytes are the primary effectors of airway reactivity which modulates airway resistance and hence ventilation. Stimulation of airway myocytes results in an increase in the cytosolic Ca(2+) concentration ([Ca(2+)](i)) and the subsequent activation of the contractile apparatus. Many contractile agonists, including acetylcholine, induce [Ca(2+)](i) increase via Ca(2+) release from the sarcoplasmic reticulum through InsP(3) receptors. Several models have been developed to explain the characteristics of InsP(3)-induced [Ca(2+)](i) responses, in particular Ca(2+) oscillations. The article reviews the modelling of the major structures implicated in intracellular Ca(2+) handling, i.e., InsP(3) receptors, SERCAs, mitochondria and Ca(2+)-binding cytosolic proteins. We developed theoretical models specifically dedicated to the airway myocyte which include the major mechanisms responsible for intracellular Ca(2+) handling identified in these cells. These biocomputations pointed out the importance of the relative proportion of InsP(3) receptor isoforms and the respective role of the different mechanisms responsible for cytosolic Ca(2+) clearance in the pattern of [Ca(2+)](i) variations. We have developed a theoretical model of membrane conductances that predicts the variations in membrane potential and extracellular Ca(2+) influx. Stimulation of this model by simulated increase in [Ca(2+)](i) predicts membrane depolarisation, but not great enough to trigger a significant opening of voltage-dependant Ca(2+) channels. This may explain why airway contraction induced by cholinergic stimulation does not greatly depend on extracellular calcium. The development of such models of airway myocytes is important for the understanding of the cellular mechanisms of airway reactivity and their possible modulation by pharmacological agents.
Collapse
Affiliation(s)
- Etienne Roux
- Laboratoire de Physiologie Cellulaire Respiratoire, INSERM E 356, Université Victor Segalen Bordeaux 2, 146 rue Léo-Saignat, 33076 Bordeaux cedex, France.
| | | | | | | |
Collapse
|