1
|
Mishra V, Adlakha N. Cross Talking Calcium, IP 3 and Buffer Dynamics Alters ATP and NADH Level in Obese and Normal Hepatocyte Cell. Cell Biochem Biophys 2024:10.1007/s12013-024-01306-9. [PMID: 38789660 DOI: 10.1007/s12013-024-01306-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/08/2024] [Indexed: 05/26/2024]
Abstract
The cross talk between calcium (Ca2+), IP3 and buffer dynamics regulate various mechanisms in hepatocyte cells. The study of independent systems of calcium, IP3, and buffer signaling provides limited information about cell dynamics. In the current study, coupled reaction-diffusion equations are used to design a cross-talk model for IP3, buffer, and calcium dynamics in a hepatocyte cell. The one-way feedback of calcium, buffer, and IP3 in ATP production, ATP degradation, and NADH production rate is incorporated into the model. Numerical simulation has been done using the Finite Element Method (FEM) along the spatial direction and the Crank-Nicolson (C-N) method along the temporal direction. The numerical results are analysed to determine the effects of alterations in processes of cross-talking dynamics of IP3, buffer, and calcium on ATP and NADH production and degradation rate of ATP in a hepatocyte cell under normal and obesity conditions. The comparative analysis of these findings unveils notable distinctions induced by obesity in calcium dynamics, ATP and NADH synthesis, and ATP degradation kinetics.
Collapse
Affiliation(s)
- Vedika Mishra
- Department of Mathematics, SVNIT, Surat, 395007, Gujarat, India.
| | - Neeru Adlakha
- Department of Mathematics, SVNIT, Surat, 395007, Gujarat, India
| |
Collapse
|
2
|
Gao H, Liu L, Zaikin A, Chen S. An efficient reduced-lattice model of IP 3R for probing Ca 2+ dynamics. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2023; 1865:184195. [PMID: 37353068 DOI: 10.1016/j.bbamem.2023.184195] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Revised: 06/15/2023] [Accepted: 06/16/2023] [Indexed: 06/25/2023]
Abstract
Numerous cellular processes are regulated by Ca2+ signals, and the endoplasmic reticulum (ER) membrane's inositol triphosphate receptor (IP3R) is critical for modulating intracellular Ca2+ dynamics. The IP3Rs are seen to be clustered in a variety of cell types. The combination of IP3Rs clustering and IP3Rs-mediated Ca2+-induced Ca2+ release results in the hierarchical organization of the Ca2+ signals, which challenges the numerical simulation given the multiple spatial and temporal scales that must be covered. The previous methods rather ignore the spatial feature of IP3Rs or fail to coordinate the conflicts between the real biological relevance and the computational cost. In this work, a general and efficient reduced-lattice model is presented for the simulation of IP3Rs-mediated multiscale Ca2+ dynamics. The model highlights biological details that make up the majority of the calcium events, including IP3Rs clustering and calcium domains, and it reduces the complexity by approximating the minor details. The model's extensibility provides fresh insights into the function of IP3Rs in producing global Ca2+ events and supports the research under more physiological circumstances. Our work contributes to a novel toolkit for modeling multiscale Ca2+ dynamics and advances knowledge of Ca2+ signals.
Collapse
Affiliation(s)
- Huayi Gao
- Britton Chance Center for Biomedical Photonics, Wuhan National Laboratory for Optoelectronics-Huazhong University of Science and Technology, Wuhan, China; MoE Key Laboratory for Biomedical Photonics, School of Engineering Sciences, Huazhong University of Science and Technology, Wuhan, China
| | - Langzhou Liu
- Britton Chance Center for Biomedical Photonics, Wuhan National Laboratory for Optoelectronics-Huazhong University of Science and Technology, Wuhan, China; MoE Key Laboratory for Biomedical Photonics, School of Engineering Sciences, Huazhong University of Science and Technology, Wuhan, China
| | - Alexey Zaikin
- Institute for Women's Health and Department of Mathematics, University College London, Gower St., London, WC1E 6BT, UK; Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Ave., 7 Nizhniy Novgorod, 603022, Russia; Centre for Analysis of Complex Systems, Sechenov First Moscow State Medical University, Moscow, 119019, Russia
| | - Shangbin Chen
- Britton Chance Center for Biomedical Photonics, Wuhan National Laboratory for Optoelectronics-Huazhong University of Science and Technology, Wuhan, China; MoE Key Laboratory for Biomedical Photonics, School of Engineering Sciences, Huazhong University of Science and Technology, Wuhan, China.
| |
Collapse
|
3
|
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
|
4
|
Takano T, Wahl AM, Huang KT, Narita T, Rugis J, Sneyd J, Yule DI. Highly localized intracellular Ca 2+ signals promote optimal salivary gland fluid secretion. eLife 2021; 10:66170. [PMID: 34240705 PMCID: PMC8352588 DOI: 10.7554/elife.66170] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2020] [Accepted: 07/08/2021] [Indexed: 11/23/2022] Open
Abstract
Salivary fluid secretion involves an intricate choreography of membrane transporters to result in the trans-epithelial movement of NaCl and water into the acinus lumen. Current models are largely based on experimental observations in enzymatically isolated cells where the Ca2+ signal invariably propagates globally and thus appears ideally suited to activate spatially separated Cl and K channels, present on the apical and basolateral plasma membrane, respectively. We monitored Ca2+ signals and salivary secretion in live mice expressing GCamp6F, following stimulation of the nerves innervating the submandibular gland. Consistent with in vitro studies, Ca2+ signals were initiated in the apical endoplasmic reticulum. In marked contrast to in vitro data, highly localized trains of Ca2+ transients that failed to fully propagate from the apical region were observed. Following stimuli optimum for secretion, large apical-basal gradients were elicited. A new mathematical model, incorporating these data was constructed to probe how salivary secretion can be optimally stimulated by apical Ca2+ signals.
Collapse
Affiliation(s)
- Takahiro Takano
- Department of Pharmacology and Physiology, University of Rochester, Rochester, United States
| | - Amanda M Wahl
- Department of Pharmacology and Physiology, University of Rochester, Rochester, United States
| | - Kai-Ting Huang
- Department of Pharmacology and Physiology, University of Rochester, Rochester, United States
| | - Takanori Narita
- Department of Veterinary Medicine, College of Bioresource Sciences, Nihon University, Fujisawa, Japan
| | - John Rugis
- Department of Mathematics, University of Auckland, Auckland, New Zealand
| | - James Sneyd
- Department of Mathematics, University of Auckland, Auckland, New Zealand
| | - David I Yule
- Department of Pharmacology and Physiology, University of Rochester, Rochester, United States
| |
Collapse
|
5
|
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
|
6
|
Cloete I, Corrêa-Velloso JC, Bartlett PJ, Kirk V, Thomas AP, Sneyd J. A Tale of two receptors. J Theor Biol 2021; 518:110629. [PMID: 33607144 DOI: 10.1016/j.jtbi.2021.110629] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2020] [Revised: 01/10/2021] [Accepted: 02/05/2021] [Indexed: 11/26/2022]
Abstract
Calcium (Ca2+) oscillations in hepatocytes have a wide dynamic range. In particular, recent experimental evidence shows that agonist stimulation of the P2Y family of receptors leads to qualitatively diverse Ca2+ oscillations. We present a new model of Ca2+ oscillations in hepatocytes based on these experiments to investigate the mechanisms controlling P2Y-activated Ca2+ oscillations. The model accounts for Ca2+ regulation of the IP3 receptor (IP3R), the positive feedback from Ca2+ on phospholipase C (PLC) and the P2Y receptor phosphorylation by protein kinase C (PKC). Furthermore, PKC is shown to control multiple cellular substrates. Utilising the model, we suggest the activity and intensity of PLC and PKC necessary to explain the qualitatively diverse Ca2+ oscillations in response to P2Y receptor activation.
Collapse
Affiliation(s)
- Ielyaas Cloete
- Department of Mathematics, University of Auckland, Auckland 1142, New Zealand
| | - Juliana C Corrêa-Velloso
- Department of Pharmacology, Physiology and Neuroscience, New Jersey Medical School Rutgers, The State University of New Jersey, Newark, NJ 07103, United States
| | - Paula J Bartlett
- Department of Pharmacology, Physiology and Neuroscience, New Jersey Medical School Rutgers, The State University of New Jersey, Newark, NJ 07103, United States
| | - Vivien Kirk
- Department of Mathematics, University of Auckland, Auckland 1142, New Zealand
| | - Andrew P Thomas
- Department of Pharmacology, Physiology and Neuroscience, New Jersey Medical School Rutgers, The State University of New Jersey, Newark, NJ 07103, United States
| | - James Sneyd
- Department of Mathematics, University of Auckland, Auckland 1142, New Zealand
| |
Collapse
|
7
|
Dual mechanisms of Ca2+ oscillations in hepatocytes. J Theor Biol 2020; 503:110390. [DOI: 10.1016/j.jtbi.2020.110390] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2020] [Revised: 06/19/2020] [Accepted: 06/24/2020] [Indexed: 11/30/2022]
|
8
|
Powell J, Falcke M, Skupin A, Bellamy TC, Kypraios T, Thul R. A Statistical View on Calcium Oscillations. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2020; 1131:799-826. [PMID: 31646535 DOI: 10.1007/978-3-030-12457-1_32] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Transient rises and falls of the intracellular calcium concentration have been observed in numerous cell types and under a plethora of conditions. There is now a growing body of evidence that these whole-cell calcium oscillations are stochastic, which poses a significant challenge for modelling. In this review, we take a closer look at recently developed statistical approaches to calcium oscillations. These models describe the timing of whole-cell calcium spikes, yet their parametrisations reflect subcellular processes. We show how non-stationary calcium spike sequences, which e.g. occur during slow depletion of intracellular calcium stores or in the presence of time-dependent stimulation, can be analysed with the help of so-called intensity functions. By utilising Bayesian concepts, we demonstrate how values of key parameters of the statistical model can be inferred from single cell calcium spike sequences and illustrate what information whole-cell statistical models can provide about the subcellular mechanistic processes that drive calcium oscillations. In particular, we find that the interspike interval distribution of HEK293 cells under constant stimulation is captured by a Gamma distribution.
Collapse
Affiliation(s)
- Jake Powell
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Martin Falcke
- Max Delbrück Centre for Molecular Medicine, Berlin, Germany.,Department of Physics, Humboldt University, Berlin, Germany
| | - Alexander Skupin
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belval, Luxembourg.,National Biomedical Computation Resource, University California San Diego, La Jolla, CA, USA
| | - Tomas C Bellamy
- School of Life Sciences, University of Nottingham, Nottingham, UK
| | - Theodore Kypraios
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Rüdiger Thul
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK.
| |
Collapse
|
9
|
Voorsluijs V, Dawson SP, De Decker Y, Dupont G. Deterministic Limit of Intracellular Calcium Spikes. PHYSICAL REVIEW LETTERS 2019; 122:088101. [PMID: 30932600 DOI: 10.1103/physrevlett.122.088101] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/10/2018] [Revised: 01/09/2019] [Indexed: 06/09/2023]
Abstract
In nonexcitable cells, global Ca^{2+} spikes emerge from the collective dynamics of clusters of Ca^{2+} channels that are coupled by diffusion. Current modeling approaches have opposed stochastic descriptions of these systems to purely deterministic models, while both paradoxically appear compatible with experimental data. Combining fully stochastic simulations and mean-field analyses, we demonstrate that these two approaches can be reconciled. Our fully stochastic model generates spike sequences that can be seen as noise-perturbed oscillations of deterministic origin, while displaying statistical properties in agreement with experimental data. These underlying deterministic oscillations arise from a phenomenological spike nucleation mechanism.
Collapse
Affiliation(s)
- V Voorsluijs
- Nonlinear Physical Chemistry Unit and Center for Nonlinear Phenomena and Complex Systems (CENOLI), Université libre de Bruxelles, Boulevard du Triomphe, C.P. 231, B-1050 Brussels, Belgium
| | - S Ponce Dawson
- Departamento de Física, FCEN-UBA and IFIBA, UBA-CONICET, Ciudad Universitaria, Pabellón I, 1428 Buenos Aires, Argentina
| | - Y De Decker
- Nonlinear Physical Chemistry Unit and Center for Nonlinear Phenomena and Complex Systems (CENOLI), Université libre de Bruxelles, Boulevard du Triomphe, C.P. 231, B-1050 Brussels, Belgium
| | - G Dupont
- Unité de Chronobiologie Théorique, Université libre de Bruxelles, Boulevard du Triomphe, C.P. 231, B-1050 Brussels, Belgium
| |
Collapse
|
10
|
Savtchenko LP, Bard L, Jensen TP, Reynolds JP, Kraev I, Medvedev N, Stewart MG, Henneberger C, Rusakov DA. Disentangling astroglial physiology with a realistic cell model in silico. Nat Commun 2018; 9:3554. [PMID: 30177844 PMCID: PMC6120909 DOI: 10.1038/s41467-018-05896-w] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2018] [Accepted: 07/10/2018] [Indexed: 12/05/2022] Open
Abstract
Electrically non-excitable astroglia take up neurotransmitters, buffer extracellular K+ and generate Ca2+ signals that release molecular regulators of neural circuitry. The underlying machinery remains enigmatic, mainly because the sponge-like astrocyte morphology has been difficult to access experimentally or explore theoretically. Here, we systematically incorporate multi-scale, tri-dimensional astroglial architecture into a realistic multi-compartmental cell model, which we constrain by empirical tests and integrate into the NEURON computational biophysical environment. This approach is implemented as a flexible astrocyte-model builder ASTRO. As a proof-of-concept, we explore an in silico astrocyte to evaluate basic cell physiology features inaccessible experimentally. Our simulations suggest that currents generated by glutamate transporters or K+ channels have negligible distant effects on membrane voltage and that individual astrocytes can successfully handle extracellular K+ hotspots. We show how intracellular Ca2+ buffers affect Ca2+ waves and why the classical Ca2+ sparks-and-puffs mechanism is theoretically compatible with common readouts of astroglial Ca2+ imaging.
Collapse
Affiliation(s)
- Leonid P Savtchenko
- UCL Institute of Neurology, University College London, London, WC1N 3BG, UK.
| | - Lucie Bard
- UCL Institute of Neurology, University College London, London, WC1N 3BG, UK
| | - Thomas P Jensen
- UCL Institute of Neurology, University College London, London, WC1N 3BG, UK
| | - James P Reynolds
- UCL Institute of Neurology, University College London, London, WC1N 3BG, UK
| | - Igor Kraev
- The Open University, Milton Keynes, MK7 6AA, UK
| | | | | | - Christian Henneberger
- UCL Institute of Neurology, University College London, London, WC1N 3BG, UK
- German Center of Neurodegenerative Diseases (DZNE), Bonn, 53127, Germany
- Institute of Cellular Neurosciences, University of Bonn Medical School, Bonn, 53127, Germany
| | - Dmitri A Rusakov
- UCL Institute of Neurology, University College London, London, WC1N 3BG, UK.
| |
Collapse
|
11
|
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
|
12
|
Thurley K, Wu LF, Altschuler SJ. Modeling Cell-to-Cell Communication Networks Using Response-Time Distributions. Cell Syst 2018; 6:355-367.e5. [PMID: 29525203 PMCID: PMC5913757 DOI: 10.1016/j.cels.2018.01.016] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2017] [Revised: 10/10/2017] [Accepted: 01/26/2018] [Indexed: 01/30/2023]
Abstract
Cell-to-cell communication networks have critical roles in coordinating diverse organismal processes, such as tissue development or immune cell response. However, compared with intracellular signal transduction networks, the function and engineering principles of cell-to-cell communication networks are far less understood. Major complications include: cells are themselves regulated by complex intracellular signaling networks; individual cells are heterogeneous; and output of any one cell can recursively become an additional input signal to other cells. Here, we make use of a framework that treats intracellular signal transduction networks as "black boxes" with characterized input-to-output response relationships. We study simple cell-to-cell communication circuit motifs and find conditions that generate bimodal responses in time, as well as mechanisms for independently controlling synchronization and delay of cell-population responses. We apply our modeling approach to explain otherwise puzzling data on cytokine secretion onset times in T cells. Our approach can be used to predict communication network structure using experimentally accessible input-to-output measurements and without detailed knowledge of intermediate steps.
Collapse
Affiliation(s)
- Kevin Thurley
- Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, CA 94158, USA,Correspondence: (K.T.), (L.F.W.), (S.J.A.)
| | - Lani F. Wu
- Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, CA 94158, USA,Correspondence: (K.T.), (L.F.W.), (S.J.A.)
| | - Steven J. Altschuler
- Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, CA 94158, USA,Correspondence: (K.T.), (L.F.W.), (S.J.A.)
| |
Collapse
|
13
|
Tilūnaitė A, Croft W, Russell N, Bellamy TC, Thul R. A Bayesian approach to modelling heterogeneous calcium responses in cell populations. PLoS Comput Biol 2017; 13:e1005794. [PMID: 28985235 PMCID: PMC5646906 DOI: 10.1371/journal.pcbi.1005794] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2017] [Revised: 10/18/2017] [Accepted: 09/27/2017] [Indexed: 12/18/2022] Open
Abstract
Calcium responses have been observed as spikes of the whole-cell calcium concentration in numerous cell types and are essential for translating extracellular stimuli into cellular responses. While there are several suggestions for how this encoding is achieved, we still lack a comprehensive theory. To achieve this goal it is necessary to reliably predict the temporal evolution of calcium spike sequences for a given stimulus. Here, we propose a modelling framework that allows us to quantitatively describe the timing of calcium spikes. Using a Bayesian approach, we show that Gaussian processes model calcium spike rates with high fidelity and perform better than standard tools such as peri-stimulus time histograms and kernel smoothing. We employ our modelling concept to analyse calcium spike sequences from dynamically-stimulated HEK293T cells. Under these conditions, different cells often experience diverse stimulus time courses, which is a situation likely to occur in vivo. This single cell variability and the concomitant small number of calcium spikes per cell pose a significant modelling challenge, but we demonstrate that Gaussian processes can successfully describe calcium spike rates in these circumstances. Our results therefore pave the way towards a statistical description of heterogeneous calcium oscillations in a dynamic environment. Upon stimulation a large number of cell types respond with transient increases of the intracellular calcium concentration, which often take the form of repetitive spikes. It is therefore believed that calcium spikes play a central role in cellular signal transduction. A critical feature of these calcium spikes is that they occur randomly, which raises the question of how we can predict the timing of calcium spikes. We here show that by using Bayesian ideas and concepts from stochastic processes, we can quantitatively compute the calcium spike rate for a given stimulus. Our analysis also demonstrates that traditional methods for spike rate estimation perform less favourably compared to a Bayesian approach when small numbers of cells are investigated. To test our methodology under conditions that closely mimic those experienced in vivo we challenged cells with agonist concentrations that vary both in space and time. We find that cells that experience similar stimulus profiles are described by similar calcium spike rates. This suggests that calcium spike rates may constitute a quantitative description of whole-cell calcium spiking that reflects both the randomness and the spatiotemporal organisation of the calcium signalling machinery.
Collapse
Affiliation(s)
- Agne Tilūnaitė
- School of Mathematical Sciences, University of Nottingham, Nottingham, England, United Kingdom
| | - Wayne Croft
- School of Life Sciences, University of Nottingham, Nottingham, England, United Kingdom
| | - Noah Russell
- Department of Electrical and Electronic Engineering, University of Nottingham, Nottingham, England, United Kingdom
| | - Tomas C Bellamy
- School of Life Sciences, University of Nottingham, Nottingham, England, United Kingdom
| | - Rüdiger Thul
- School of Mathematical Sciences, University of Nottingham, Nottingham, England, United Kingdom
| |
Collapse
|
14
|
Cao P, Falcke M, Sneyd J. Mapping Interpuff Interval Distribution to the Properties of Inositol Trisphosphate Receptors. Biophys J 2017; 112:2138-2146. [PMID: 28538151 DOI: 10.1016/j.bpj.2017.03.019] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2017] [Revised: 03/14/2017] [Accepted: 03/24/2017] [Indexed: 01/24/2023] Open
Abstract
Tightly clustered inositol trisphosphate receptors (IP3Rs) control localized Ca2+ liberation from the endoplasmic reticulum to generate repetitive Ca2+ puffs. Distributions of the interpuff interval (IPI), i.e., the waiting time between successive puffs, are found to be well characterized by a probability density function involving only two parameters, λ and ξ, which represent the basal rate of puff generation and the recovery rate from refractoriness, respectively. However, how the two parameters depend on the kinetic parameters of single IP3Rs in a cluster is still unclear. In this article, using a stochastic puff model and a single-channel data-based IP3R model, we establish the dependencies of λ and ξ on two important IP3R model parameters, IP3 concentration ([IP3]) and the recovery rate from Ca2+ inhibition (rlow). By varying [IP3] and rlow in physiologically plausible ranges, we find that the ξ-λ plane is comprised of only two disjoint regions, a biologically impermissible region and a region where each parameter set (ξ, λ) can be caused by using two different combinations of [IP3] and rlow. The two combinations utilize very different mechanisms to maintain the same IPI distribution, and the mechanistic difference provides a way of identifying IP3R kinetic parameters by observing properties of the IPI.
Collapse
Affiliation(s)
- Pengxing Cao
- School of Mathematics and Statistics, The University of Melbourne, Melbourne, Victoria, Australia
| | - 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
|
15
|
|
16
|
Chen Y, Qi H, Li X, Cai M, Chen X, Liu W, Shuai J. Suppressing effect of Ca^{2+} blips on puff amplitudes by inhibiting channels to prevent recovery. Phys Rev E 2016; 94:022411. [PMID: 27627339 DOI: 10.1103/physreve.94.022411] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2016] [Indexed: 11/07/2022]
Abstract
As local signals, calcium puffs arise from the concerted opening of a few nearby inositol 1,4,5-trisphospate receptor channels to release Ca^{2+} ions from the endoplasmic reticulum. Although Ca^{2+} puffs have been well studied, little is known about the modulation of cytosolic basal Ca^{2+} concentration ([Ca^{2+}]_{Basal}) on puff dynamics. In this paper we consider a puff model to study how the statistical properties of puffs are modulated by [Ca^{2+}]_{Basal}. The puff frequency and lifetime trivially increase with the increasing [Ca^{2+}]_{Basal}, but an unexpected result is that the puff amplitude and the maximum open-channel number of the puff show decreasing relationship with the increasing [Ca^{2+}]_{Basal}. The underlying dynamics is related not only to the increasing puff frequency which gives a shorter recovery time, but also to the increasing frequency of blips with only one channel open. We indicate that Ca^{2+} blips cause the channels to be inhibited and prevent their recovery during interpuff intervals, resulting in the suppressing effect on puff amplitudes. With increasing [Ca^{2+}]_{Basal}, more blips occur to cause more channels to be inhibited, leaving fewer channels available for puff events. This study shows that the blips may play relevant functions in global Ca^{2+} waves through modulating puff dynamics.
Collapse
Affiliation(s)
- Yuan Chen
- Department of Physics, College of Physics Science and Technology, Xiamen University, Xiamen 361005, China
| | - Hong Qi
- Complex Systems Research Center, Shanxi University, Taiyuan 030006, China
| | - Xiang Li
- Department of Physics, College of Physics Science and Technology, Xiamen University, Xiamen 361005, China
| | - Meichun Cai
- Department of Physics, College of Physics Science and Technology, Xiamen University, Xiamen 361005, China
| | - Xingqiang Chen
- Department of Physics, College of Physics Science and Technology, Xiamen University, Xiamen 361005, China
| | - Wen Liu
- Department of Physics, College of Physics Science and Technology, Xiamen University, Xiamen 361005, China
| | - Jianwei Shuai
- Department of Physics, College of Physics Science and Technology, Xiamen University, Xiamen 361005, China.,State Key Laboratory of Cellular Stress Biology, Innovation Center for Cell Signaling Network, Xiamen University, Xiamen 361102, China
| |
Collapse
|
17
|
Damped oscillations in a multiple delayed feedback NF-κB signaling module. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2015; 44:677-84. [PMID: 26290058 DOI: 10.1007/s00249-015-1066-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2015] [Revised: 06/10/2015] [Accepted: 07/28/2015] [Indexed: 10/23/2022]
Abstract
The NF-[Formula: see text]B signaling system regulates a number of cellular processes. Recent studies with simplified models found a damped function of the dual delayed feedback NF-κB signaling module. We use a computational model to investigate how multiple delayed feedback aids achieving damping oscillation in the system and how internal noise can influence the damping function. A curve-fitting method (CFM) is introduced to quantify the damped oscillation. Our results show that (1) the structure of multiple delayed feedback, containing double or triple significantly delayed feedback, determines achieving damped oscillation. (2) Internal noise could aid the system to achieve damped oscillation under almost all circumstances.
Collapse
|
18
|
Ca2+-activation kinetics modulate successive puff/spark amplitude, duration and inter-event-interval correlations in a Langevin model of stochastic Ca2+ release. Math Biosci 2015; 264:101-7. [DOI: 10.1016/j.mbs.2015.03.012] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2014] [Revised: 03/16/2015] [Accepted: 03/27/2015] [Indexed: 01/12/2023]
|
19
|
Abstract
In a normal human life span, the heart beats about 2 to 3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.
Collapse
Affiliation(s)
- Zhilin Qu
- Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, California 90095, USA
- Correspondence to: Zhilin Qu, PhD, Department of Medicine, Division of Cardiology, David Geffen School of Medicine at UCLA, A2-237 CHS, 650 Charles E. Young Drive South, Los Angeles, CA 90095, Tel: 310-794-6050, Fax: 310-206-9133,
| | - Gang Hu
- Department of Physics, Beijing Normal University, Beijing 100875, China
| | - Alan Garfinkel
- Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, California 90095, USA
- Department of Integrative Biology and Physiology, University of California, Los Angeles, California 90095, USA
| | - James N. Weiss
- Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, California 90095, USA
- Department of Physiology, David Geffen School of Medicine, University of California, Los Angeles, California 90095, USA
| |
Collapse
|
20
|
Dickinson GD, Parker I. Factors determining the recruitment of inositol trisphosphate receptor channels during calcium puffs. Biophys J 2014; 105:2474-84. [PMID: 24314078 DOI: 10.1016/j.bpj.2013.10.028] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2013] [Accepted: 10/28/2013] [Indexed: 11/29/2022] Open
Abstract
Puffs are localized, transient elevations in cytosolic Ca(2+) that serve both as the building blocks of global cellular Ca(2+) signals and as local signals in their own right. They arise from clustered inositol 1,4,5-trisphosphate receptor/channels (IP3Rs), whose openings are coordinated by Ca(2+)-induced Ca(2+) release (CICR). We utilized total internal reflection fluorescence imaging of Ca(2+) signals in neuroblastoma cells with single-channel resolution to elucidate the mechanisms determining the triggering, amplitudes, kinetics, and spatial spread of puffs. We find that any given channel in a cluster has a mean probability of ∼66% of opening following opening of an initial "trigger" channel, and the probability of puff triggering thus increases steeply with increasing number of channels in a cluster (cluster size). Mean puff amplitudes scale with cluster size, but individual amplitudes vary widely, even at sites of similar cluster size, displaying similar proportions of events involving any given number of the channels in the cluster. Stochastic variation in numbers of Ca(2+)-inhibited IP3Rs likely contributes to the variability of amplitudes of repeated puffs at a site but the amplitudes of successive puffs were uncorrelated, even though we observed statistical correlations between interpuff intervals and puff amplitudes. Initial puffs evoked following photorelease of IP3-which would not be subject to earlier Ca(2+)-inhibition-also showed wide variability, indicating that mechanisms such as stochastic variation in IP3 binding and channel recruitment by CICR further determine puff amplitudes. The mean termination time of puffs lengthened with increasing puff amplitude size, consistent with independent closings of channels after a given mean open time, but we found no correlation of termination time with cluster size independent of puff amplitude. The spatial extent of puffs increased with their amplitude, and puffs of similar size were of similar width, independent of cluster size.
Collapse
Affiliation(s)
- George D Dickinson
- Department of Neurobiology and Behavior, University of California, Irvine, CA.
| | | |
Collapse
|
21
|
Ellefsen KL, Settle B, Parker I, Smith IF. An algorithm for automated detection, localization and measurement of local calcium signals from camera-based imaging. Cell Calcium 2014; 56:147-56. [PMID: 25047761 DOI: 10.1016/j.ceca.2014.06.003] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2014] [Revised: 06/10/2014] [Accepted: 06/15/2014] [Indexed: 10/25/2022]
Abstract
Local Ca(2+) transients such as puffs and sparks form the building blocks of cellular Ca(2+) signaling in numerous cell types. They have traditionally been studied by linescan confocal microscopy, but advances in TIRF microscopy together with improved electron-multiplied CCD (EMCCD) cameras now enable rapid (>500 frames s(-1)) imaging of subcellular Ca(2+) signals with high spatial resolution in two dimensions. This approach yields vastly more information (ca. 1 Gb min(-1)) than linescan imaging, rendering visual identification and analysis of local events imaged both laborious and subject to user bias. Here we describe a routine to rapidly automate identification and analysis of local Ca(2+) events. This features an intuitive graphical user-interfaces and runs under Matlab and the open-source Python software. The underlying algorithm features spatial and temporal noise filtering to reliably detect even small events in the presence of noisy and fluctuating baselines; localizes sites of Ca(2+) release with sub-pixel resolution; facilitates user review and editing of data; and outputs time-sequences of fluorescence ratio signals for identified event sites along with Excel-compatible tables listing amplitudes and kinetics of events.
Collapse
Affiliation(s)
- Kyle L Ellefsen
- Department of Neurobiology & Behavior, University of California, Irvine, CA 92697, USA
| | - Brett Settle
- Department of Neurobiology & Behavior, University of California, Irvine, CA 92697, USA
| | - Ian Parker
- Department of Neurobiology & Behavior, University of California, Irvine, CA 92697, USA; Department of Physiology & Biophysics, University of California, Irvine, CA 92697, USA
| | - Ian F Smith
- Department of Neurobiology & Behavior, University of California, Irvine, CA 92697, USA.
| |
Collapse
|
22
|
Thurley K, Tovey SC, Moenke G, Prince VL, Meena A, Thomas AP, Skupin A, Taylor CW, Falcke M. Reliable encoding of stimulus intensities within random sequences of intracellular Ca2+ spikes. Sci Signal 2014; 7:ra59. [PMID: 24962706 PMCID: PMC4092318 DOI: 10.1126/scisignal.2005237] [Citation(s) in RCA: 83] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Ca(2+) is a ubiquitous intracellular messenger that regulates diverse cellular activities. Extracellular stimuli often evoke sequences of intracellular Ca(2+) spikes, and spike frequency may encode stimulus intensity. However, the timing of spikes within a cell is random because each interspike interval has a large stochastic component. In human embryonic kidney (HEK) 293 cells and rat primary hepatocytes, we found that the average interspike interval also varied between individual cells. To evaluate how individual cells reliably encoded stimuli when Ca(2+) spikes exhibited such unpredictability, we combined Ca(2+) imaging of single cells with mathematical analyses of the Ca(2+) spikes evoked by receptors that stimulate formation of inositol 1,4,5-trisphosphate (IP3). This analysis revealed that signal-to-noise ratios were improved by slow recovery from feedback inhibition of Ca(2+) spiking operating at the whole-cell level and that they were robust against perturbations of the signaling pathway. Despite variability in the frequency of Ca(2+) spikes between cells, steps in stimulus intensity caused the stochastic period of the interspike interval to change by the same factor in all cells. These fold changes reliably encoded changes in stimulus intensity, and they resulted in an exponential dependence of average interspike interval on stimulation strength. We conclude that Ca(2+) spikes enable reliable signaling in a cell population despite randomness and cell-to-cell variability, because global feedback reduces noise, and changes in stimulus intensity are represented by fold changes in the stochastic period of the interspike interval.
Collapse
Affiliation(s)
- Kevin Thurley
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Robert Rössle Straße 10, Berlin 13125, Germany. Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK. Institute for Theoretical Biology, Charité-Universitätsmedizin Berlin, Invalidenstraße 43, Berlin 10115, Germany
| | - Stephen C Tovey
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK
| | - Gregor Moenke
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Robert Rössle Straße 10, Berlin 13125, Germany
| | - Victoria L Prince
- Department of Pharmacology and Physiology, New Jersey Medical School, Rutgers, The State University of New Jersey, 185 South Orange Avenue, Newark, NJ 07103, USA
| | - Abha Meena
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK
| | - Andrew P Thomas
- Department of Pharmacology and Physiology, New Jersey Medical School, Rutgers, The State University of New Jersey, 185 South Orange Avenue, Newark, NJ 07103, USA
| | - Alexander Skupin
- Luxembourg Centre for Systems Biomedicine, 7 Avenue des Hauts Fourneaux, Esch sur Alzette 4362, Luxembourg. National Center for Microscopy and Imaging Research, University of California, San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Colin W Taylor
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK.
| | - Martin Falcke
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine, Robert Rössle Straße 10, Berlin 13125, Germany. Department of Physics, Humboldt University Berlin, Newtonstraße 15, Berlin 12489, Germany.
| |
Collapse
|
23
|
Koumura T, Urakubo H, Ohashi K, Fujii M, Kuroda S. Stochasticity in Ca2+ increase in spines enables robust and sensitive information coding. PLoS One 2014; 9:e99040. [PMID: 24932482 PMCID: PMC4059641 DOI: 10.1371/journal.pone.0099040] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2014] [Accepted: 04/18/2014] [Indexed: 11/19/2022] Open
Abstract
A dendritic spine is a very small structure (∼0.1 µm3) of a neuron that processes input timing information. Why are spines so small? Here, we provide functional reasons; the size of spines is optimal for information coding. Spines code input timing information by the probability of Ca2+ increases, which makes robust and sensitive information coding possible. We created a stochastic simulation model of input timing-dependent Ca2+ increases in a cerebellar Purkinje cell's spine. Spines used probability coding of Ca2+ increases rather than amplitude coding for input timing detection via stochastic facilitation by utilizing the small number of molecules in a spine volume, where information per volume appeared optimal. Probability coding of Ca2+ increases in a spine volume was more robust against input fluctuation and more sensitive to input numbers than amplitude coding of Ca2+ increases in a cell volume. Thus, stochasticity is a strategy by which neurons robustly and sensitively code information.
Collapse
Affiliation(s)
- Takuya Koumura
- Undergraduate Department of Bioinformatics and Systems Biology, University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Hidetoshi Urakubo
- Department of Biophysics and Biochemistry, Graduate School of Science, University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Kaoru Ohashi
- Department of Biophysics and Biochemistry, Graduate School of Science, University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Masashi Fujii
- Department of Biophysics and Biochemistry, Graduate School of Science, University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Shinya Kuroda
- Undergraduate Department of Bioinformatics and Systems Biology, University of Tokyo, Bunkyo-ku, Tokyo, Japan
- Department of Biophysics and Biochemistry, Graduate School of Science, University of Tokyo, Bunkyo-ku, Tokyo, Japan
- * E-mail:
| |
Collapse
|
24
|
Thul R. Translating intracellular calcium signaling into models. Cold Spring Harb Protoc 2014; 2014:2014/5/pdb.top066266. [PMID: 24786496 DOI: 10.1101/pdb.top066266] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
The rich experimental data on intracellular calcium has put theoreticians in an ideal position to derive models of intracellular calcium signaling. Over the last 25 years, a large number of modeling frameworks have been suggested. Here, I will review some of the milestones of intracellular calcium modeling with a special emphasis on calcium-induced calcium release (CICR) through inositol-1,4,5-trisphosphate and ryanodine receptors. I will highlight key features of CICR and how they are represented in models as well as the challenges that theoreticians face when translating our current understanding of calcium signals into equations. The selected examples demonstrate that a successful model provides mechanistic insights into the molecular machinery of the Ca²⁺ signaling toolbox and determines the contribution of local Ca²⁺ release to global Ca²⁺ patterns, which at the moment cannot be resolved experimentally.
Collapse
Affiliation(s)
- Rüdiger Thul
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| |
Collapse
|
25
|
Bartlett PJ, Gaspers LD, Pierobon N, Thomas AP. Calcium-dependent regulation of glucose homeostasis in the liver. Cell Calcium 2014; 55:306-16. [PMID: 24630174 DOI: 10.1016/j.ceca.2014.02.007] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2014] [Revised: 02/07/2014] [Accepted: 02/08/2014] [Indexed: 02/09/2023]
Abstract
A major role of the liver is to integrate multiple signals to maintain normal blood glucose levels. The balance between glucose storage and mobilization is primarily regulated by the counteracting effects of insulin and glucagon. However, numerous signals converge in the liver to ensure energy demand matches the physiological status of the organism. Many circulating hormones regulate glycogenolysis, gluconeogenesis and mitochondrial metabolism by calcium-dependent signaling mechanisms that manifest as cytosolic Ca(2+) oscillations. Stimulus-strength is encoded in the Ca(2+) oscillation frequency, and also by the range of intercellular Ca(2+) wave propagation in the intact liver. In this article, we describe how Ca(2+) oscillations and waves can regulate glucose output and oxidative metabolism in the intact liver; how multiple stimuli are decoded though Ca(2+) signaling at the organ level, and the implications of Ca(2+) signal dysregulation in diseases such as metabolic syndrome and non-alcoholic fatty liver disease.
Collapse
Affiliation(s)
- Paula J Bartlett
- Department of Pharmacology and Physiology, New Jersey Medical School, Rutgers, The State University of New Jersey, Newark, NJ, USA.
| | - Lawrence D Gaspers
- Department of Pharmacology and Physiology, New Jersey Medical School, Rutgers, The State University of New Jersey, Newark, NJ, USA
| | - Nicola Pierobon
- Department of Pharmacology and Physiology, New Jersey Medical School, Rutgers, The State University of New Jersey, Newark, NJ, USA
| | - Andrew P Thomas
- Department of Pharmacology and Physiology, New Jersey Medical School, Rutgers, The State University of New Jersey, Newark, NJ, USA
| |
Collapse
|
26
|
Nivala M, Ko CY, Nivala M, Weiss JN, Qu Z. The emergence of subcellular pacemaker sites for calcium waves and oscillations. J Physiol 2013; 591:5305-20. [PMID: 24042497 DOI: 10.1113/jphysiol.2013.259960] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Calcium (Ca(2+)) waves generating oscillatory Ca(2+) signals are widely observed in biological cells. Experimental studies have shown that under certain conditions, initiation of Ca(2+) waves is random in space and time, while under other conditions, waves occur repetitively from preferred locations (pacemaker sites) from which they entrain the whole cell. In this study, we use computer simulations to investigate the self-organization of Ca(2+) sparks into pacemaker sites generating Ca(2+) oscillations. In both ventricular myocyte experiments and computer simulations of a heterogeneous Ca(2+) release unit (CRU) network model, we show that Ca(2+) waves occur randomly in space and time when the Ca(2+) level is low, but as the Ca(2+) level increases, waves occur repetitively from the same sites. Our analysis indicates that this transition to entrainment can be attributed to the fact that random Ca(2+) sparks self-organize into Ca(2+) oscillations differently at low and high Ca(2+) levels. At low Ca(2+), the whole cell Ca(2+) oscillation frequency of the coupled CRU system is much slower than that of an isolated single CRU. Compared to a single CRU, the distribution of interspike intervals (ISIs) of the coupled CRU network exhibits a greater variation, and its ISI distribution is asymmetric with respect to the peak, exhibiting a fat tail. At high Ca(2+), however, the coupled CRU network has a faster frequency and lesser ISI variation compared to an individual CRU. The ISI distribution of the coupled network no longer exhibits a fat tail and is well-approximated by a Gaussian distribution. This same Ca(2+) oscillation behaviour can also be achieved by varying the number of ryanodine receptors per CRU or the distance between CRUs. Using these results, we develop a theory for the entrainment of random oscillators which provides a unified explanation for the experimental observations underlying the emergence of pacemaker sites and Ca(2+) oscillations.
Collapse
Affiliation(s)
- Michael Nivala
- Z. Qu: Department of Medicine, David Geffen School of Medicine at UCLA, A2-237 CHS, 650 Charles E. Young Drive South, Los Angeles, CA 90095.
| | | | | | | | | |
Collapse
|
27
|
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
|
28
|
Flegg MB, Rüdiger S, Erban R. Diffusive spatio-temporal noise in a first-passage time model for intracellular calcium release. J Chem Phys 2013; 138:154103. [DOI: 10.1063/1.4796417] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
|
29
|
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
|
30
|
Dickinson GD, Swaminathan D, Parker I. The probability of triggering calcium puffs is linearly related to the number of inositol trisphosphate receptors in a cluster. Biophys J 2012; 102:1826-36. [PMID: 22768938 DOI: 10.1016/j.bpj.2012.03.029] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2011] [Revised: 03/12/2012] [Accepted: 03/14/2012] [Indexed: 10/28/2022] Open
Abstract
Puffs are local Ca(2+) signals that arise by Ca(2+) liberation from the endoplasmic reticulum through concerted opening of tightly clustered inositol trisphosphate receptor/channels (IP(3)R). They serve both local signaling functions and trigger global Ca(2+) waves. The numbers of functional IP(3)R within clusters differ appreciably between different puff sites, and we investigated how the probability of puff occurrence varies with cluster size. We imaged puffs in SH-SY5Y cells using total internal fluorescence microscopy, and estimated cluster sizes from the magnitude of the largest puff observed at each site relative to the signal from a single channel. We find that the initial triggering rate of puffs following photorelease of IP(3), and the average frequency of subsequent repetitive puffs, vary about linearly with cluster size. These data accord well with stochastic simulations in which opening of any individual IP(3)R channel within a cluster triggers a puff via Ca(2+)-induced Ca(2+) release. An important consequence is that the signaling power of a puff site (average amount of Ca(2+) released per puff × puff frequency) varies about the square of cluster size, implying that large clusters contribute disproportionately to cellular signaling and, because of their higher puff frequency, preferentially act as pacemakers to initiate Ca(2+) waves.
Collapse
Affiliation(s)
- George D Dickinson
- Department of Neurobiology and Behavior, University of California, Irvine, California, USA.
| | | | | |
Collapse
|
31
|
Lopez L, Piegari E, Sigaut L, Ponce Dawson S. Intracellular calcium signals display an avalanche-like behavior over multiple lengthscales. Front Physiol 2012; 3:350. [PMID: 22969730 PMCID: PMC3432517 DOI: 10.3389/fphys.2012.00350] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2012] [Accepted: 08/15/2012] [Indexed: 12/01/2022] Open
Abstract
Many natural phenomena display “self-organized criticality” (SOC), (Bak et al., 1987). This refers to spatially extended systems for which patterns of activity characterized by different lengthscales can occur with a probability density that follows a power law with pattern size. Differently from power laws at phase transitions, systems displaying SOC do not need the tuning of an external parameter. Here we analyze intracellular calcium (Ca2+) signals, a key component of the signaling toolkit of almost any cell type. Ca2+ signals can either be spatially restricted (local) or propagate throughout the cell (global). Different models have suggested that the transition from local to global signals is similar to that of directed percolation. Directed percolation has been associated, in turn, to the appearance of SOC. In this paper we discuss these issues within the framework of simple models of Ca2+ signal propagation. We also analyze the size distribution of local signals (“puffs”) observed in immature Xenopus Laevis oocytes. The puff amplitude distribution obtained from observed local signals is not Gaussian with a noticeable fraction of large size events. The experimental distribution of puff areas in the spatio-temporal record of the image has a long tail that is approximately log-normal. The distribution can also be fitted with a power law relationship albeit with a smaller goodness of fit. The power law behavior is encountered within a simple model that includes some coupling among individual signals for a wide range of parameter values. An analysis of the model shows that a global elevation of the Ca2+ concentration plays a major role in determining whether the puff size distribution is long-tailed or not. This suggests that Ca2+-clearing from the cytosol is key to determine whether IP3-mediated Ca2+ signals can display a SOC-like behavior or not.
Collapse
Affiliation(s)
- Lucía Lopez
- Departamento de Física FCEN-UBA and IFIBA, Ciudad Universitaria, Pabellón I Buenos Aires, Argentina
| | | | | | | |
Collapse
|
32
|
Nivala M, Ko CY, Nivala M, Weiss JN, Qu Z. Criticality in intracellular calcium signaling in cardiac myocytes. Biophys J 2012; 102:2433-42. [PMID: 22713558 DOI: 10.1016/j.bpj.2012.05.001] [Citation(s) in RCA: 70] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2012] [Revised: 04/30/2012] [Accepted: 05/01/2012] [Indexed: 11/18/2022] Open
Abstract
Calcium (Ca) is a ubiquitous second messenger that regulates many biological functions. The elementary events of local Ca signaling are Ca sparks, which occur randomly in time and space, and integrate to produce global signaling events such as intra- and intercellular Ca waves and whole-cell Ca oscillations. Despite extensive experimental characterization in many systems, the transition from local random to global synchronous events is still poorly understood. Here we show that criticality, a ubiquitous dynamical phenomenon in nature, is responsible for the transition from local to global Ca signaling. We demonstrate this first in a computational model of Ca signaling in a cardiac myocyte and then experimentally in mouse ventricular myocytes, complemented by a theoretical agent-based model to delineate the underlying dynamics. We show that the interaction between the Ca release units via Ca-induced Ca release causes self-organization of Ca spark clusters. When the coupling between Ca release units is weak, the cluster-size distribution is exponential. As the interactions become strong, the cluster-size distribution changes to a power-law distribution, which is characteristic of criticality in thermodynamic and complex nonlinear systems, and facilitates the formation and propagation of Ca waves and whole-cell Ca oscillations. Our findings illustrate how criticality is harnessed by a biological cell to regulate Ca signaling via self-organization of random subcellular events into cellular-scale oscillations, and provide a general theoretical framework for the transition from local Ca signaling to global Ca signaling in biological cells.
Collapse
Affiliation(s)
- Michael Nivala
- Cardiology Division, Department of Medicine, David Geffen School of Medicine, University of California, Los Angeles, California, USA
| | | | | | | | | |
Collapse
|
33
|
Termination of Ca²+ release for clustered IP₃R channels. PLoS Comput Biol 2012; 8:e1002485. [PMID: 22693433 PMCID: PMC3364945 DOI: 10.1371/journal.pcbi.1002485] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2011] [Accepted: 03/07/2012] [Indexed: 01/17/2023] Open
Abstract
In many cell types, release of calcium ions is controlled by inositol 1,4,5-trisphosphate () receptor channels. Elevations in concentration after intracellular release through receptors (R) can either propagate in the form of waves spreading through the entire cell or produce spatially localized puffs. The appearance of waves and puffs is thought to implicate random initial openings of one or a few channels and subsequent activation of neighboring channels because of an “autocatalytic” feedback. It is much less clear, however, what determines the further time course of release, particularly since the lifetime is very different for waves (several seconds) and puffs (around 100 ms). Here we study the lifetime of signals and their dependence on residual microdomains. Our general idea is that microdomains are dynamical and mediate the effect of other physiological processes. Specifically, we focus on the mechanism by which binding proteins (buffers) alter the lifetime of signals. We use stochastic simulations of channel gating coupled to a coarse-grained description for the concentration. To describe the concentration in a phenomenological way, we here introduce a differential equation, which reflects the buffer characteristics by a few effective parameters. This non-stationary model for microdomains gives deep insight into the dynamical differences between puffs and waves. It provides a novel explanation for the different lifetimes of puffs and waves and suggests that puffs are terminated by inhibition while unbinding is responsible for termination of waves. Thus our analysis hints at an additional role of and shows how cells can make use of the full complexity in R gating behavior to achieve different signals. Calcium signals are important for a host of cellular processes such as neurotransmitter release, cell contraction and gene expression. While the principles of activation and spreading of calcium signals have been largely understood, it is much less clear how their spatio-temporal appearance is shaped. This issue is of high relevance since the spatio-temporal signature is thought to carry the information content. In our paper we study the dynamical mechanisms that determine the time course of calcium release from receptor channels. We use a stochastic channel description combined with a recently developed model for the distribution of released calcium in a microdomain. The simulations uncover a complex control mechanism, which allows for the tuning of release from short frequent puffs to extended and less frequent wave-like release. Unexpectedly, the model predicts that for wave-like release the dissociation of from the receptors leads to termination of the calcium signal. This effect relies on a well-known gating property of R channels, which earlier has been regarded as superfluous in studies for groups of channels. Our results also provide a missing link to understand cellular response to calcium-binding proteins and present a novel mechanism for information processing by R channels.
Collapse
|
34
|
Skupin A, Thurley K. Calcium signaling: from single channels to pathways. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2012; 740:531-51. [PMID: 22453959 DOI: 10.1007/978-94-007-2888-2_24] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Ca(2+) is not only one of the most versatile and ubiquitous second messengers but also a well-established representative example of cell signaling. The identification of most key elements involved in Ca(2+) signaling enables a mechanistic and quantitative understanding of this particular pathway. Cellular behavior relies in general on the orchestration of molecular behavior leading to reliable cellular responses that allow for regulation and adaptation. Ca(2+) signaling uses a hierarchical organization to transform single molecule behavior into cell wide signals. We have recently shown experimentally that this organization carries single channel signatures onto the whole cell level and renders Ca(2+) oscillations stochastic. Here, we briefly review the co-evolution of experimental and theoretical studies in Ca(2+) -signaling and show how dynamic bottom-up modeling can be used to address -biological questions and illuminate biological principles of cell signaling.
Collapse
Affiliation(s)
- Alexander Skupin
- Luxembourg Centre of Systems Biomedicine, University Luxembourg, Luxembourg.
| | | |
Collapse
|
35
|
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
|