1
|
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
|
2
|
Falcke M, Moein M, Tilūnaitė A, Thul R, Skupin A. On the phase space structure of IP 3 induced Ca 2+ signalling and concepts for predictive modeling. CHAOS (WOODBURY, N.Y.) 2018; 28:045115. [PMID: 31906671 DOI: 10.1063/1.5021073] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The correspondence between mathematical structures and experimental systems is the basis of the generalizability of results found with specific systems and is the basis of the predictive power of theoretical physics. While physicists have confidence in this correspondence, it is less recognized in cellular biophysics. On the one hand, the complex organization of cellular dynamics involving a plethora of interacting molecules and the basic observation of cell variability seem to question its possibility. The practical difficulties of deriving the equations describing cellular behaviour from first principles support these doubts. On the other hand, ignoring such a correspondence would severely limit the possibility of predictive quantitative theory in biophysics. Additionally, the existence of functional modules (like pathways) across cell types suggests also the existence of mathematical structures with comparable universality. Only a few cellular systems have been sufficiently investigated in a variety of cell types to follow up these basic questions. IP3 induced Ca2+signalling is one of them, and the mathematical structure corresponding to it is subject of ongoing discussion. We review the system's general properties observed in a variety of cell types. They are captured by a reaction diffusion system. We discuss the phase space structure of its local dynamics. The spiking regime corresponds to noisy excitability. Models focussing on different aspects can be derived starting from this phase space structure. We discuss how the initial assumptions on the set of stochastic variables and phase space structure shape the predictions of parameter dependencies of the mathematical models resulting from the derivation.
Collapse
Affiliation(s)
- Martin Falcke
- Max Delbrück Centre for Molecular Medicine, Robert Rössler Strasse 10, 13125 Berlin, Germany and Department of Physics, Humboldt University, Newtonstr. 15, 12489 Berlin, Germany
| | - Mahsa Moein
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 7, Rue de Swing, Belval L-4367, Luxembourg
| | - Agne Tilūnaitė
- Systems Biology Laboratory, School of Mathematics and Statistics, University of Melbourne, Parkville, Victoria 3010, Australia
| | - Rüdiger Thul
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Alexander Skupin
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 7, Rue de Swing, Belval L-4367, Luxembourg
| |
Collapse
|
3
|
Gernert R, Emary C, Klapp SHL. Waiting time distribution for continuous stochastic systems. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:062115. [PMID: 25615052 DOI: 10.1103/physreve.90.062115] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/09/2014] [Indexed: 06/04/2023]
Abstract
The waiting time distribution (WTD) is a common tool for analyzing discrete stochastic processes in classical and quantum systems. However, there are many physical examples where the dynamics is continuous and only approximately discrete, or where it is favourable to discuss the dynamics on a discretized and a continuous level in parallel. An example is the hindered motion of particles through potential landscapes with barriers. In the present paper we propose a consistent generalization of the WTD from the discrete case to situations where the particles perform continuous barrier crossing characterized by a finite duration. To this end, we introduce a recipe to calculate the WTD from the Fokker-Planck (Smoluchowski) equation. In contrast to the closely related first passage time distribution (FPTD), which is frequently used to describe continuous processes, the WTD contains information about the direction of motion. As an application, we consider the paradigmatic example of an overdamped particle diffusing through a washboard potential. To verify the approach and to elucidate its numerical implications, we compare the WTD defined via the Smoluchowski equation with data from direct simulation of the underlying Langevin equation and find full consistency provided that the jumps in the Langevin approach are defined properly. Moreover, for sufficiently large energy barriers, the WTD defined via the Smoluchowski equation becomes consistent with that resulting from the analytical solution of a (two-state) master equation model for the short-time dynamics developed previously by us [Phys. Rev. E 86, 061135 (2012)]. Thus, our approach "interpolates" between these two types of stochastic motion. We illustrate our approach for both symmetric systems and systems under constant force.
Collapse
Affiliation(s)
- Robert Gernert
- Institut für Theoretische Physik, Sekr. EW 7-1, Technische Universität Berlin, Hardenbergstrasse 36, D-10623 Berlin, Germany
| | - Clive Emary
- Department of Physics and Mathematics, University of Hull, Kingston-upon-Hull, HU6 7RX, United Kingdom
| | - Sabine H L Klapp
- Institut für Theoretische Physik, Sekr. EW 7-1, Technische Universität Berlin, Hardenbergstrasse 36, D-10623 Berlin, Germany
| |
Collapse
|
4
|
Thul R. Time to blip--stochastic simulation of single channel opening. Cold Spring Harb Protoc 2014; 2014:2014/5/pdb.prot073239. [PMID: 24786499 DOI: 10.1101/pdb.prot073239] [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/25/2022]
Abstract
The stochastic dynamics of the inositol-1,4,5-trisphosphate (IP3) receptor (IP3R) is key to understanding a wide range of observed calcium (Ca²⁺) signals. The stochastic nature results from the constant binding and unbinding of Ca²⁺ and IP3 to and from their respective binding sites and is especially important in the initiation of a Ca²⁺ puff (i.e., the release of Ca²⁺ through a cluster of IP3Rs). Once the first IP3R opens, the Ca²⁺ concentration rises significantly around the ion channel, increasing the open probability for neighboring IP3Rs. This opening may trigger the activation of further receptors, giving rise to a Ca²⁺ puff. In this protocol, we determine the time that it takes for a single IP3R to open from rest. We explicitly take into account the tetrameric structure of the IP3R and the fact that multiple subunits must be active before the channel opens. We develop code for a stochastic simulation of the IP3R and simulate it using the software package MATLAB. This protocol shows the basic form of a stochastic simulation algorithm and may serve as a starting point to investigate more complex gating dynamics.
Collapse
Affiliation(s)
- Rüdiger Thul
- School of Mathematical Sciences, University of Nottingham, Nottingham, NG7 2RD, United Kingdom
| |
Collapse
|
5
|
Christian N, Skupin A, Morante S, Jansen K, Rossi G, Ebenhöh O. Mesoscopic behavior from microscopic Markov dynamics and its application to calcium release channels. J Theor Biol 2013; 343:102-12. [PMID: 24270093 DOI: 10.1016/j.jtbi.2013.11.010] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2013] [Revised: 09/03/2013] [Accepted: 11/13/2013] [Indexed: 11/19/2022]
Abstract
A major challenge in biology is to understand how molecular processes determine phenotypic features. We address this fundamental problem in a class of model systems by developing a general mathematical framework that allows the calculation of mesoscopic properties from the knowledge of microscopic Markovian transition probabilities. We show how exact analytic formulae for the first and second moments of resident time distributions in mesostates can be derived from microscopic resident times and transition probabilities even for systems with a large number of microstates. We apply our formalism to models of the inositol trisphosphate receptor, which plays a key role in generating calcium signals triggering a wide variety of cellular responses. We demonstrate how experimentally accessible quantities, such as opening and closing times and the coefficient of variation of inter-spike intervals, and other, more elaborated, quantities can be analytically calculated from the underlying microscopic Markovian dynamics. A virtue of our approach is that we do not need to follow the detailed time evolution of the whole system, as we derive the relevant properties of its steady state without having to take into account the often extremely complicated transient features. We emphasize that our formulae fully agree with results obtained by stochastic simulations and approaches based on a full determination of the microscopic system's time evolution. We also illustrate how experiments can be devised to discriminate between alternative molecular models of the inositol trisphosphate receptor. The developed approach is applicable to any system described by a Markov process and, owing to the analytic nature of the resulting formulae, provides an easy way to characterize also rare events that are of particular importance to understand the intermittency properties of complex dynamic systems.
Collapse
Affiliation(s)
- Nils Christian
- University of Aberdeen, Department of Physics, Meston Walk, Aberdeen AB24 3UE, UK; University Luxembourg, Luxembourg Centre for Systems Biomedicine, 7, Avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg.
| | - Alexander Skupin
- University Luxembourg, Luxembourg Centre for Systems Biomedicine, 7, Avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg; Institute for Systems Biology, 401 Terry Ave N, Seattle, WA 98109, USA
| | - Silvia Morante
- Università di Roma Tor Vergata and INFN, Sezione di Roma 2, Via della Ricerca Scientifica, I-00133 Roma, Italy
| | - Karl Jansen
- NIC/DESY Zeuthen, Platanenallee 6, D-15738 Zeuthen, Germany
| | - Giancarlo Rossi
- Università di Roma Tor Vergata and INFN, Sezione di Roma 2, Via della Ricerca Scientifica, I-00133 Roma, Italy
| | - Oliver Ebenhöh
- University of Aberdeen, Department of Physics, Meston Walk, Aberdeen AB24 3UE, UK; Helmholtz Centre Potsdam, GFZ German Research Centre for Geosciences, 14473 Potsdam, Germany; Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich-Heine-University, Universitätsstraße 1, D-40225 Düsseldorf, Germany.
| |
Collapse
|
6
|
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
|
7
|
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
|
8
|
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
|
9
|
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
|
10
|
Derivation of Ca2+ signals from puff properties reveals that pathway function is robust against cell variability but sensitive for control. Proc Natl Acad Sci U S A 2010; 108:427-32. [PMID: 21173273 DOI: 10.1073/pnas.1008435108] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Ca(2+) is a universal second messenger in eukaryotic cells transmitting information through sequences of concentration spikes. A prominent mechanism to generate these spikes involves Ca(2+) release from the endoplasmic reticulum Ca(2+) store via inositol 1,4,5-trisphosphate (IP(3))-sensitive channels. Puffs are elemental events of IP(3)-induced Ca(2+) release through single clusters of channels. Intracellular Ca(2+) dynamics are a stochastic system, but a complete stochastic theory has not been developed yet. We formulate the theory in terms of interpuff interval and puff duration distributions because, unlike the properties of individual channels, they can be measured in vivo. Our theory reproduces the typical spectrum of Ca(2+) signals like puffs, spiking, and bursting in analytically treatable test cases as well as in more realistic simulations. We find conditions for spiking and calculate interspike interval (ISI) distributions. Signal form, average ISI and ISI distributions depend sensitively on the details of cluster properties and their spatial arrangement. In contrast to that, the relation between the average and the standard deviation of ISIs does not depend on cluster properties and cluster arrangement and is robust with respect to cell variability. It is controlled by the global feedback processes in the Ca(2+) signaling pathway (e.g., via IP(3)-3-kinase or endoplasmic reticulum depletion). That relation is essential for pathway function because it ensures frequency encoding despite the randomness of ISIs and determines the maximal spike train information content. Hence, we find a division of tasks between global feedbacks and local cluster properties that guarantees robustness of function while maintaining sensitivity of control of the average ISI.
Collapse
|
11
|
Tania N, Keener JP. Calsequestrin mediates changes in spontaneous calcium release profiles. J Theor Biol 2010; 265:359-76. [PMID: 20648970 DOI: 10.1016/j.jtbi.2010.05.028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
Calsequestrin (CSQ) is the primary calcium buffer within the sarcoplasmic reticulum (SR) of cardiac cells. It has also been identified as a regulator of Ryanodine receptor (RyR) calcium release channels by serving as a SR luminal sensor. When calsequestrin is free and unbound to calcium, it can bind to RyR and desensitize the channel from cytoplasmic calcium activation. In this paper, we study the role of CSQ as a buffer and RyR luminal sensor using a mechanistic model of RyR-CSQ interaction. By using various asymptotic approximations and mean first exit time calculation, we derive a minimal model of a calcium release unit which includes CSQ dependence. Using this model, we then analyze the effect of changing CSQ expression on the calcium release profile and the rate of spontaneous calcium release. We show that because of its buffering capability, increasing CSQ increases the spark duration and size. However, because of luminal sensing effects, increasing CSQ depresses the basal spark rate and increases the critical SR level for calcium release termination. Finally, we show that with increased bulk cytoplasmic calcium concentration, the CRU model exhibits deterministic oscillations.
Collapse
Affiliation(s)
- Nessy Tania
- Department of Mathematics, University of Utah, 155 S. 1400 E. Room 233, Salt Lake City, UT 84112, USA.
| | | |
Collapse
|