1
|
Szischik CL, Reves Szemere J, Balderrama R, Sánchez de la Vega C, Ventura AC. Transient frequency preference responses in cell signaling systems. NPJ Syst Biol Appl 2024; 10:86. [PMID: 39128915 PMCID: PMC11317535 DOI: 10.1038/s41540-024-00413-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 07/29/2024] [Indexed: 08/13/2024] Open
Abstract
Ligand-receptor systems, covalent modification cycles, and transcriptional networks are the fundamental components of cell signaling and gene expression systems. While their behavior in reaching a steady-state regime under step-like stimulation is well understood, their response under repetitive stimulation, particularly at early time stages is poorly characterized. Yet, early-stage responses to external inputs are arguably as informative as late-stage ones. In simple systems, a periodic stimulation elicits an initial transient response, followed by periodic behavior. Transient responses are relevant when the stimulation has a limited time span, or when the stimulated component's timescale is slow as compared to the timescales of the downstream processes, in which case the latter processes may be capturing only those transients. In this study, we analyze the frequency response of simple motifs at different time stages. We use dose-conserved pulsatile input signals and consider different metrics versus frequency curves. We show that in ligand-receptor systems, there is a frequency preference response in some specific metrics during the transient stages, which is not present in the periodic regime. We suggest this is a general system-level mechanism that cells may use to filter input signals that have consequences for higher order circuits. In addition, we evaluate how the described behavior in isolated motifs is reflected in similar types of responses in cascades and pathways of which they are a part. Our studies suggest that transient frequency preferences are important dynamic features of cell signaling and gene expression systems, which have been overlooked.
Collapse
Affiliation(s)
- Candela L Szischik
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física. Ciudad Universitaria, 1428, Buenos Aires, Argentina
- Instituto de Fisiología, Biología Molecular y Neurociencias (IFIBYNE UBA-CONICET), Consejo Nacional de Investigaciones Científicas y Técnicas of Argentina-Universidad de Buenos Aires, 1428, Buenos Aires, Argentina
| | - Juliana Reves Szemere
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física. Ciudad Universitaria, 1428, Buenos Aires, Argentina
- Instituto de Fisiología, Biología Molecular y Neurociencias (IFIBYNE UBA-CONICET), Consejo Nacional de Investigaciones Científicas y Técnicas of Argentina-Universidad de Buenos Aires, 1428, Buenos Aires, Argentina
- Universidad Pedagógica Nacional and Universidad Nacional de La Pampa, Facultad de Ciencias Exactas y Naturales, Departamento de Física, Santa Rosa, Argentina
| | - Rocío Balderrama
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Matemática. Ciudad Universitaria, Buenos Aires, Argentina
- Instituto de Investigaciones Matemáticas Luis A. Santaló (IMAS - CONICET), Consejo Nacional de Investigaciones Científicas y Técnicas of Argentina, Buenos Aires, Argentina
| | - Constanza Sánchez de la Vega
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Matemática. Ciudad Universitaria, Buenos Aires, Argentina
- Instituto de Cálculo, FCEyN, CONICET-UBA, Buenos Aires, Argentina
| | - Alejandra C Ventura
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física. Ciudad Universitaria, 1428, Buenos Aires, Argentina.
- Instituto de Fisiología, Biología Molecular y Neurociencias (IFIBYNE UBA-CONICET), Consejo Nacional de Investigaciones Científicas y Técnicas of Argentina-Universidad de Buenos Aires, 1428, Buenos Aires, Argentina.
| |
Collapse
|
2
|
Voorsluijs V, Avanzini F, Falasco G, Esposito M, Skupin A. Calcium oscillations optimize the energetic efficiency of mitochondrial metabolism. iScience 2024; 27:109078. [PMID: 38375217 PMCID: PMC10875125 DOI: 10.1016/j.isci.2024.109078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2023] [Revised: 12/26/2023] [Accepted: 01/26/2024] [Indexed: 02/21/2024] Open
Abstract
Energy transduction is central to living organisms, but the impact of enzyme regulation and signaling on its thermodynamic efficiency is generally overlooked. Here, we analyze the efficiency of ATP production by the tricarboxylic acid cycle and oxidative phosphorylation, which generate most of the chemical energy in eukaryotes. Calcium signaling regulates this pathway and can affect its energetic output, but the concrete energetic impact of this cross-talk remains elusive. Calcium enhances ATP production by activating key enzymes of the tricarboxylic acid cycle while calcium homeostasis is ATP-dependent. We propose a detailed kinetic model describing the calcium-mitochondria cross-talk and analyze it using nonequilibrium thermodynamics: after identifying the effective reactions driving mitochondrial metabolism out of equilibrium, we quantify the mitochondrial thermodynamic efficiency for different conditions. Calcium oscillations, triggered by extracellular stimulation or energy deficiency, boost the thermodynamic efficiency of mitochondrial metabolism, suggesting a compensatory role of calcium signaling in mitochondrial bioenergetics.
Collapse
Affiliation(s)
- Valérie Voorsluijs
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 avenue du Swing, 4367 Belvaux, Luxembourg
- Complex Systems and Statistical Mechanics, Department of Physics and Materials Science, University of Luxembourg, 162 A avenue de la Faïencerie, 1511 Luxembourg, Luxembourg
| | - Francesco Avanzini
- Complex Systems and Statistical Mechanics, Department of Physics and Materials Science, University of Luxembourg, 162 A avenue de la Faïencerie, 1511 Luxembourg, Luxembourg
- Department of Chemical Sciences, University of Padova, 1 Via F. Marzolo, 35131 Padova, Italy
| | - Gianmaria Falasco
- Complex Systems and Statistical Mechanics, Department of Physics and Materials Science, University of Luxembourg, 162 A avenue de la Faïencerie, 1511 Luxembourg, Luxembourg
- Department of Physics and Astronomy, University of Padova, 8 Via F. Marzolo, 35131 Padova, Italy
| | - Massimiliano Esposito
- Complex Systems and Statistical Mechanics, Department of Physics and Materials Science, University of Luxembourg, 162 A avenue de la Faïencerie, 1511 Luxembourg, Luxembourg
| | - Alexander Skupin
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 avenue du Swing, 4367 Belvaux, Luxembourg
- Department of Physics and Materials Science, University of Luxembourg, 162 A avenue de la Faïencerie, 1511 Luxembourg, Luxembourg
- Department of Neuroscience, University of California, San Diego, 9500 Gilman Drive, San Diego, CA 92093, USA
| |
Collapse
|
3
|
Ramlow L, Falcke M, Lindner B. An integrate-and-fire approach to Ca 2+ signaling. Part II: Cumulative refractoriness. Biophys J 2023; 122:4710-4729. [PMID: 37981761 PMCID: PMC10754692 DOI: 10.1016/j.bpj.2023.11.015] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Revised: 10/20/2023] [Accepted: 11/15/2023] [Indexed: 11/21/2023] Open
Abstract
Inositol 1,4,5-trisphosphate-induced Ca2+ signaling is a second messenger system used by almost all eukaryotic cells. The agonist concentration stimulating Ca2+ signals is encoded in the frequency of a Ca2+ concentration spike sequence. When a cell is stimulated, the interspike intervals (ISIs) often show a distinct transient during which they gradually increase, a system property we refer to as cumulative refractoriness. We extend a previously published stochastic model to include the Ca2+ concentration in the intracellular Ca2+ store as a slow adaptation variable. This model can reproduce both stationary and transient statistics of experimentally observed ISI sequences. We derive approximate expressions for the mean and coefficient of variation of the stationary ISIs. We also consider the response to the onset of a constant stimulus and estimate the length of the transient and the strength of the adaptation of the ISI. We show that the adaptation sets the coefficient of variation in agreement with current ideas derived from experiments. Moreover, we explain why, despite a pronounced transient behavior, ISI correlations can be weak, as often observed in experiments. Finally, we fit our model to reproduce the transient statistics of experimentally observed ISI sequences in stimulated HEK cells. The fitted model is able to qualitatively reproduce the relationship between the stationary interval correlations and the number of transient intervals, as well as the strength of the ISI adaptation. We also find positive correlations in the experimental sequence that cannot be explained by our model.
Collapse
Affiliation(s)
- Lukas Ramlow
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany; Department of Physics, Humboldt University Berlin, Berlin, Germany; Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Martin Falcke
- Department of Physics, Humboldt University Berlin, Berlin, Germany; Max Delbrück Center for Molecular Medicine, Berlin, Germany.
| | - Benjamin Lindner
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany; Department of Physics, Humboldt University Berlin, Berlin, Germany
| |
Collapse
|
4
|
Ho KYL, An K, Carr RL, Dvoskin AD, Ou AYJ, Vogl W, Tanentzapf G. Maintenance of hematopoietic stem cell niche homeostasis requires gap junction-mediated calcium signaling. Proc Natl Acad Sci U S A 2023; 120:e2303018120. [PMID: 37903259 PMCID: PMC10636368 DOI: 10.1073/pnas.2303018120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2023] [Accepted: 09/11/2023] [Indexed: 11/01/2023] Open
Abstract
Regulation of stem cells requires coordination of the cells that make up the stem cell niche. Here, we describe a mechanism that allows communication between niche cells to coordinate their activity and shape the signaling environment surrounding resident stem cells. Using the Drosophila hematopoietic organ, the lymph gland, we show that cells of the hematopoietic niche, the posterior signaling center (PSC), communicate using gap junctions (GJs) and form a signaling network. This network allows PSC cells to exchange Ca2+ signals repetitively which regulate the hematopoietic niche. Disruption of Ca2+ signaling in the PSC or the GJ-mediated network connecting niche cells causes dysregulation of the PSC and blood progenitor differentiation. Analysis of PSC-derived cell signaling shows that the Hedgehog pathway acts downstream of GJ-mediated Ca2+ signaling to modulate the niche microenvironment. These data show that GJ-mediated communication between hematopoietic niche cells maintains their homeostasis and consequently controls blood progenitor behavior.
Collapse
Affiliation(s)
- Kevin Y. L. Ho
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BCV6T 1Z3, Canada
| | - Kevin An
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BCV6T 1Z3, Canada
| | - Rosalyn L. Carr
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BCV6T 1Z3, Canada
- School of Biomedical Engineering, University of British Columbia, Vancouver, BCV6T 1Z3, Canada
- British Columbia Children’s Hospital Research Institute, British Columbia Children’s Hospital, Vancouver, BCV5Z 4H4, Canada
| | - Alexandra D. Dvoskin
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BCV6T 1Z3, Canada
| | - Annie Y. J. Ou
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BCV6T 1Z3, Canada
- School of Kinesiology, University of British Columbia, Vancouver, BCV6T 1Z1, Canada
| | - Wayne Vogl
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BCV6T 1Z3, Canada
| | - Guy Tanentzapf
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BCV6T 1Z3, Canada
| |
Collapse
|
5
|
Paknejad N, Sapuru V, Hite RK. Structural titration reveals Ca 2+-dependent conformational landscape of the IP 3 receptor. Nat Commun 2023; 14:6897. [PMID: 37898605 PMCID: PMC10613215 DOI: 10.1038/s41467-023-42707-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 10/18/2023] [Indexed: 10/30/2023] Open
Abstract
Inositol 1,4,5-trisphosphate receptors (IP3Rs) are endoplasmic reticulum Ca2+ channels whose biphasic dependence on cytosolic Ca2+ gives rise to Ca2+ oscillations that regulate fertilization, cell division and cell death. Despite the critical roles of IP3R-mediated Ca2+ responses, the structural underpinnings of the biphasic Ca2+ dependence that underlies Ca2+ oscillations are incompletely understood. Here, we collect cryo-EM images of an IP3R with Ca2+ concentrations spanning five orders of magnitude. Unbiased image analysis reveals that Ca2+ binding does not explicitly induce conformational changes but rather biases a complex conformational landscape consisting of resting, preactivated, activated, and inhibited states. Using particle counts as a proxy for relative conformational free energy, we demonstrate that Ca2+ binding at a high-affinity site allows IP3Rs to activate by escaping a low-energy resting state through an ensemble of preactivated states. At high Ca2+ concentrations, IP3Rs preferentially enter an inhibited state stabilized by a second, low-affinity Ca2+ binding site. Together, these studies provide a mechanistic basis for the biphasic Ca2+-dependence of IP3R channel activity.
Collapse
Affiliation(s)
- Navid Paknejad
- Structural Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
- Physiology, Biophysics, and Systems Biology (PBSB) Program, Weill Cornell Graduate School of Biomedical Sciences, 1300 York Avenue, New York, NY, 10065, USA
| | - Vinay Sapuru
- Structural Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
- Physiology, Biophysics, and Systems Biology (PBSB) Program, Weill Cornell Graduate School of Biomedical Sciences, 1300 York Avenue, New York, NY, 10065, USA
| | - Richard K Hite
- Structural Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA.
| |
Collapse
|
6
|
Friedhoff VN, Lindner B, Falcke M. Modeling IP 3-induced Ca 2+ signaling based on its interspike interval statistics. Biophys J 2023; 122:2818-2831. [PMID: 37312455 PMCID: PMC10398346 DOI: 10.1016/j.bpj.2023.06.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Revised: 04/24/2023] [Accepted: 06/06/2023] [Indexed: 06/15/2023] Open
Abstract
Inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ signaling is a second messenger system used by almost all eukaryotic cells. Recent research demonstrated randomness of Ca2+ signaling on all structural levels. We compile eight general properties of Ca2+ spiking common to all cell types investigated and suggest a theory of Ca2+ spiking starting from the random behavior of IP3 receptor channel clusters mediating the release of Ca2+ from the endoplasmic reticulum capturing all general properties and pathway-specific behavior. Spike generation begins after the absolute refractory period of the previous spike. According to its hierarchical spreading from initiating channel openings to cell level, we describe it as a first passage process from none to all clusters open while the cell recovers from the inhibition which terminated the previous spike. Our theory reproduces the exponential stimulation response relation of the average interspike interval Tav and its robustness properties, random spike timing with a linear moment relation between Tav and the interspike interval SD and its robustness properties, sensitive dependency of Tav on diffusion properties, and nonoscillatory local dynamics. We explain large cell variability of Tav observed in experiments by variability of channel cluster coupling by Ca2+-induced Ca2+ release, the number of clusters, and IP3 pathway component expression levels. We predict the relation between puff probability and agonist concentration and [IP3] and agonist concentration. Differences of spike behavior between cell types and stimulating agonists are explained by the different types of negative feedback terminating spikes. In summary, the hierarchical random character of spike generation explains all of the identified general properties.
Collapse
Affiliation(s)
- Victor Nicolai Friedhoff
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany; Department of Physics, Humboldt University, Berlin, Germany
| | - Benjamin Lindner
- Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany; Department of Physics, Humboldt University, Berlin, Germany
| | - Martin Falcke
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany; Department of Physics, Humboldt University, Berlin, Germany.
| |
Collapse
|
7
|
Smith HA, Thillaiappan NB, Rossi AM. IP 3 receptors: An "elementary" journey from structure to signals. Cell Calcium 2023; 113:102761. [PMID: 37271052 DOI: 10.1016/j.ceca.2023.102761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 05/15/2023] [Accepted: 05/18/2023] [Indexed: 06/06/2023]
Abstract
Inositol 1,4,5-trisphosphate receptors (IP3Rs) are large tetrameric channels which sit mostly in the membrane of the endoplasmic reticulum (ER) and mediate Ca2+ release from intracellular stores in response to extracellular stimuli in almost all cells. Dual regulation of IP3Rs by IP3 and Ca2+ itself, upstream "licensing", and the arrangement of IP3Rs into small clusters in the ER membrane, allow IP3Rs to generate spatially and temporally diverse Ca2+ signals. The characteristic biphasic regulation of IP3Rs by cytosolic Ca2+ concentration underpins regenerative Ca2+ signals by Ca2+-induced Ca2+-release, while also preventing uncontrolled explosive Ca2+ release. In this way, cells can harness a simple ion such as Ca2+ as a near-universal intracellular messenger to regulate diverse cellular functions, including those with conflicting outcomes such as cell survival and cell death. High-resolution structures of the IP3R bound to IP3 and Ca2+ in different combinations have together started to unravel the workings of this giant channel. Here we discuss, in the context of recently published structures, how the tight regulation of IP3Rs and their cellular geography lead to generation of "elementary" local Ca2+ signals known as Ca2+ "puffs", which form the fundamental bottleneck through which all IP3-mediated cytosolic Ca2+ signals must first pass.
Collapse
Affiliation(s)
- Holly A Smith
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, United Kingdom
| | | | - Ana M Rossi
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, United Kingdom.
| |
Collapse
|
8
|
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
|
9
|
Abstract
Inositol 1,4,5-trisphosphate (IP3) plays a key role in calcium signaling. After stimulation, it diffuses from the plasma membrane where it is produced to the endoplasmic reticulum where its receptors are localized. Based on in vitro measurements, IP3 was long thought to be a global messenger characterized by a diffusion coefficient of ~ 280 μm2s-1. However, in vivo observations revealed that this value does not match with the timing of localized Ca2+ increases induced by the confined release of a non-metabolizable IP3 analog. A theoretical analysis of these data concluded that in intact cells diffusion of IP3 is strongly hindered, leading to a 30-fold reduction of the diffusion coefficient. Here, we performed a new computational analysis of the same observations using a stochastic model of Ca2+ puffs. Our simulations concluded that the value of the effective IP3 diffusion coefficient is close to 100 μm2s-1. Such moderate reduction with respect to in vitro estimations quantitatively agrees with a buffering effect by non-fully bound inactive IP3 receptors. The model also reveals that IP3 spreading is not much affected by the endoplasmic reticulum, which represents an obstacle to the free displacement of molecules, but can be significantly increased in cells displaying elongated, 1-dimensional like geometries.
Collapse
|
10
|
Mehari FT, Miller M, Pick R, Bader A, Pekayvaz K, Napoli M, Uhl B, Reichel CA, Sperandio M, Walzog B, Schulz C, Massberg S, Stark K. Intravital calcium imaging in myeloid leukocytes identifies calcium frequency spectra as indicators of functional states. Sci Signal 2022; 15:eabe6909. [PMID: 35881691 DOI: 10.1126/scisignal.abe6909] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
The assessment of leukocyte activation in vivo is mainly based on surrogate parameters, such as cell shape changes and migration patterns. Consequently, additional parameters are required to dissect the complex spatiotemporal activation of leukocytes during inflammation. Here, we showed that intravital microscopy of myeloid leukocyte Ca2+ signals with Ca2+ reporter mouse strains combined with bioinformatic signal analysis provided a tool to assess their activation in vivo. We demonstrated by two-photon microscopy that tissue-resident macrophages reacted to sterile inflammation in the cremaster muscle with Ca2+ transients in a distinct spatiotemporal pattern. Moreover, through high-resolution, intravital spinning disk confocal microscopy, we identified the intracellular Ca2+ signaling patterns of neutrophils during the migration cascade in vivo. These patterns were modulated by the Ca2+ channel Orai1 and Gαi-coupled GPCRs, whose effects were evident through analysis of the range of frequencies of the Ca2+ signal (frequency spectra), which provided insights into the complex patterns of leukocyte Ca2+ oscillations. Together, these findings establish Ca2+ frequency spectra as an additional dimension to assess leukocyte activation and migration during inflammation in vivo.
Collapse
Affiliation(s)
- Fitsumbirhan T Mehari
- Medizinische Klinik und Poliklinik I, University Hospital, LMU Munich, 81377 Munich, Germany.,Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany
| | - Meike Miller
- Medizinische Klinik und Poliklinik I, University Hospital, LMU Munich, 81377 Munich, Germany.,Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany
| | - Robert Pick
- Institute for Cardiovascular Physiology and Pathophysiology, Walter Brendel Centre for Experimental Medicine, Biomedical Center (BMC), LMU Munich, 82152 Planegg, Germany
| | - Almke Bader
- Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany.,Institute for Cardiovascular Physiology and Pathophysiology, Walter Brendel Centre for Experimental Medicine, Biomedical Center (BMC), LMU Munich, 82152 Planegg, Germany
| | - Kami Pekayvaz
- Medizinische Klinik und Poliklinik I, University Hospital, LMU Munich, 81377 Munich, Germany.,Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany.,German Center for Cardiovascular Research (DZHK), partner site Munich Heart Alliance, 80802 Munich, Germany
| | - Matteo Napoli
- Institute for Cardiovascular Physiology and Pathophysiology, Walter Brendel Centre for Experimental Medicine, Biomedical Center (BMC), LMU Munich, 82152 Planegg, Germany
| | - Bernd Uhl
- Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany.,Department of Otorhinolaryngology, University Hospital, LMU Munich, 81377, Munich, Germany
| | - Christoph A Reichel
- Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany.,Department of Otorhinolaryngology, University Hospital, LMU Munich, 81377, Munich, Germany
| | - Markus Sperandio
- Institute for Cardiovascular Physiology and Pathophysiology, Walter Brendel Centre for Experimental Medicine, Biomedical Center (BMC), LMU Munich, 82152 Planegg, Germany
| | - Barbara Walzog
- Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany.,Institute for Cardiovascular Physiology and Pathophysiology, Walter Brendel Centre for Experimental Medicine, Biomedical Center (BMC), LMU Munich, 82152 Planegg, Germany
| | - Christian Schulz
- Medizinische Klinik und Poliklinik I, University Hospital, LMU Munich, 81377 Munich, Germany.,Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany.,German Center for Cardiovascular Research (DZHK), partner site Munich Heart Alliance, 80802 Munich, Germany
| | - Steffen Massberg
- Medizinische Klinik und Poliklinik I, University Hospital, LMU Munich, 81377 Munich, Germany.,Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany.,German Center for Cardiovascular Research (DZHK), partner site Munich Heart Alliance, 80802 Munich, Germany
| | - Konstantin Stark
- Medizinische Klinik und Poliklinik I, University Hospital, LMU Munich, 81377 Munich, Germany.,Walter Brendel Centre of Experimental Medicine, Faculty of Medicine, LMU Munich, 81377 Munich, Germany.,German Center for Cardiovascular Research (DZHK), partner site Munich Heart Alliance, 80802 Munich, Germany
| |
Collapse
|
11
|
Zhang Z, Hu G, Zhang Y, Qu Z. Kramers Rate Theory of Pacemaker Dynamics in Noisy Excitable Media. PHYSICAL REVIEW LETTERS 2022; 129:048101. [PMID: 35939013 PMCID: PMC11323706 DOI: 10.1103/physrevlett.129.048101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/06/2022] [Accepted: 06/29/2022] [Indexed: 06/01/2023]
Abstract
Rhythmic activities, which are usually driven by pacemakers, are common in biological systems. In noisy excitable media, pacemakers are self-organized firing clusters, but the underlying dynamics remains to be elucidated. Here we develop a Kramers rate theory of coupled cells to describe the firing properties of pacemakers and their dependence on coupling strength and system size and dimension. The theory captures accurately the simulation results of tissue models with stochastic Hodgkin-Huxley equations except when transitions from pacemakers to spiral waves occur under weak coupling.
Collapse
Affiliation(s)
- Zhaoyang Zhang
- Department of Physics, School of Physical Science and Technology, Ningbo University, Ningbo, Zhejiang 315211, China
| | - Gang Hu
- Department of Physics, Beijing Normal University, Beijing 100875, China
| | - Yuhao Zhang
- Department of Physics, School of Physical Science and Technology, Ningbo University, Ningbo, Zhejiang 315211, China
| | - Zhilin Qu
- Department of Medicine, University of California, Los Angeles, California 90095, USA
| |
Collapse
|
12
|
Bashkirtseva I, Ryashko L. How noise can generate calcium spike-type oscillations in deterministic equilibrium modes. Phys Rev E 2022; 105:054404. [PMID: 35706230 DOI: 10.1103/physreve.105.054404] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Accepted: 04/20/2022] [Indexed: 06/15/2023]
Abstract
Stochastic excitability of spiking oscillatory regimes in the calcium kinetics is studied on the basis of the Li-Rinzel conceptual model. The probabilistic mechanisms of the noise-induced generation of large-amplitude oscillations in parametric zones, where the original deterministic model has only stable equilibria, are investigated numerically and analytically. A parametric statistical description of interspike intervals is curried out and the phenomenon of coherence resonance is discussed. For the analytical study of the stochastic excitement, the confidence domain method using a stochastic sensitivity technique is applied. In this analysis, a key role of mutual arrangement of the confidence ellipses and separatrices detaching the sub- and supercritical regions is demonstrated. It is shown that in the Li-Rinzel model such separatrices are the stable manifolds of the saddle equilibria and the transient semiattractors.
Collapse
Affiliation(s)
- Irina Bashkirtseva
- Institute of Natural Sciences and Mathematics, Ural Federal University, Lenina 51, Ekaterinburg, Russia
| | - Lev Ryashko
- Institute of Natural Sciences and Mathematics, Ural Federal University, Lenina 51, Ekaterinburg, Russia
| |
Collapse
|
13
|
Schwabe D, Falcke M. On the relation between input and output distributions of scRNA-seq experiments. Bioinformatics 2022; 38:1336-1343. [PMID: 34908126 DOI: 10.1093/bioinformatics/btab841] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Revised: 12/01/2021] [Accepted: 12/12/2021] [Indexed: 01/05/2023] Open
Abstract
MOTIVATION Single-cell RNA sequencing determines RNA copy numbers per cell for a given gene. However, technical noise poses the question how observed distributions (output) are connected to their cellular distributions (input). RESULTS We model a single-cell RNA sequencing setup consisting of PCR amplification and sequencing, and derive probability distribution functions for the output distribution given an input distribution. We provide copy number distributions arising from single transcripts during PCR amplification with exact expressions for mean and variance. We prove that the coefficient of variation of the output of sequencing is always larger than that of the input distribution. Experimental data reveals the variance and mean of the input distribution to obey characteristic relations, which we specifically determine for a HeLa dataset. We can calculate as many moments of the input distribution as are known of the output distribution (up to all). This, in principle, completely determines the input from the output distribution. AVAILABILITY AND IMPLEMENTATION Source code freely available at https://github.com/danielschw188/InputOutputSCRNASeq. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Daniel Schwabe
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, 13125 Berlin, Germany
| | - Martin Falcke
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, 13125 Berlin, Germany.,Department of Physics, Humboldt University Berlin, 12489 Berlin, Germany
| |
Collapse
|
14
|
Ho KYL, Khadilkar RJ, Carr RL, Tanentzapf G. A gap-junction-mediated, calcium-signaling network controls blood progenitor fate decisions in hematopoiesis. Curr Biol 2021; 31:4697-4712.e6. [PMID: 34480855 DOI: 10.1016/j.cub.2021.08.027] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Revised: 05/28/2021] [Accepted: 08/06/2021] [Indexed: 11/24/2022]
Abstract
Stem cell homeostasis requires coordinated fate decisions among stem cells that are often widely distributed within a tissue at varying distances from their stem cell niche. This requires a mechanism to ensure robust fate decisions within a population of stem cells. Here, we show that, in the Drosophila hematopoietic organ, the lymph gland (LG), gap junctions form a network that coordinates fate decisions between blood progenitors. Using live imaging of calcium signaling in intact LGs, we find that blood progenitors are connected through a signaling network. Blocking gap junction function disrupts this network, alters the pattern of encoded calcium signals, and leads to loss of progenitors and precocious blood cell differentiation. Ectopic and uniform activation of the calcium-signaling mediator CaMKII restores progenitor homeostasis when gap junctions are disrupted. Overall, these data show that gap junctions equilibrate cell signals between blood progenitors to coordinate fate decisions and maintain hematopoietic homeostasis.
Collapse
Affiliation(s)
- Kevin Y L Ho
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BC V6T 1Z3, Canada
| | - Rohan J Khadilkar
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BC V6T 1Z3, Canada; Advanced Centre for Treatment, Research and Education in Cancer-Tata Memorial Centre (ACTREC-TMC), Kharghar, Navi Mumbai, Maharashtra 410210, India
| | - Rosalyn L Carr
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BC V6T 1Z3, Canada; School of Biomedical Engineering, University of British Columbia, Vancouver, BC V6T 1Z3, Canada
| | - Guy Tanentzapf
- Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BC V6T 1Z3, Canada.
| |
Collapse
|
15
|
Skinner DJ, Dunkel J. Estimating Entropy Production from Waiting Time Distributions. PHYSICAL REVIEW LETTERS 2021; 127:198101. [PMID: 34797138 DOI: 10.1103/physrevlett.127.198101] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/23/2021] [Accepted: 08/17/2021] [Indexed: 06/13/2023]
Abstract
Living systems operate far from thermal equilibrium by converting the chemical potential of ATP into mechanical work to achieve growth, replication, or locomotion. Given time series observations of intra-, inter-, or multicellular processes, a key challenge is to detect nonequilibrium behavior and quantify the rate of free energy consumption. Obtaining reliable bounds on energy consumption and entropy production directly from experimental data remains difficult in practice, as many degrees of freedom typically are hidden to the observer, so that the accessible coarse-grained dynamics may not obviously violate detailed balance. Here, we introduce a novel method for bounding the entropy production of physical and living systems which uses only the waiting time statistics of hidden Markov processes and, hence, can be directly applied to experimental data. By determining a universal limiting curve, we infer entropy production bounds from experimental data for gene regulatory networks, mammalian behavioral dynamics, and numerous other biological processes. Further considering the asymptotic limit of increasingly precise biological timers, we estimate the necessary entropic cost of heartbeat regulation in humans, dogs, and mice.
Collapse
Affiliation(s)
- Dominic J Skinner
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, USA
| | - Jörn Dunkel
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, USA
| |
Collapse
|
16
|
High rates of calcium-free diffusion in the cytosol of living cells. Biophys J 2021; 120:3960-3972. [PMID: 34454909 DOI: 10.1016/j.bpj.2021.08.019] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Revised: 06/07/2021] [Accepted: 08/11/2021] [Indexed: 11/20/2022] Open
Abstract
Calcium (Ca2+) is a universal second messenger that participates in the regulation of innumerous physiological processes. The way in which local elevations of the cytosolic Ca2+ concentration spread in space and time is key for the versatility of the signals. Ca2+ diffusion in the cytosol is hindered by its interaction with proteins that act as buffers. Depending on the concentrations and the kinetics of the interactions, there is a large range of values at which Ca2+ diffusion can proceed. Having reliable estimates of this range, particularly of its highest end, which corresponds to the ions free diffusion, is key to understand how the signals propagate. In this work, we present the first experimental results with which the Ca2+-free diffusion coefficient is directly quantified in the cytosol of living cells. By means of fluorescence correlation spectroscopy experiments performed in Xenopus laevis oocytes and in cells of Saccharomyces cerevisiae, we show that the ions can freely diffuse in the cytosol at a higher rate than previously thought.
Collapse
|
17
|
KRAP tethers IP 3 receptors to actin and licenses them to evoke cytosolic Ca 2+ signals. Nat Commun 2021; 12:4514. [PMID: 34301929 PMCID: PMC8302619 DOI: 10.1038/s41467-021-24739-9] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Accepted: 07/01/2021] [Indexed: 02/06/2023] Open
Abstract
Regulation of IP3 receptors (IP3Rs) by IP3 and Ca2+ allows regenerative Ca2+ signals, the smallest being Ca2+ puffs, which arise from coordinated openings of a few clustered IP3Rs. Cells express thousands of mostly mobile IP3Rs, yet Ca2+ puffs occur at a few immobile IP3R clusters. By imaging cells with endogenous IP3Rs tagged with EGFP, we show that KRas-induced actin-interacting protein (KRAP) tethers IP3Rs to actin beneath the plasma membrane. Loss of KRAP abolishes Ca2+ puffs and the global increases in cytosolic Ca2+ concentration evoked by more intense stimulation. Over-expressing KRAP immobilizes additional IP3R clusters and results in more Ca2+ puffs and larger global Ca2+ signals. Endogenous KRAP determines which IP3Rs will respond: it tethers IP3R clusters to actin alongside sites where store-operated Ca2+ entry occurs, licenses IP3Rs to evoke Ca2+ puffs and global cytosolic Ca2+ signals, implicates the actin cytoskeleton in IP3R regulation and may allow local activation of Ca2+ entry.
Collapse
|
18
|
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
|
19
|
Nienałtowski K, Rigby RE, Walczak J, Zakrzewska KE, Głów E, Rehwinkel J, Komorowski M. Fractional response analysis reveals logarithmic cytokine responses in cellular populations. Nat Commun 2021; 12:4175. [PMID: 34234126 PMCID: PMC8263596 DOI: 10.1038/s41467-021-24449-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2020] [Accepted: 06/17/2021] [Indexed: 01/10/2023] Open
Abstract
Although we can now measure single-cell signaling responses with multivariate, high-throughput techniques our ability to interpret such measurements is still limited. Even interpretation of dose–response based on single-cell data is not straightforward: signaling responses can differ significantly between cells, encompass multiple signaling effectors, and have dynamic character. Here, we use probabilistic modeling and information-theory to introduce fractional response analysis (FRA), which quantifies changes in fractions of cells with given response levels. FRA can be universally performed for heterogeneous, multivariate, and dynamic measurements and, as we demonstrate, quantifies otherwise hidden patterns in single-cell data. In particular, we show that fractional responses to type I interferon in human peripheral blood mononuclear cells are very similar across different cell types, despite significant differences in mean or median responses and degrees of cell-to-cell heterogeneity. Further, we demonstrate that fractional responses to cytokines scale linearly with the log of the cytokine dose, which uncovers that heterogeneous cellular populations are sensitive to fold-changes in the dose, as opposed to additive changes. Our ability to interpret single-cell multivariate signaling responses is still limited. Here the authors introduce fractional response analysis (FRA), involving fractional cell counting, capable of deconvoluting heterogeneous multivariate responses of cellular populations.
Collapse
Affiliation(s)
- Karol Nienałtowski
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Rachel E Rigby
- Medical Research Council Human Immunology Unit, Medical Research Council Weatherall Institute of Molecular Medicine, Radcliffe Department of Medicine, University of Oxford, Oxford, UK
| | - Jarosław Walczak
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Karolina E Zakrzewska
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Edyta Głów
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Jan Rehwinkel
- Medical Research Council Human Immunology Unit, Medical Research Council Weatherall Institute of Molecular Medicine, Radcliffe Department of Medicine, University of Oxford, Oxford, UK
| | - Michał Komorowski
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland.
| |
Collapse
|
20
|
Abstract
Living systems maintain or increase local order by working against the second law of thermodynamics. Thermodynamic consistency is restored as they consume free energy, thereby increasing the net entropy of their environment. Recently introduced estimators for the entropy production rate have provided major insights into the efficiency of important cellular processes. In experiments, however, many degrees of freedom typically remain hidden to the observer, and, in these cases, existing methods are not optimal. Here, by reformulating the problem within an optimization framework, we are able to infer improved bounds on the rate of entropy production from partial measurements of biological systems. Our approach yields provably optimal estimates given certain measurable transition statistics. In contrast to prevailing methods, the improved estimator reveals nonzero entropy production rates even when nonequilibrium processes appear time symmetric and therefore may pretend to obey detailed balance. We demonstrate the broad applicability of this framework by providing improved bounds on the energy consumption rates in a diverse range of biological systems including bacterial flagella motors, growing microtubules, and calcium oscillations within human embryonic kidney cells.
Collapse
Affiliation(s)
- Dominic J Skinner
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Jörn Dunkel
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 02139
| |
Collapse
|
21
|
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
|
22
|
Azpeitia E, Balanzario EP, Wagner A. Signaling pathways have an inherent need for noise to acquire information. BMC Bioinformatics 2020; 21:462. [PMID: 33066727 PMCID: PMC7568421 DOI: 10.1186/s12859-020-03778-x] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2020] [Accepted: 09/23/2020] [Indexed: 11/16/2022] Open
Abstract
BACKGROUND All living systems acquire information about their environment. At the cellular level, they do so through signaling pathways. Such pathways rely on reversible binding interactions between molecules that detect and transmit the presence of an extracellular cue or signal to the cell's interior. These interactions are inherently stochastic and thus noisy. On the one hand, noise can cause a signaling pathway to produce the same response for different stimuli, which reduces the amount of information a pathway acquires. On the other hand, in processes such as stochastic resonance, noise can improve the detection of weak stimuli and thus the acquisition of information. It is not clear whether the kinetic parameters that determine a pathway's operation cause noise to reduce or increase the acquisition of information. RESULTS We analyze how the kinetic properties of the reversible binding interactions used by signaling pathways affect the relationship between noise, the response to a signal, and information acquisition. Our results show that, under a wide range of biologically sensible parameter values, a noisy dynamic of reversible binding interactions is necessary to produce distinct responses to different stimuli. As a consequence, noise is indispensable for the acquisition of information in signaling pathways. CONCLUSIONS Our observations go beyond previous work by showing that noise plays a positive role in signaling pathways, demonstrating that noise is essential when such pathways acquire information.
Collapse
Affiliation(s)
- Eugenio Azpeitia
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Zurich, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
- Centro de Ciencias Matemáticas, Universidad Nacional Autónoma de México, Morelia, Mexico
| | - Eugenio P Balanzario
- Centro de Ciencias Matemáticas, Universidad Nacional Autónoma de México, Morelia, Mexico
| | - Andreas Wagner
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Zurich, Switzerland.
- Swiss Institute of Bioinformatics, Lausanne, Switzerland.
- The Santa Fe Institute, Santa Fe, NM, USA.
| |
Collapse
|
23
|
Dhyani V, Gare S, Gupta RK, Swain S, Venkatesh K, Giri L. GPCR mediated control of calcium dynamics: A systems perspective. Cell Signal 2020; 74:109717. [PMID: 32711109 PMCID: PMC7375278 DOI: 10.1016/j.cellsig.2020.109717] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Revised: 07/14/2020] [Accepted: 07/15/2020] [Indexed: 02/09/2023]
Abstract
G-protein coupled receptor (GPCR) mediated calcium (Ca2+)-signaling transduction remains crucial in designing drugs for various complex diseases including neurodegeneration, chronic heart failure as well as respiratory diseases. Although there are several reviews detailing various aspects of Ca2+-signaling such as the role of IP3 receptors and Ca2+-induced-Ca2+-release, none of them provide an integrated view of the mathematical descriptions of GPCR signal transduction and investigations on dose-response curves. This article is the first study in reviewing the network structures underlying GPCR signal transduction that control downstream [Cac2+]-oscillations. The central theme of this paper is to present the biochemical pathways, as well as molecular mechanisms underlying the GPCR-mediated Ca2+-dynamics in order to facilitate a better understanding of how agonist concentration is encoded in Ca2+-signals for Gαq, Gαs, and Gαi/o signaling pathways. Moreover, we present the GPCR targeting drugs that are relevant for treating cardiac, respiratory, and neuro-diseases. The current paper presents the ODE formulation for various models along with the detailed schematics of signaling networks. To provide a systems perspective, we present the network motifs that can provide readers an insight into the complex and intriguing science of agonist-mediated Ca2+-dynamics. One of the features of this review is to pinpoint the interplay between positive and negative feedback loops that are involved in controlling intracellular [Cac2+]-oscillations. Furthermore, we review several examples of dose-response curves obtained from [Cac2+]-spiking for various GPCR pathways. This paper is expected to be useful for pharmacologists and computational biologists for designing clinical applications of GPCR targeting drugs through modulation of Ca2+-dynamics.
Collapse
Affiliation(s)
- Vaibhav Dhyani
- Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Sangareddy, Telangana, India
| | - Suman Gare
- Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Sangareddy, Telangana, India
| | - Rishikesh Kumar Gupta
- Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Sangareddy, Telangana, India
| | - Sarpras Swain
- Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Sangareddy, Telangana, India
| | - K.V. Venkatesh
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Powai, Mumbai, Maharashtra, India
| | - Lopamudra Giri
- Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Sangareddy, Telangana, India.
| |
Collapse
|
24
|
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]
|
25
|
Maity A, Wollman R. Information transmission from NFkB signaling dynamics to gene expression. PLoS Comput Biol 2020; 16:e1008011. [PMID: 32797040 PMCID: PMC7478807 DOI: 10.1371/journal.pcbi.1008011] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2019] [Revised: 09/08/2020] [Accepted: 06/02/2020] [Indexed: 02/06/2023] Open
Abstract
The dynamic signal encoding paradigm suggests that information flows from the extracellular environment into specific signaling patterns (encoding) that are then read by downstream effectors to control cellular behavior. Previous work empirically quantified the information content of dynamic signaling patterns. However, whether this information can be faithfully transmitted to the gene expression level is unclear. Here we used NFkB signaling as a model to understand the accuracy of information transmission from signaling dynamics into gene expression. Using a detailed mathematical model, we simulated realistic NFkB signaling patterns with different degrees of variability. The NFkB patterns were used as an input to a simple gene expression model. Analysis of information transmission between ligand and NFkB and ligand and gene expression allows us to determine information loss in transmission between receptors to dynamic signaling patterns and between signaling dynamics to gene expression. Information loss could occur due to biochemical noise or due to a lack of specificity. We found that noise-free gene expression has very little information loss suggesting that gene expression can preserve specificity in NFkB patterns. As expected, the addition of noise to the gene expression model results in information loss. Interestingly, this effect can be mitigated by a specific choice of parameters that can substantially reduce information loss due to biochemical noise during gene expression. Overall our results show that the cellular capacity for information transmission from dynamic signaling patterns to gene expression can be high enough to preserve ligand specificity and thereby the accuracy of cellular response to environmental cues.
Collapse
Affiliation(s)
- Alok Maity
- Institute for Quantitative and Computational Biosciences, University of California, Los Angeles, California, United States of America
| | - Roy Wollman
- Institute for Quantitative and Computational Biosciences, University of California, Los Angeles, California, United States of America
- Departments of Integrative Biology and Physiology and Chemistry and Biochemistry, University of California UCLA, California, United States of America
- * E-mail:
| |
Collapse
|
26
|
Liu X, Parhi KK. Molecular and DNA Artificial Neural Networks via Fractional Coding. IEEE TRANSACTIONS ON BIOMEDICAL CIRCUITS AND SYSTEMS 2020; 14:490-503. [PMID: 32149654 DOI: 10.1109/tbcas.2020.2979485] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
This article considers implementation of artificial neural networks (ANNs) using molecular computing and DNA based on fractional coding. Prior work had addressed molecular two-layer ANNs with binary inputs and arbitrary weights. In prior work using fractional coding, a simple molecular perceptron that computes sigmoid of scaled weighted sum of the inputs was presented where the inputs and the weights lie between [-1,1]. Even for computing the perceptron, the prior approach suffers from two major limitations. First, it cannot compute the sigmoid of the weighted sum, but only the sigmoid of the scaled weighted sum. Second, many machine learning applications require the coefficients to be arbitrarily positive and negative numbers that are not bounded between [-1,1]; such numbers cannot be handled by the prior perceptron using fractional coding. This paper makes four contributions. First molecular perceptrons that can handle arbitrary weights and can compute sigmoid of the weighted sums are presented. Thus, these molecular perceptrons are ideal for regression applications and multi-layer ANNs. A new molecular divider is introduced and is used to compute sigmoid(ax) where . Second, based on fractional coding, a molecular artificial neural network (ANN) with one hidden layer is presented. Third, a trained ANN classifier with one hidden layer from seizure prediction application from electroencephalogram is mapped to molecular reactions and DNA and their performances are presented. Fourth, molecular activation functions for rectified linear unit (ReLU) and softmax are also presented.
Collapse
|
27
|
Bartlett PJ, Cloete I, Sneyd J, Thomas AP. IP 3-Dependent Ca 2+ Oscillations Switch into a Dual Oscillator Mechanism in the Presence of PLC-Linked Hormones. iScience 2020; 23:101062. [PMID: 32353764 PMCID: PMC7191650 DOI: 10.1016/j.isci.2020.101062] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2019] [Revised: 12/11/2019] [Accepted: 04/09/2020] [Indexed: 11/28/2022] Open
Abstract
Ca2+ oscillations that depend on inositol-1,4,5-trisphosphate (IP3) have been ascribed to biphasic Ca2+ regulation of the IP3 receptor (IP3R) or feedback mechanisms controlling IP3 levels in different cell types. IP3 uncaging in hepatocytes elicits Ca2+ transients that are often localized at the subcellular level and increase in magnitude with stimulus strength. However, this does not reproduce the broad baseline-separated global Ca2+ oscillations elicited by vasopressin. Addition of hormone to cells activated by IP3 uncaging initiates a qualitative transition from high-frequency spatially disorganized Ca2+ transients, to low-frequency, oscillatory Ca2+ waves that propagate throughout the cell. A mathematical model with dual coupled oscillators that integrates Ca2+-induced Ca2+ release at the IP3R and mutual feedback mechanisms of cross-coupling between Ca2+ and IP3 reproduces this behavior. Thus, multiple Ca2+ oscillation modes can coexist in the same cell, and hormonal stimulation can switch from the simpler to the more complex to yield robust signaling.
Collapse
Affiliation(s)
- Paula J Bartlett
- Department of Pharmacology, Physiology and Neuroscience, New Jersey Medical School Rutgers, The State University of New Jersey, Newark, NJ 07103, USA
| | - Ielyaas Cloete
- Department of Mathematics, The University of Auckland, Auckland, New Zealand
| | - James Sneyd
- Department of Mathematics, The University of Auckland, Auckland, 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, USA.
| |
Collapse
|
28
|
Papp B, Launay S, Gélébart P, Arbabian A, Enyedi A, Brouland JP, Carosella ED, Adle-Biassette H. Endoplasmic Reticulum Calcium Pumps and Tumor Cell Differentiation. Int J Mol Sci 2020; 21:ijms21093351. [PMID: 32397400 PMCID: PMC7247589 DOI: 10.3390/ijms21093351] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2020] [Revised: 05/01/2020] [Accepted: 05/02/2020] [Indexed: 12/21/2022] Open
Abstract
Endoplasmic reticulum (ER) calcium homeostasis plays an essential role in cellular calcium signaling, intra-ER protein chaperoning and maturation, as well as in the interaction of the ER with other organelles. Calcium is accumulated in the ER by sarco/endoplasmic reticulum calcium ATPases (SERCA enzymes) that generate by active, ATP-dependent transport, a several thousand-fold calcium ion concentration gradient between the cytosol (low nanomolar) and the ER lumen (high micromolar). SERCA enzymes are coded by three genes that by alternative splicing give rise to several isoforms, which can display isoform-specific calcium transport characteristics. SERCA expression levels and isoenzyme composition vary according to cell type, and this constitutes a mechanism whereby ER calcium homeostasis is adapted to the signaling and metabolic needs of the cell, depending on its phenotype, its state of activation and differentiation. As reviewed here, in several normal epithelial cell types including bronchial, mammary, gastric, colonic and choroid plexus epithelium, as well as in mature cells of hematopoietic origin such as pumps are simultaneously expressed, whereas in corresponding tumors and leukemias SERCA3 expression is selectively down-regulated. SERCA3 expression is restored during the pharmacologically induced differentiation of various cancer and leukemia cell types. SERCA3 is a useful marker for the study of cell differentiation, and the loss of SERCA3 expression constitutes a previously unrecognized example of the remodeling of calcium homeostasis in tumors.
Collapse
Affiliation(s)
- Bela Papp
- Institut National de la Santé et de la Recherche Médicale, UMR U976, Institut Saint-Louis, 75010 Paris, France
- Institut de Recherche Saint-Louis, Hôpital Saint-Louis, Université de Paris, 75010 Paris, France
- CEA, DRF-Institut Francois Jacob, Department of Hemato-Immunology Research, Hôpital Saint-Louis, 75010 Paris, France;
- Correspondence: or
| | - Sophie Launay
- EA481, UFR Santé, Université de Bourgogne Franche-Comté, 25000 Besançon, France;
| | - Pascal Gélébart
- Department of Clinical Science-Hematology Section, Haukeland University Hospital, University of Bergen, 5021 Bergen, Norway;
| | - Atousa Arbabian
- Laboratoire d’Innovation Vaccins, Institut Pasteur de Paris, 75015 Paris, France;
| | - Agnes Enyedi
- Second Department of Pathology, Semmelweis University, 1091 Budapest, Hungary;
| | - Jean-Philippe Brouland
- Institut Universitaire de Pathologie, Centre Hospitalier Universitaire Vaudois, 1011 Lausanne, Switzerland;
| | - Edgardo D. Carosella
- CEA, DRF-Institut Francois Jacob, Department of Hemato-Immunology Research, Hôpital Saint-Louis, 75010 Paris, France;
| | - Homa Adle-Biassette
- AP-HP, Service d’Anatomie et Cytologie Pathologiques, Hôpital Lariboisière, 75010 Paris, France;
- Université de Paris, NeuroDiderot, Inserm UMR 1141, 75019 Paris, France
| |
Collapse
|
29
|
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
|
30
|
Moein M, Grzyb K, Gonçalves Martins T, Komoto S, Peri F, Crawford AD, Fouquier d'Herouel A, Skupin A. CaSiAn: a Calcium Signaling Analyzer tool. Bioinformatics 2019; 34:3052-3054. [PMID: 29668830 PMCID: PMC6129310 DOI: 10.1093/bioinformatics/bty281] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2017] [Accepted: 04/11/2018] [Indexed: 12/04/2022] Open
Abstract
Summary Ca2+ is a central second messenger in eukaryotic cells that regulates many cellular processes. Recently, we have indicated that typical Ca2+ signals are not purely oscillatory as widely assumed, but exhibit stochastic spiking with cell type and pathway specific characteristics. Here, we present the Calcium Signaling Analyzer (CaSiAn), an open source software tool that allows for quantifying these signal characteristics including individual spike properties and time course statistics in a semi-automated manner. CaSiAn provides an intuitive graphical user interface allowing experimentalists to easily process a large amount of Ca2+ signals, interactively tune peak detection, revise statistical measures and access the quantified signal properties as excel or text files. Availability and implementation CaSiAn is implemented in Java and available on Github (https://github.com/mmahsa/CaSiAn) as well as on the project page (http://r3lab.uni.lu/web/casa). Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Mahsa Moein
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux L-4367, Luxembourg
| | - Kamil Grzyb
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux L-4367, Luxembourg
| | | | - Shinya Komoto
- Okinawa Institute of Science and Technology Graduate University, Onna, Japan
| | - Francesca Peri
- Developmental Biology Unit, European Molecular Biology Laboratory (EMBL), Heidelberg, Germany
| | - Alexander D Crawford
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux L-4367, Luxembourg
| | | | - Alexander Skupin
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux L-4367, Luxembourg.,University California San Diego, La Jolla, CA, USA
| |
Collapse
|
31
|
McCarron JG, Wilson C, Heathcote HR, Zhang X, Buckley C, Lee MD. Heterogeneity and emergent behaviour in the vascular endothelium. Curr Opin Pharmacol 2019; 45:23-32. [PMID: 31005824 PMCID: PMC6700393 DOI: 10.1016/j.coph.2019.03.008] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2018] [Accepted: 03/18/2019] [Indexed: 12/16/2022]
Abstract
The endothelium is the single layer of cells lining all blood vessels, and it is a remarkable cardiovascular control centre. Each endothelial cell has only a small number (on average six) of interconnected neighbours. Yet this arrangement produces a large repertoire of behaviours, capable of controlling numerous cardiovascular functions in a flexible and dynamic way. The endothelium regulates the delivery of nutrients and removal of waste by regulating blood flow and vascular permeability. The endothelium regulates blood clotting, responses to infection and inflammation, the formation of new blood vessels, and remodelling of the blood vessel wall. To carry out these roles, the endothelium autonomously interprets a complex environment crammed with signals from hormones, neurotransmitters, pericytes, smooth muscle cells, various blood cells, viral or bacterial infection and proinflammatory cytokines. It is generally assumed that the endothelium responds to these instructions with coordinated responses in a homogeneous population of endothelial cells. Here, we highlight evidence that shows that neighbouring endothelial cells are highly heterogeneous and display different sensitivities to various activators. Cells with various sensitivities process different extracellular signals into distinct streams of information in parallel, like a vast switchboard. Communication occurs among cells and new ‘emergent’ signals are generated that are non-linear composites of the inputs. Emergent signals cannot be predicted or deduced from the properties of individual cells. Heterogeneity and emergent behaviour bestow capabilities on the endothelial collective that far exceed those of individual cells. The implications of heterogeneity and emergent behaviour for understanding vascular disease and drug discovery are discussed.
Collapse
Affiliation(s)
- John G McCarron
- Strathclyde Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral Street, Glasgow G4 0RE, UK.
| | - Calum Wilson
- Strathclyde Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral Street, Glasgow G4 0RE, UK
| | - Helen R Heathcote
- Strathclyde Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral Street, Glasgow G4 0RE, UK
| | - Xun Zhang
- Strathclyde Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral Street, Glasgow G4 0RE, UK
| | - Charlotte Buckley
- Strathclyde Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral Street, Glasgow G4 0RE, UK
| | - Matthew D Lee
- Strathclyde Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral Street, Glasgow G4 0RE, UK
| |
Collapse
|
32
|
Prole DL, Taylor CW. Structure and Function of IP 3 Receptors. Cold Spring Harb Perspect Biol 2019; 11:cshperspect.a035063. [PMID: 30745293 DOI: 10.1101/cshperspect.a035063] [Citation(s) in RCA: 100] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Inositol 1,4,5-trisphosphate receptors (IP3Rs), by releasing Ca2+ from the endoplasmic reticulum (ER) of animal cells, allow Ca2+ to be redistributed from the ER to the cytosol or other organelles, and they initiate store-operated Ca2+ entry (SOCE). For all three IP3R subtypes, binding of IP3 primes them to bind Ca2+, which then triggers channel opening. We are now close to understanding the structural basis of IP3R activation. Ca2+-induced Ca2+ release regulated by IP3 allows IP3Rs to regeneratively propagate Ca2+ signals. The smallest of these regenerative events is a Ca2+ puff, which arises from the nearly simultaneous opening of a small cluster of IP3Rs. Ca2+ puffs are the basic building blocks for all IP3-evoked Ca2+ signals, but only some IP3 clusters, namely those parked alongside the ER-plasma membrane junctions where SOCE occurs, are licensed to respond. The location of these licensed IP3Rs may allow them to selectively regulate SOCE.
Collapse
Affiliation(s)
- David L Prole
- Department of Pharmacology, University of Cambridge, Cambridge CB2 1PD, United Kingdom
| | - Colin W Taylor
- Department of Pharmacology, University of Cambridge, Cambridge CB2 1PD, United Kingdom
| |
Collapse
|
33
|
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
|
34
|
Sumit M, Jovic A, Neubig RR, Takayama S, Linderman JJ. A Two-Pulse Cellular Stimulation Test Elucidates Variability and Mechanisms in Signaling Pathways. Biophys J 2019; 116:962-973. [PMID: 30782397 DOI: 10.1016/j.bpj.2019.01.022] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2018] [Revised: 01/14/2019] [Accepted: 01/18/2019] [Indexed: 12/14/2022] Open
Abstract
Mammalian cells respond in a variable manner when provided with physiological pulses of ligand, such as low concentrations of acetylcholine present for just tens of seconds or TNFα for just tens of minutes. For a two-pulse stimulation, some cells respond to both pulses, some do not respond, and yet others respond to only one or the other pulse. Are these different response patterns the result of the small number of ligands being able to only stochastically activate the pathway at random times or an output pattern from a deterministic algorithm responding differently to different stimulation intervals? If the response is deterministic in nature, what parameters determine whether a response is generated or skipped? To answer these questions, we developed a two-pulse test that utilizes different rest periods between stimulation pulses. This "rest-period test" revealed that cells skip responses predictably as the rest period is shortened. By combining these experimental results with a mathematical model of the pathway, we further obtained mechanistic insight into potential sources of response variability. Our analysis indicates that in both intracellular calcium and NFκB signaling, response variability is consistent with extrinsic noise (cell-to-cell variability in protein levels), a short-term memory of stimulation, and high Hill coefficient processes. Furthermore, these results support recent works that have emphasized the role of deterministic processes for explaining apparently stochastic cellular response variability and indicate that even weak stimulations likely guide mammalian cells to appropriate fates rather than leaving outcomes to chance. We envision that the rest-period test can be applied to other signaling pathways to extract mechanistic insight.
Collapse
Affiliation(s)
- Madhuresh Sumit
- Biophysics Graduate Program, University of Michigan, Ann Arbor, Michigan
| | - Andreja Jovic
- Program in Molecular Pharmacology, Memorial Sloan-Kettering Cancer Center, New York, New York
| | - Richard R Neubig
- Department of Pharmacology and Toxicology, Michigan State University, East Lansing, Michigan
| | - Shuichi Takayama
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, Michigan; Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory School of Medicine, Atlanta, Georgia.
| | - Jennifer J Linderman
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, Michigan; Department of Chemical Engineering, University of Michigan, Ann Arbor, Michigan.
| |
Collapse
|
35
|
Wacquier B, Voorsluijs V, Combettes L, Dupont G. Coding and decoding of oscillatory Ca 2+ signals. Semin Cell Dev Biol 2019; 94:11-19. [PMID: 30659886 DOI: 10.1016/j.semcdb.2019.01.008] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2018] [Revised: 01/11/2019] [Accepted: 01/14/2019] [Indexed: 01/08/2023]
Abstract
About 30 years after their first observation, Ca2+ oscillations are now recognised as a universal mechanism of signal transduction. These oscillations are driven by periodic cycles of release and uptake of Ca2+ between the cytoplasm and the endoplasmic reticulum. Their frequency often increases with the level of stimulation, which can be decoded by some molecules. However, it is becoming increasingly evident that the widespread core oscillatory mechanism is modulated in many ways, depending on the cell type and on the physiological conditions. Interplay with inositol 1,4,5-trisphosphate metabolism and with other Ca2+ stores as the extracellular medium or mitochondria can much affect the properties of these oscillations. In many cases, these finely tuned characteristics of Ca2+ oscillations impact the physiological response that is triggered by the signal. Moreover, oscillations are intrinsically irregular. This randomness can also be exploited by the cell. In this review, we discuss evidences of these additional manifestations of the versatility of Ca2+ signalling.
Collapse
Affiliation(s)
- Benjamin Wacquier
- Unit of Theoretical Chronobiology, Faculté des Sciences, Université Libre de Bruxelles (ULB), Brussels, Belgium
| | - Valérie Voorsluijs
- Nonlinear Physical Chemistry Unit & Center for Nonlinear Phenomena and Complex Systems (CENOLI), Université Libre de Bruxelles (ULB), Brussels, Belgium
| | | | - Geneviève Dupont
- Unit of Theoretical Chronobiology, Faculté des Sciences, Université Libre de Bruxelles (ULB), Brussels, Belgium.
| |
Collapse
|
36
|
Rossi AM, Taylor CW. IP3 receptors – lessons from analyses ex cellula. J Cell Sci 2018; 132:132/4/jcs222463. [DOI: 10.1242/jcs.222463] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
ABSTRACT
Inositol 1,4,5-trisphosphate receptors (IP3Rs) are widely expressed intracellular channels that release Ca2+ from the endoplasmic reticulum (ER). We review how studies of IP3Rs removed from their intracellular environment (‘ex cellula’), alongside similar analyses of ryanodine receptors, have contributed to understanding IP3R behaviour. Analyses of permeabilized cells have demonstrated that the ER is the major intracellular Ca2+ store, and that IP3 stimulates Ca2+ release from this store. Radioligand binding confirmed that the 4,5-phosphates of IP3 are essential for activating IP3Rs, and facilitated IP3R purification and cloning, which paved the way for structural analyses. Reconstitution of IP3Rs into lipid bilayers and patch-clamp recording from the nuclear envelope have established that IP3Rs have a large conductance and select weakly between Ca2+ and other cations. Structural analyses are now revealing how IP3 binding to the N-terminus of the tetrameric IP3R opens the pore ∼7 nm away from the IP3-binding core (IBC). Communication between the IBC and pore passes through a nexus of interleaved domains contributed by structures associated with the pore and cytosolic domains, which together contribute to a Ca2+-binding site. These structural analyses provide evidence to support the suggestion that IP3 gates IP3Rs by first stimulating Ca2+ binding, which leads to pore opening and Ca2+ release.
Collapse
Affiliation(s)
- Ana M. Rossi
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK
| | - Colin W. Taylor
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK
| |
Collapse
|
37
|
Thillaiappan NB, Chakraborty P, Hasan G, Taylor CW. IP 3 receptors and Ca 2+ entry. BIOCHIMICA ET BIOPHYSICA ACTA-MOLECULAR CELL RESEARCH 2018; 1866:1092-1100. [PMID: 30448464 DOI: 10.1016/j.bbamcr.2018.11.007] [Citation(s) in RCA: 48] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/19/2018] [Revised: 11/07/2018] [Accepted: 11/08/2018] [Indexed: 12/23/2022]
Abstract
Inositol 1,4,5-trisphosphate receptors (IP3R) are the most widely expressed intracellular Ca2+ release channels. Their activation by IP3 and Ca2+ allows Ca2+ to pass rapidly from the ER lumen to the cytosol. The resulting increase in cytosolic [Ca2+] may directly regulate cytosolic effectors or fuel Ca2+ uptake by other organelles, while the decrease in ER luminal [Ca2+] stimulates store-operated Ca2+ entry (SOCE). We are close to understanding the structural basis of both IP3R activation, and the interactions between the ER Ca2+-sensor, STIM, and the plasma membrane Ca2+ channel, Orai, that lead to SOCE. IP3Rs are the usual means through which extracellular stimuli, through ER Ca2+ release, stimulate SOCE. Here, we review evidence that the IP3Rs most likely to respond to IP3 are optimally placed to allow regulation of SOCE. We also consider evidence that IP3Rs may regulate SOCE downstream of their ability to deplete ER Ca2+ stores. Finally, we review evidence that IP3Rs in the plasma membrane can also directly mediate Ca2+ entry in some cells.
Collapse
Affiliation(s)
| | - Pragnya Chakraborty
- Department of Pharmacology, Tennis Court Road, Cambridge CB2 1PD, United Kingdom; National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bellary Road, Bangalore 560065, India
| | - Gaiti Hasan
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bellary Road, Bangalore 560065, India
| | - Colin W Taylor
- Department of Pharmacology, Tennis Court Road, Cambridge CB2 1PD, United Kingdom.
| |
Collapse
|
38
|
Abstract
One challenge in biology is to make sense of the complexity of biological networks. A good system to approach this is signaling pathways, whose well-characterized molecular details allow us to relate the internal processes of each pathway to their input-output behavior. In this study, we analyzed mathematical models of three metazoan signaling pathways: the canonical Wnt, MAPK/ERK, and Tgfβ pathways. We find an unexpected convergence: the three pathways behave in some physiological contexts as linear signal transmitters. Testing the results experimentally, we present direct measurements of linear input-output behavior in the Wnt and ERK pathways. Analytics from each model further reveal that linearity arises through different means in each pathway, which we tested experimentally in the Wnt and ERK pathways. Linearity is a desired property in engineering where it facilitates fidelity and superposition in signal transmission. Our findings illustrate how cells tune different complex networks to converge on the same behavior.
Collapse
Affiliation(s)
- Harry Nunns
- Division of Biology and Biological EngineeringCalifornia Institute of TechnologyPasadenaUnited States
| | - Lea Goentoro
- Division of Biology and Biological EngineeringCalifornia Institute of TechnologyPasadenaUnited States
| |
Collapse
|
39
|
Hannanta-anan P, Chow BY. Optogenetic Inhibition of Gα q Protein Signaling Reduces Calcium Oscillation Stochasticity. ACS Synth Biol 2018; 7:1488-1495. [PMID: 29792810 PMCID: PMC6311707 DOI: 10.1021/acssynbio.8b00065] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
As fast terminators of G-protein coupled receptor (GPCR) signaling, regulators of G-protein signaling (RGS) serve critical roles in fine-tuning second messenger levels and, consequently, cellular responses to external stimuli. Here, we report the creation of an optogenetic RGS2 (opto-RGS2) that suppresses agonist-evoked calcium oscillations by the inactivation of Gαq protein. In this system, cryptochrome-mediated heterodimerization of the catalytic RGS2-box with its N-terminal amphipathic helix reconstitutes a functional membrane-localized complex that can dynamically suppress store-operated release of calcium. Engineered opto-RGS2 cell lines were used to establish the role of RGS2 as a key inhibitory feedback regulator of the stochasticity of the Gαq-mediated calcium spike timing. RGS2 reduced the stochasticity of carbachol-stimulated calcium oscillations, and the feedback inhibition was coupled to the global calcium elevation by calmodulin/RGS2 interactions. The identification of a critical negative feedback circuit exemplifies the utility of optogenetic approaches for interrogating RGS/GPCR biology and calcium encoding principles through temporally precise molecular gain-of-function.
Collapse
Affiliation(s)
| | - Brian Y. Chow
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA 19104
| |
Collapse
|
40
|
Salehi SA, Liu X, Riedel MD, Parhi KK. Computing Mathematical Functions using DNA via Fractional Coding. Sci Rep 2018; 8:8312. [PMID: 29844537 PMCID: PMC5974329 DOI: 10.1038/s41598-018-26709-6] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2017] [Accepted: 05/18/2018] [Indexed: 12/22/2022] Open
Abstract
This paper discusses the implementation of mathematical functions such as exponentials, trigonometric functions, the sigmoid function and the perceptron function with molecular reactions in general, and DNA strand displacement reactions in particular. The molecular constructs for these functions are predicated on a novel representation for input and output values: a fractional encoding, in which values are represented by the relative concentrations of two molecular types, denoted as type-1 and type-0. This representation is inspired by a technique from digital electronic design, termed stochastic logic, in which values are represented by the probability of 1's in a stream of randomly generated 0's and 1's. Research in the electronic realm has shown that a variety of complex functions can be computed with remarkably simple circuitry with this stochastic approach. This paper demonstrates how stochastic electronic designs can be translated to molecular circuits. It presents molecular implementations of mathematical functions that are considerably more complex than any shown to date. All designs are validated using mass-action simulations of the chemical kinetics of DNA strand displacement reactions.
Collapse
Affiliation(s)
- Sayed Ahmad Salehi
- Department of Electrical and Computer Engineering, University of Minnesota, 200 Union St. S.E., Minneapolis, MN, 55455, USA
| | - Xingyi Liu
- Department of Electrical and Computer Engineering, University of Minnesota, 200 Union St. S.E., Minneapolis, MN, 55455, USA
| | - Marc D Riedel
- Department of Electrical and Computer Engineering, University of Minnesota, 200 Union St. S.E., Minneapolis, MN, 55455, USA
| | - Keshab K Parhi
- Department of Electrical and Computer Engineering, University of Minnesota, 200 Union St. S.E., Minneapolis, MN, 55455, USA.
| |
Collapse
|
41
|
Falcke M, Friedhoff VN. The stretch to stray on time: Resonant length of random walks in a transient. CHAOS (WOODBURY, N.Y.) 2018; 28:053117. [PMID: 29857685 DOI: 10.1063/1.5023164] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
First-passage times in random walks have a vast number of diverse applications in physics, chemistry, biology, and finance. In general, environmental conditions for a stochastic process are not constant on the time scale of the average first-passage time or control might be applied to reduce noise. We investigate moments of the first-passage time distribution under an exponential transient describing relaxation of environmental conditions. We solve the Laplace-transformed (generalized) master equation analytically using a novel method that is applicable to general state schemes. The first-passage time from one end to the other of a linear chain of states is our application for the solutions. The dependence of its average on the relaxation rate obeys a power law for slow transients. The exponent ν depends on the chain length N like ν=-N/(N+1) to leading order. Slow transients substantially reduce the noise of first-passage times expressed as the coefficient of variation (CV), even if the average first-passage time is much longer than the transient. The CV has a pronounced minimum for some lengths, which we call resonant lengths. These results also suggest a simple and efficient noise control strategy and are closely related to the timing of repetitive excitations, coherence resonance, and information transmission by noisy excitable systems. A resonant number of steps from the inhibited state to the excitation threshold and slow recovery from negative feedback provide optimal timing noise reduction and information transmission.
Collapse
Affiliation(s)
- Martin Falcke
- Max Delbrück Center for Molecular Medicine, Robert Rössle Str. 10, 13125 Berlin, Germany
| | | |
Collapse
|
42
|
Gonze D, Gérard C, Wacquier B, Woller A, Tosenberger A, Goldbeter A, Dupont G. Modeling-Based Investigation of the Effect of Noise in Cellular Systems. Front Mol Biosci 2018; 5:34. [PMID: 29707543 PMCID: PMC5907451 DOI: 10.3389/fmolb.2018.00034] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2018] [Accepted: 03/26/2018] [Indexed: 12/14/2022] Open
Abstract
Noise is pervasive in cellular biology and inevitably affects the dynamics of cellular processes. Biological systems have developed regulatory mechanisms to ensure robustness with respect to noise or to take advantage of stochasticity. We review here, through a couple of selected examples, some insights on possible robustness factors and constructive roles of noise provided by computational modeling. In particular, we focus on (1) factors that likely contribute to the robustness of oscillatory processes such as the circadian clocks and the cell cycle, (2) how reliable coding/decoding of calcium-mediated signaling could be achieved in presence of noise and, in some cases, enhanced through stochastic resonance, and (3) how embryonic cell differentiation processes can exploit stochasticity to create heterogeneity in a population of identical cells.
Collapse
Affiliation(s)
- Didier Gonze
- Unité de Chronobiologie Théorique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| | - Claude Gérard
- de Duve Institute (LPAD Group), Université Catholique de Louvain, Brussels, Belgium
| | - Benjamin Wacquier
- Unité de Chronobiologie Théorique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| | - Aurore Woller
- Unité de Chronobiologie Théorique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| | - Alen Tosenberger
- Unité de Chronobiologie Théorique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| | - Albert Goldbeter
- Unité de Chronobiologie Théorique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| | - Geneviève Dupont
- Unité de Chronobiologie Théorique, Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium
| |
Collapse
|
43
|
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
|
44
|
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: 34] [Impact Index Per Article: 5.7] [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
|
45
|
Watson U, Jain R, Asthana S, Saini DK. Spatiotemporal Modulation of ERK Activation by GPCRs. INTERNATIONAL REVIEW OF CELL AND MOLECULAR BIOLOGY 2018; 338:111-140. [DOI: 10.1016/bs.ircmb.2018.02.004] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
|
46
|
Ca 2+ signals initiate at immobile IP 3 receptors adjacent to ER-plasma membrane junctions. Nat Commun 2017; 8:1505. [PMID: 29138405 PMCID: PMC5686115 DOI: 10.1038/s41467-017-01644-8] [Citation(s) in RCA: 102] [Impact Index Per Article: 14.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2016] [Accepted: 10/06/2017] [Indexed: 11/08/2022] Open
Abstract
IP3 receptors (IP3Rs) release Ca2+ from the ER when they bind IP3 and Ca2+. The spatial organization of IP3Rs determines both the propagation of Ca2+ signals between IP3Rs and the selective regulation of cellular responses. Here we use gene editing to fluorescently tag endogenous IP3Rs, and super-resolution microscopy to determine the geography of IP3Rs and Ca2+ signals within living cells. We show that native IP3Rs cluster within ER membranes. Most IP3R clusters are mobile, moved by diffusion and microtubule motors. Ca2+ signals are generated by a small population of immobile IP3Rs. These IP3Rs are licensed to respond, but they do not readily mix with mobile IP3Rs. The licensed IP3Rs reside alongside ER-plasma membrane junctions where STIM1, which regulates store-operated Ca2+ entry, accumulates after depletion of Ca2+ stores. IP3Rs tethered close to ER-plasma membrane junctions are licensed to respond and optimally placed to be activated by endogenous IP3 and to regulate Ca2+ entry. IP3 receptors mediate Ca2+ release from the endoplasmic reticulum. Here the authors show that only a small fraction of IP3 receptors initiate Ca2+ signals; these immobile IP3 receptors adjacent to the plasma membrane are optimally placed to control STIM1-dependent Ca2+ entry.
Collapse
|
47
|
Grundmann M, Kostenis E. Temporal Bias: Time-Encoded Dynamic GPCR Signaling. Trends Pharmacol Sci 2017; 38:1110-1124. [PMID: 29074251 DOI: 10.1016/j.tips.2017.09.004] [Citation(s) in RCA: 56] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2017] [Revised: 09/24/2017] [Accepted: 09/25/2017] [Indexed: 12/21/2022]
Abstract
Evidence suggests that cells can time-encode signals for secure transport and perception of information, and it appears that this dynamic signaling is a common principle of nature to code information in time. G-protein-coupled receptor (GPCR) signaling networks are no exception as their composition and signal transduction appear temporally flexible. In this review, we discuss the potential mechanisms by which GPCRs code biological information in time to create 'temporal bias.' We highlight dynamic signaling patterns from the second messenger to the receptor-ligand level and shed light on the dynamics of G-protein cycles, the kinetics of ligand-receptor interaction, and the occurrence of distinct signaling waves within the cell. A dynamic feature such as temporal bias adds to the complexity of GPCR signaling bias and gives rise to the question whether this trait could be exploited to gain control over time-encoded cell physiology.
Collapse
Affiliation(s)
- Manuel Grundmann
- Molecular-, Cellular- and Pharmacobiology Section, Institute for Pharmaceutical Biology, University of Bonn, Nussallee 6, 53115 Bonn, Germany; Kidney Disease Research, Bayer Pharma AG, Aprather Weg 18a, 42113 Wuppertal, Germany
| | - Evi Kostenis
- Molecular-, Cellular- and Pharmacobiology Section, Institute for Pharmaceutical Biology, University of Bonn, Nussallee 6, 53115 Bonn, Germany.
| |
Collapse
|
48
|
Abstract
One driving motivation of systems biology is the search for general principles that govern the design of biological systems. But questions often arise as to what kind of general principles biology could have. Concepts from engineering such as robustness and modularity are indeed becoming a regular way of describing biological systems. Another source of potential general principles is the emerging similarities found in processes across biological hierarchies. In this piece, I describe several emerging cross-hierarchy similarities. Identification of more cross-hierarchy principles, and understanding the implications these convergence have on the construction of biological systems, I believe, present exciting challenges for systems biology in the decades to come.
Collapse
|
49
|
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
|
50
|
Keebler MV, Taylor CW. Endogenous signalling pathways and caged IP 3 evoke Ca 2+ puffs at the same abundant immobile intracellular sites. J Cell Sci 2017; 130:3728-3739. [PMID: 28893841 PMCID: PMC5702060 DOI: 10.1242/jcs.208520] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Accepted: 09/04/2017] [Indexed: 12/18/2022] Open
Abstract
The building blocks of intracellular Ca2+ signals evoked by inositol 1,4,5-trisphosphate receptors (IP3Rs) are Ca2+ puffs, transient focal increases in Ca2+ concentration that reflect the opening of small clusters of IP3Rs. We use total internal reflection fluorescence microscopy and automated analyses to detect Ca2+ puffs evoked by photolysis of caged IP3 or activation of endogenous muscarinic receptors with carbachol in human embryonic kidney 293 cells. Ca2+ puffs evoked by carbachol initiated at an estimated 65±7 sites/cell, and the sites remained immobile for many minutes. Photolysis of caged IP3 evoked Ca2+ puffs at a similar number of sites (100±35). Increasing the carbachol concentration increased the frequency of Ca2+ puffs without unmasking additional Ca2+ release sites. By measuring responses to sequential stimulation with carbachol or photolysed caged IP3, we established that the two stimuli evoked Ca2+ puffs at the same sites. We conclude that IP3-evoked Ca2+ puffs initiate at numerous immobile sites and the sites become more likely to fire as the IP3 concentration increases; there is no evidence that endogenous signalling pathways selectively deliver IP3 to specific sites. Summary: Ca2+ puffs are the building blocks for IP3-evoked Ca2+ signals. Ca2+ puffs evoked by caged IP3 or via endogenous signalling pathways initiate at the same fixed intracellular sites.
Collapse
Affiliation(s)
- Michael V Keebler
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK
| | - Colin W Taylor
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, UK
| |
Collapse
|