A mathematical model of agonist-induced propagation of calcium waves in astrocytes

Analysis of Network Models with Neuron-Astrocyte Interactions

Neural networks, composed of many neurons and governed by complex interactions between them, are a widely accepted formalism for modeling and exploring global dynamics and emergent properties in brain systems. In the past decades, experimental evidence of computationally relevant neuron-astrocyte interactions, as well as the astrocytic modulation of global neural dynamics, have accumulated. These findings motivated advances in computational glioscience and inspired several models integrating mechanisms of neuron-astrocyte interactions into the standard neural network formalism. These models were developed to study, for example, synchronization, information transfer, synaptic plasticity, and hyperexcitability, as well as classification tasks and hardware implementations. We here focus on network models of at least two neurons interacting bidirectionally with at least two astrocytes that include explicitly modeled astrocytic calcium dynamics. In this study, we analyze the evolution of these models and the biophysical, biochemical, cellular, and network mechanisms used to construct them. Based on our analysis, we propose how to systematically describe and categorize interaction schemes between cells in neuron-astrocyte networks. We additionally study the models in view of the existing experimental data and present future perspectives. Our analysis is an important first step towards understanding astrocytic contribution to brain functions. However, more advances are needed to collect comprehensive data about astrocyte morphology and physiology in vivo and to better integrate them in data-driven computational models. Broadening the discussion about theoretical approaches and expanding the computational tools is necessary to better understand astrocytes' roles in brain functions.

Local and CNS-Wide Astrocyte Intracellular Calcium Signaling Attenuation In Vivo with CalExflox Mice

Astrocytes exist throughout the CNS and affect neural circuits and behavior through intracellular Ca2+ signaling. Studying the function(s) of astrocyte Ca2+ signaling has proven difficult because of the paucity of tools to achieve selective attenuation. Based on recent studies, we generated and used male and female knock-in mice for Cre-dependent expression of mCherry-tagged hPMCA2w/b to attenuate astrocyte Ca2+ signaling in genetically defined cells in vivo (CalExflox mice for Calcium Extrusion). We characterized CalExflox mice following local AAV-Cre microinjections into the striatum and found reduced astrocyte Ca2+ signaling (∼90%) accompanied with repetitive self-grooming behavior. We also crossed CalExflox mice to astrocyte-specific Aldh1l1-Cre/ERT2 mice to achieve inducible global CNS-wide Ca2+ signaling attenuation. Within 6 d of induction in the bigenic mice, we observed significantly altered ambulation in the open field, disrupted motor coordination and gait, and premature lethality. Furthermore, with histologic, imaging, and transcriptomic analyses, we identified cellular and molecular alterations in the cerebellum following mCherry-tagged hPMCA2w/b expression. Our data show that expression of mCherry-tagged hPMCA2w/b with CalExflox mice throughout the CNS resulted in substantial attenuation of astrocyte Ca2+ signaling and significant behavioral alterations in adult mice. We interpreted these findings candidly in relation to the ability of CalEx to attenuate astrocyte Ca2+ signaling, with regards to additional mechanistic interpretations of the data, and their relation to past studies that reduced astrocyte Ca2+ signaling throughout the CNS. The data and resources provide complementary ways to interrogate the function(s) of astrocytes in multiple experimental scenarios.SIGNIFICANCE STATEMENT Astrocytes represent a significant fraction of all brain cells and tile the entire central nervous system. Unlike neurons, astrocytes lack propagated electrical signals. Instead, astrocytes are proposed to use diverse and dynamic intracellular Ca2+ signals to communicate with other cells. An open question concerns if and how astrocyte Ca2+ signaling regulates behavior in adult mice. We approached this problem by generating a new transgenic mouse line to achieve inducible astrocyte Ca2+ signaling attenuation in vivo We report our data with this mouse line and we interpret the findings candidly in relation to past studies and within the framework of different mechanistic interpretations.

Chain Modeling of Molecular Communications for Body Area Network

Molecular communications provide an attractive opportunity to precisely regulate biological signaling in nano-medicine applications of body area networks. In this paper, we utilize molecular communication tools to interpret how neural signals are generated in response to external stimuli. First, we propose a chain model of molecular communication system by considering three types of biological signaling through different communication media. Second, communication models of hormonal signaling, Ca 2 + signaling and neural signaling are developed based on existing knowledge. Third, an amplify-and-forward relaying mechanism is proposed to connect different types of signaling. Simulation results demonstrate that the proposed communication system facilitates the information exchange between the neural system and nano-machines, and suggests that proper adjustment can optimize the communication system performance.

Computational Models for Calcium-Mediated Astrocyte Functions

The computational neuroscience field has heavily concentrated on the modeling of neuronal functions, largely ignoring other brain cells, including one type of glial cell, the astrocytes. Despite the short history of modeling astrocytic functions, we were delighted about the hundreds of models developed so far to study the role of astrocytes, most often in calcium dynamics, synchronization, information transfer, and plasticity in vitro, but also in vascular events, hyperexcitability, and homeostasis. Our goal here is to present the state-of-the-art in computational modeling of astrocytes in order to facilitate better understanding of the functions and dynamics of astrocytes in the brain. Due to the large number of models, we concentrated on a hundred models that include biophysical descriptions for calcium signaling and dynamics in astrocytes. We categorized the models into four groups: single astrocyte models, astrocyte network models, neuron-astrocyte synapse models, and neuron-astrocyte network models to ease their use in future modeling projects. We characterized the models based on which earlier models were used for building the models and which type of biological entities were described in the astrocyte models. Features of the models were compared and contrasted so that similarities and differences were more readily apparent. We discovered that most of the models were basically generated from a small set of previously published models with small variations. However, neither citations to all the previous models with similar core structure nor explanations of what was built on top of the previous models were provided, which made it possible, in some cases, to have the same models published several times without an explicit intention to make new predictions about the roles of astrocytes in brain functions. Furthermore, only a few of the models are available online which makes it difficult to reproduce the simulation results and further develop the models. Thus, we would like to emphasize that only via reproducible research are we able to build better computational models for astrocytes, which truly advance science. Our study is the first to characterize in detail the biophysical and biochemical mechanisms that have been modeled for astrocytes.

Reproducibility and Comparability of Computational Models for Astrocyte Calcium Excitability

The scientific community across all disciplines faces the same challenges of ensuring accessibility, reproducibility, and efficient comparability of scientific results. Computational neuroscience is a rapidly developing field, where reproducibility and comparability of research results have gained increasing interest over the past years. As the number of computational models of brain functions is increasing, we chose to address reproducibility using four previously published computational models of astrocyte excitability as an example. Although not conventionally taken into account when modeling neuronal systems, astrocytes have been shown to take part in a variety of in vitro and in vivo phenomena including synaptic transmission. Two of the selected astrocyte models describe spontaneous calcium excitability, and the other two neurotransmitter-evoked calcium excitability. We specifically addressed how well the original simulation results can be reproduced with a reimplementation of the models. Additionally, we studied how well the selected models can be reused and whether they are comparable in other stimulation conditions and research settings. Unexpectedly, we found out that three of the model publications did not give all the necessary information required to reimplement the models. In addition, we were able to reproduce the original results of only one of the models completely based on the information given in the original publications and in the errata. We actually found errors in the equations provided by two of the model publications; after modifying the equations accordingly, the original results were reproduced more accurately. Even though the selected models were developed to describe the same biological event, namely astrocyte calcium excitability, the models behaved quite differently compared to one another. Our findings on a specific set of published astrocyte models stress the importance of proper validation of the models against experimental wet-lab data from astrocytes as well as the careful review process of models. A variety of aspects of model development could be improved, including the presentation of models in publications and databases. Specifically, all necessary mathematical equations, as well as parameter values, initial values of variables, and stimuli used should be given precisely for successful reproduction of scientific results.

Astrocyte-neuron interactions: from experimental research-based models to translational medicine

In this chapter, we review the principal astrocyte functions and the interactions between neurons and astrocytes. We then address how the experimentally observed functions have been verified in computational models and review recent experimental literature on astrocyte-neuron interactions. Benefits of computational neuroscience work are highlighted through selected studies with neurons and astrocytes by analyzing the existing models qualitatively and assessing the relevance of these models to experimental data. Common strategies to mathematical modeling and computer simulation in neuroscience are summarized for the nontechnical reader. The astrocyte-neuron interactions are then further illustrated by examples of some neurological and neurodegenerative diseases, where the miscommunication between glia and neurons is found to be increasingly important.

Combined computational and experimental approaches to understanding the Ca(2+) regulatory network in neurons

Ca(2+) is a ubiquitous signaling ion that regulates a variety of neuronal functions by binding to and altering the state of effector proteins. Spatial relationships and temporal dynamics of Ca(2+) elevations determine many cellular responses of neurons to chemical and electrical stimulation. There is a wealth of information regarding the properties and distribution of Ca(2+) channels, pumps, exchangers, and buffers that participate in Ca(2+) regulation. At the same time, new imaging techniques permit characterization of evoked Ca(2+) signals with increasing spatial and temporal resolution. However, understanding the mechanistic link between functional properties of Ca(2+) handling proteins and the stimulus-evoked Ca(2+) signals they orchestrate requires consideration of the way Ca(2+) handling mechanisms operate together as a system in native cells. A wide array of biophysical modeling approaches is available for studying this problem and can be used in a variety of ways. Models can be useful to explain the behavior of complex systems, to evaluate the role of individual Ca(2+) handling mechanisms, to extract valuable parameters, and to generate predictions that can be validated experimentally. In this review, we discuss recent advances in understanding the underlying mechanisms of Ca(2+) signaling in neurons via mathematical modeling. We emphasize the value of developing realistic models based on experimentally validated descriptions of Ca(2+) transport and buffering that can be tested and refined through new experiments to develop increasingly accurate biophysical descriptions of Ca(2+) signaling in neurons.

Recent advances in brain physiology and cognitive processing

The discovery of participation of astrocytes as active elements in glutamatergic tripartite synapses (composed by functional units of two neurons and one astrocyte) has led to the construction of models of cognitive functioning in the human brain, focusing on associative learning, sensory integration, conscious processing and memory formation/retrieval. We have modelled human cognitive functions by means of an ensemble of functional units (tripartite synapses) connected by gap junctions that link distributed astrocytes, allowing the formation of intra- and intercellular calcium waves that putatively mediate large-scale cognitive information processing. The model contains a diagram of molecular mechanisms present in tripartite synapses and contributes to explain the physiological bases of cognitive functions. It can be potentially expanded to explain emotional functions and psychiatric phenomena.

Astrocytes and human cognition: modeling information integration and modulation of neuronal activity

Recent research focusing on the participation of astrocytes in glutamatergic tripartite synapses has revealed mechanisms that support cognitive functions common to human and other mammalian species, such as learning, perception, conscious integration, memory formation/retrieval and the control of voluntary behavior. Astrocytes can modulate neuronal activity by means of release of glutamate, d-serine, adenosine triphosphate and other signaling molecules, contributing to sustain, reinforce or depress pre- and post-synaptic membranes. We review molecular mechanisms present in tripartite synapses and model the cognitive role of astrocytes. Single protoplasmic astrocytes operate as a "Local Hub", integrating information patterns from neuronal and glial populations. Two mechanisms, here modeled as the "domino" and "carousel" effects, contribute to the formation of intercellular calcium waves. As waves propagate through gap junctions and reach other types of astrocytes (interlaminar, polarized, fibrous and varicose projection), the active astroglial network functions as a "Master Hub" that integrates results of distributed processing from several brain areas and supports conscious states. Response of this network would define the effect exerted on neuronal plasticity (membrane potentiation or depression), behavior and psychosomatic processes. Theoretical results of our modeling can contribute to the development of new experimental research programs to test cognitive functions of astrocytes.

Generation and propagation of gastric slow waves

1. Mechanisms underlying the generation and propagation of gastrointestinal slow wave depolarizations have long been controversial. The present review aims to collate present knowledge on this subject with specific reference to slow waves in gastric smooth muscle. 2. At present, there is strong agreement that interstitial cells of Cajal (ICC) are the pacemaker cells that generate slow waves. What has been less clear is the relative role of primary types of ICC, including the network in the myenteric plexus (ICC-MY) and the intramuscular network (ICC-IM). It is concluded that both ICC-MY and ICC-IM are likely to serve a major role in slow wave generation and propagation. 3. There has been long-standing controversy as to how slow waves 'propagate' circumferentially and down the gastrointestinal tract. Two mechanisms have been proposed, one being action potential (AP)-like conduction and the other phase wave-based 'propagation' resulting from an interaction of coupled oscillators. Studies made on single bundle gastric strips indicate that both mechanisms apply with relative dominance depending on conditions; the phase wave mechanism is dominant under circumstances of rhythmically generating slow waves and the AP-like propagation is dominant when the system is perturbed. 4. The phase wave mechanism (termed Ca(2+) phase wave) uses cyclical Ca(2+) release as the oscillator, with coupling between oscillators mediated by several factors, including: (i) store-induced depolarization; (ii) resultant electrical current flow/depolarization through the pacemaker cell network; and (iii) depolarization-induced increase in excitability of downstream Ca(2+) stores. An analogy is provided by pendulums in an array coupled together by a network of springs. These, when randomly activated, entrain to swing at the same frequency but with a relative delay along the row giving the impression of a propagating wave. 5. The AP-like mechanism (termed voltage-accelerated Ca(2+) wave) propagates sequentially like a conducting AP. However, it is different in that it depends on regenerative store Ca(2+) release and resultant depolarization rather than regenerative activation of voltage-dependent channels in the cell membrane. 6. The applicability of these mechanisms to describing propagation in large intact gastrointestinal tissues, where voltage-dependent Ca(2+) entry is also likely to be functional, is discussed.

On the role of synchrony for neuron-astrocyte interactions and perceptual conscious processing

Recent research on brain correlates of cognitive processes revealed the occurrence of global synchronization during conscious processing of sensory stimuli. In spite of technological progress in brain imaging, an explanation of the computational role of synchrony is still a highly controversial issue. In this study, we depart from an analysis of the usage of blood-oxygen-level-dependent functional magnetic resonance imaging for the study of cognitive processing, leading to the identification of evoked local field potentials as the vehicle for sensory patterns that compose conscious episodes. Assuming the "astrocentric hypothesis" formulated by James M. Robertson (astrocytes being the final stage of conscious processing), we propose that the role of global synchrony in perceptual conscious processing is to induce the transfer of information patterns embodied in local field potentials to astrocytic calcium waves, further suggesting that these waves are responsible for the "binding" of spatially distributed patterns into unitary conscious episodes.

What is the role of astrocyte calcium in neurophysiology?

Astrocytes comprise approximately half of the volume of the adult mammalian brain and are the primary neuronal structural and trophic supportive elements. Astrocytes are organized into distinct nonoverlapping domains and extend elaborate and dense fine processes that interact intimately with synapses and cerebrovasculature. The recognition in the mid 1990s that astrocytes undergo elevations in intracellular calcium concentration following activation of G protein-coupled receptors by synaptically released neurotransmitters demonstrated not only that astrocytes display a form of excitability but also that astrocytes may be active participants in brain information processing. The roles that astrocytic calcium elevations play in neurophysiology and especially in modulation of neuronal activity have been intensely researched in recent years. This review will summarize the current understanding of the function of astrocytic calcium signaling in neurophysiological processes and discuss areas where the role of astrocytes remains controversial and will therefore benefit from further study.

Selective stimulation of astrocyte calcium in situ does not affect neuronal excitatory synaptic activity

Astrocytes are considered the third component of the synapse, responding to neurotransmitter release from synaptic terminals and releasing gliotransmitters--including glutamate--in a Ca(2+)-dependent manner to affect neuronal synaptic activity. Many studies reporting astrocyte-driven neuronal activity have evoked astrocyte Ca(2+) increases by application of endogenous ligands that directly activate neuronal receptors, making astrocyte contribution to neuronal effect(s) difficult to determine. We have made transgenic mice that express a Gq-coupled receptor only in astrocytes to evoke astrocyte Ca(2+) increases using an agonist that does not bind endogenous receptors in brain. By recording from CA1 pyramidal cells in acute hippocampal slices from these mice, we demonstrate that widespread Ca(2+) elevations in 80%-90% of stratum radiatum astrocytes do not increase neuronal Ca(2+), produce neuronal slow inward currents, or affect excitatory synaptic activity. Our findings call into question the developing consensus that Ca(2+)-dependent glutamate release by astrocytes directly affects neuronal synaptic activity in situ.

Pacemaking through Ca2+ stores interacting as coupled oscillators via membrane depolarization

This study presents an investigation of pacemaker mechanisms underlying lymphatic vasomotion. We tested the hypothesis that active inositol 1,4,5-trisphosphate receptor (IP(3)R)-operated Ca(2+) stores interact as coupled oscillators to produce near-synchronous Ca(2+) release events and associated pacemaker potentials, this driving action potentials and constrictions of lymphatic smooth muscle. Application of endothelin 1 (ET-1), an agonist known to enhance synthesis of IP(3), to quiescent lymphatic smooth muscle syncytia first enhanced spontaneous Ca(2+) transients and/or intracellular Ca(2+) waves. Larger near-synchronous Ca(2+) transients then occurred leading to global synchronous Ca(2+) transients associated with action potentials and resultant vasomotion. In contrast, blockade of L-type Ca(2+) channels with nifedipine prevented ET-1 from inducing near-synchronous Ca(2+) transients and resultant action potentials, leaving only asynchronous Ca(2+) transients and local Ca(2+) waves. These data were well simulated by a model of lymphatic smooth muscle with: 1), oscillatory Ca(2+) release from IP(3)R-operated Ca(2+) stores, which causes depolarization; 2), L-type Ca(2+) channels; and 3), gap junctions between cells. Stimulation of the stores caused global pacemaker activity through coupled oscillator-based entrainment of the stores. Membrane potential changes and positive feedback by L-type Ca(2+) channels to produce more store activity were fundamental to this process providing long-range electrochemical coupling between the Ca(2+) store oscillators. We conclude that lymphatic pacemaking is mediated by coupled oscillator-based interactions between active Ca(2+) stores. These are weakly coupled by inter- and intracellular diffusion of store activators and strongly coupled by membrane potential. Ca(2+) store-based pacemaking is predicted for cellular systems where: 1), oscillatory Ca(2+) release induces depolarization; 2), membrane depolarization provides positive feedback to induce further store Ca(2+) release; and 3), cells are interconnected. These conditions are met in a surprisingly large number of cellular systems including gastrointestinal, lymphatic, urethral, and vascular tissues, and in heart pacemaker cells.

A quantitative model of purinergic junctional transmission of calcium waves in astrocyte networks

A principal means of transmitting intracellular calcium (Ca2+) waves at junctions between astrocytes involves the release of the chemical transmitter adenosine triphosphate (ATP). A model of this process is presented in which activation of purinergic P2Y receptors by ATP triggers the release of ATP, in an autocrine manner, as well as concomitantly increasing intracellular Ca2+. The dependence of the temporal characteristics of the Ca2+ wave are shown to critically depend on the dissociation constant (K(R)) for ATP binding to the P2Y receptor type. Incorporating this model astrocyte into networks of these cells successfully accounts for many of the properties of propagating Ca2+ waves, such as the dependence of velocity on the type of P2Y receptor and the time-lag of the Ca2+ wave behind the ATP wave. In addition, the conditions under which Ca2+ waves may jump from one set of astrocytes across an astrocyte-free lane to another set of astrocytes are quantitatively accounted for by the model. The properties of purinergic transmission at astrocyte junctions may determine many of the characteristics of Ca2+ propagation in networks of these cells.

Role of voltage-dependent modulation of store Ca2+ release in synchronization of Ca2+ oscillations

Slow waves are rhythmic depolarizations that underlie mechanical activity of many smooth muscles. Slow waves result through rhythmic Ca(2+) release from intracellular Ca(2+) stores through inositol 1,4,5-trisphosphate (IP(3)) sensitive receptors and Ca(2+)-induced Ca(2+) release. Ca(2+) oscillations are transformed into membrane depolarizations by generation of a Ca(2+)-activated inward current. Importantly, the store Ca(2+) oscillations that underlie slow waves are entrained across many cells over large distances. It has been shown that IP(3) receptor-mediated Ca(2+) release is enhanced by membrane depolarization. Previous studies have implicated diffusion of Ca(2+) or the second messenger IP(3) across gap junctions in synchronization of Ca(2+) oscillations. In this study, a novel mechanism of Ca(2+) store entrainment through depolarization-induced IP(3) receptor-mediated Ca(2+) release is investigated. This mechanism is significantly different from chemical coupling-based mechanisms, as membrane potential has a coupling effect over distances several orders of magnitude greater than either diffusion of Ca(2+) or IP(3) through gap junctions. It is shown that electrical coupling acting through voltage-dependent modulation of store Ca(2+) release is able to synchronize oscillations of cells even when cells are widely separated and have different intrinsic frequencies of oscillation.

Intercellular calcium signalling in cultured renal epithelia: a theoretical study of synchronization mode and pacemaker activity

We investigate a two-dimensional lattice model representation of intercellular Ca2+ signalling in a population of epithelial cells coupled by gap junctions. The model is based on and compared with Ca2+ imaging data from globally bradykinin-stimulated MDCK-I (Madin-Darby canine kidney)-I cell layers. We study large-scale synchronization of relevance to our laboratory experiments. The system is found to express a wealth of dynamics, including quasiperiodic, chaotic and multiply-periodic behaviour for intermediate couplings. We take a particular interest in understanding the role of "pacemaker cells" in the synchronization process. It has been hypothesized that a few highly hormone-sensitive cells control the collective frequency of oscillation, which is close to the natural frequencies (without coupling) of these cells. The model behaviour is consistent with the conjectures of the pacemaker cell hypothesis near the critical coupling where the cells lock onto a single frequency. However, the simulations predict that the frequency in globally connected systems decreases with increasing coupling. It is found that a pacemaker is not defined by its natural frequency alone, but that other intrinsic or local factors must be considered. Inclusion of partly sensitized cells that do not oscillate autonomously in the cell layer increases the coupling necessary for global synchronization. For not excessively high coupling, these cells oscillate irregularly and with distinctive lower frequencies. In summary, the present study shows that the frequency of synchronized oscillations is not dictated by one or few fast-responding cells. The collective frequency is the result of a two-way communication between the phase-advanced pacemaker and its environment.

Ca2+ phase waves: a basis for cellular pacemaking and long-range synchronicity in the guinea-pig gastric pylorus

Ca2+ imaging and multiple microelectrode recording procedures were used to investigate a slow wave-like electrical rhythmicity in single bundle strips from the circular muscle layer of the guinea-pig gastric pylorus. The 'slow waves' (SWs) consisted of a pacemaker and regenerative component, with both potentials composed of more elementary events variously termed spontaneous transient depolarizations (STDs) or unitary potentials. STDs and SW pacemaker and regenerative potentials exhibited associated local and distributed Ca2+ transients, respectively. Ca2+ transients were often larger in cellular regions that exhibited higher basal Ca2+ indicator-associated fluorescence, typical of regions likely to contain intramuscular interstitial cells of Cajal (ICCIM). The emergence of rhythmicity arose through entrainment of STDs resulting in pacemaker Ca2+ transients and potentials, events that exhibited considerable spatial synchronicity. Application of ACh to strips exhibiting weak rhythmicity caused marked enhancement of SW synchronicity. SWs and underlying Ca2+ increases exhibited very high 'apparent conduction velocities' ('CVs') orders of magnitude greater than for sequentially conducting Ca2+ waves. Central interruption of either intercellular connectivity or inositol 1,4,5-trisphosphate receptor (IP3R)-mediated store Ca2+ release in strips caused SWs at the two ends to run independently of each other, consistent with a coupled oscillator-based mechanism. Central inhibition of stores required much wider regions of blockade than inhibition of connectivity indicating that stores were voltage-coupled. Simulations, made using a conventional store array model but now including depolarization coupled to IP3R-mediated Ca2+ release, predicted the experimental findings. The linkage between membrane voltage and Ca2+ release provides a means for stores to interact as strongly coupled oscillators, resulting in the emergence of Ca2+ phase waves and associated pacemaker potentials. This distributed pacemaker triggers regenerative Ca2+ release and resultant SWs.

Selection of intracellular calcium patterns in a model with clustered Ca2+ release channels

A two-dimensional model is proposed for intracellular Ca2+ waves, which incorporates both the discrete nature of Ca2+ release sites in the endoplasmic reticulum membrane and the stochastic dynamics of the clustered inositol 1,4,5-triphosphate (IP3) receptors. Depending on the Ca2+ diffusion coefficient and concentration of IP3, various spontaneous Ca2+ patterns, such as calcium puffs, local waves, abortive waves, global oscillation, and tide waves, can be observed. We further investigate the speed of the global waves as a function of the IP3 concentration and the Ca2+ diffusion coefficient and under what conditions the spatially averaged Ca2+ response can be described by a simple set of ordinary differential equations.

Astrocytes in adult rat brain express type 2 inositol 1,4,5-trisphosphate receptors

Astrocytes respond to neuronal activity by propagating Ca(2+) waves elicited through the inositol 1,4,5-trisphosphate pathway. We have previously shown that wave propagation is supported by specialized Ca(2+) release sites, where a number of proteins, including inositol 1,4,5-trisphosphate receptors (IP(3)R), occur together in patches. The specific IP(3)R isoform expressed by astrocytes in situ in rat brain is unknown. In the present report, we use isoform-specific antibodies to localize immunohistochemically the IP(3)R subtype expressed in astrocytes in rat brain sections. Astrocytes were identified using antibodies against the astrocyte-specific markers, S-100 beta, or GFAP. Dual indirect immunohistochemistry showed that astrocytes in all regions of adult rat brain express only IP(3)R2. High-resolution analysis showed that hippocampal astrocytes are endowed with a highly branched network of processes that bear fine hair-like extensions containing punctate patches of IP(3)R2 staining in intimate contact with synapses. Such an organization is reminiscent of signaling microdomains found in cultured glial cells. Similarly, Bergmann glial cell processes in the cerebellum also contained fine hair-like processes containing IP(3)R2 staining. The IP(3)R2-containing fine terminal branches of astrocyte processes in both brain regions were found juxtaposed to presynaptic terminals containing synaptophysin as well as PSD 95-containing postsynaptic densities. Corpus callosum astrocytes had an elongated morphology with IP(3)R2 studded processes extending along fiber tracts. Our data suggest that PLC-mediated Ca(2+) signaling in astrocytes in rat brain occurs predominantly through IP(3)R2 ion channels. Furthermore, the anatomical arrangement of the terminal astrocytic branches containing IP(3)R2 ensheathing synapses is ideal for supporting glial monitoring of neuronal activity.

Sparks and puffs in oligodendrocyte progenitors: cross talk between ryanodine receptors and inositol trisphosphate receptors

Investigating how calcium release from the endoplasmic reticulum (ER) is triggered and coordinated is crucial to our understanding of how oligodendrocyte progenitor cells (OPs) develop into myelinating cells. Sparks and puffs represent highly localized Ca(2+) release from the ER through ryanodine receptors (RyRs) and inositol trisphosphate receptors (IP(3)Rs), respectively. To study whether sparks or puffs trigger Ca(2+) waves in OPs, we performed rapid high-resolution line scan recordings in fluo-4-loaded OP processes. We found spontaneous and evoked sparks and puffs, and we have identified functional cross talk between IP(3)Rs and RyRs. Local events evoked using the IP(3)-linked agonist methacholine (MeCh) showed significantly different morphology compared with events evoked using the caffeine analog 3,7-dimethyl-1-propargylxanthine (DMPX). Pretreatment with MeCh potentiated DMPX-evoked events, whereas inhibition of RyRs potentiated events evoked by low concentrations of MeCh. Furthermore, activation of IP(3)Rs but not RyRs was critical for Ca(2+) wave initiation. Using immunocytochemistry, we show OPs express the specific Ca(2+) release channel subtypes RyR3 and IP(3)R2 in patches along OP processes. RyRs are coexpressed with IP(3)Rs in some patches, but IP(3)Rs are also found alone. This differential distribution pattern may underlie the differences in local and global Ca(2+) signals mediated by these two receptors. Thus, in OPs, interactions between IP(3)Rs and RyRs determine the spatial and temporal characteristics of calcium signaling, from microdomains to intracellular waves.

Müller cell Ca2+ waves evoked by purinergic receptor agonists in slices of rat retina

We have measured agonist evoked Ca2+ waves in Müller cells in situ within freshly isolated retinal slices. Using an eye cup dye loading procedure we were able to preferentially fill Müller glial cells in retinal slices with calcium green. Fluorescence microscopy revealed that bath perfusion of slices with purinergic agonists elicits Ca2+ waves in Müller cells, which propagate along their processes. These Ca2+ signals were insensitive to tetrodotoxin (TTX, 1.0 microM) pretreatment. Cells were readily identified as Müller cells by their unique morphology and by subsequent immunocytochemical labeling with glial fibrillary acidic protein antibodies. While cells never exhibited spontaneous Ca2+ oscillations, purinoreceptor agonists, ATP, 2 MeSATP, ADP, 2 MeSADP, and adenosine readily elicited Ca2+ waves. These waves persisted in the absence of [Ca2+]o but were abolished by thapsigargin pretreatment, suggesting that the purinergic agonists tested act by releasing Ca2+ from intracellular Ca2+ stores. The rank order of potency of different purines and pyrimidines for inducing Ca2+ signals was 2 MeSATP = 2MeSADP > ADP > ATP >> alphabetameATP = uridine triphosphate (UTP) > uridine diphosphate (UDP). The Ca2+ signals evoked by ATP, ADP, and 2 MeSATP were inhibited by reactive blue (100 microM) and suramin (200 microM), and the adenosine induced signals were abolished only by 3,7-dimethyl-1-propargylxanthine (200 microM) and not by 1,3-dipropyl-8-(2-amino-4-chlorophenyl)-xanthine) or 8-cyclopentyl-1,3-dipropylxanthine at the same concentration. Based on these pharmacological characteristics and the dose-response relationships for ATP, 2 MeSATP, 2 MeSADP, ADP, and adenosine, we concluded that Müller cells express the P1A2 and P2Y1 subtypes of purinoceptors. Analysis of Ca2+ responses showed that, similar to glial cells in culture, wave propagation occurred by regenerative amplification at specialized Ca2+ release sites (wave amplification sites), where the rate of Ca2+ release was significantly enhanced. These data suggest that Müller cells in the retina may participate in signaling, and this may serve as an extra-neuronal signaling pathway.

Lymphatic vasomotion

1. Experimental findings in the past decade have greatly advanced present understanding of electrical/mechanical rhythmicities in smooth muscle, including vasomotion. Lymphatic vessels show strong vasomotor activity and have provided a key experimental system to study these processes. 2. Evidence from lymphatic vessels, blood vessels and other smooth muscles indicates that rhythmical contractions arise through a Ca2+ store-controlled pacemaker mechanism, which can function to cause smooth muscle constriction. 3. Such a model fits with observations that vasomotion can be near synchronous over large vessel lengths involving many cells. 4. The alternative interpretation that smooth muscle rhythmicities are generated by a cardiac-like electrical pacemaker mechanism has not been substantiated in any smooth muscle preparation under normal physiological conditions. However, elements of this latter mechanism are likely to be present at least in some smooth muscles, serving to modulate pacemaking.

An image-based model of calcium waves in differentiated neuroblastoma cells

Calcium waves produced by bradykinin-induced inositol-1,4, 5-trisphosphate (InsP(3))-mediated release from endoplasmic reticulum (ER) have been imaged in N1E-115 neuroblastoma cells. A model of this process was built using the "virtual cell," a general computational system for integrating experimental image, biochemical, and electrophysiological data. The model geometry was based on a cell for which the calcium wave had been experimentally recorded. The distributions of the relevant cellular components [InsP(3) receptor (InsP(3)R)], sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA) pumps, bradykinin receptors, and ER] were based on 3D confocal immunofluorescence images. Wherever possible, known biochemical and electrophysiological data were used to constrain the model. The simulation closely matched the spatial and temporal characteristics of the experimental calcium wave. Predictions on different patterns of calcium signals after InsP(3) uncaging or for different cell geometries were confirmed experimentally, thus helping to validate the model. Models in which the spatial distributions of key components are altered suggest that initiation of the wave in the center of the neurite derives from an interplay of soma-biased ER distribution and InsP(3) generation biased toward the neurite. Simulations demonstrate that mobile buffers (like the indicator fura-2) significantly delay initiation and lower the amplitude of the wave. Analysis of the role played by calcium diffusion indicated that the speed of the wave is only slightly dependent on the ability of calcium to diffuse to and activate neighboring InsP(3) receptor sites.

Role of calcium stores and membrane voltage in the generation of slow wave action potentials in guinea-pig gastric pylorus

1. Intracellular recordings made in single bundle strips of a visceral smooth muscle revealed rhythmic spontaneous membrane depolarizations termed slow waves (SWs). These exhibited 'pacemaker' and 'regenerative' components composed of summations of more elementary events termed spontaneous transient depolarizations (STDs). 2. STDs and SWs persisted in the presence of tetrodotoxin, nifedipine and ryanodine, and upon brief exposure to Ca2+-free Cd2+-containing solutions; they were enhanced by ACh and blocked by BAPTA AM, cyclopiazonic acid and caffeine. 3. SWs were also inhibited in heparin-loaded strips. SWs were observed over a wide range of membrane potentials (e.g. -80 to -45 mV) with increased frequencies at more depolarized potentials. 4. Regular spontaneous SW activity in this preparation began after 1-3 h superfusion of the tissue with physiological saline following the dissection procedure. Membrane depolarization applied before the onset of this activity induced bursts of STD-like events (termed the 'initial' response) which, when larger than threshold levels initiated regenerative responses. The combined initial-regenerative waveform was termed the SW-like action potential. 5. Voltage-induced responses exhibited large variable latencies (typical range 0.3-4 s), refractory periods of approximately 11 s and a pharmacology that was indistinguishable from those of STDs and spontaneous SWs. 6. The data indicate that SWs arise through more elementary inositol 1,4,5-trisphosphate (IP3) receptor-induced Ca2+ release events which rhythmically synchronize to trigger regenerative Ca2+ release and induce inward current across the plasmalemma. The finding that action potentials, which were indistinguishable from SWs, could be evoked by depolarization suggests that membrane potential modulates IP3 production. Voltage feedback on intracellular IP3-sensitive Ca2+ release is likely to have a major influence on the generation and propagation of SWs.

The local control of cytosolic Ca2+ as a propagator of CNS communication--integration of mitochondrial transport mechanisms and cellular responses

Ca2+ signals propagate in wave form along individual cells of the central nervous system (CNS) and through networks of connected cells of neuronal and multiple glial cell types. In order for wave fronts to convey information, signaling mechanisms are required that allow waves to propagate reproducibly and without decrement in signal strength over long distances. CNS Ca2+ waves are under specific integrated local control, made possible by interactions at local subcellular microdomains between endoplasmic reticulum and mitochondria. Active mitochondria located near the mouth of inositol trisphosphate receptor (InsP3R) channel clusters in glia take up Ca2+, which may prevent a buildup of Ca2+ around the InsP3R channel, thereby decreasing the rate of Ca2+-induced receptor inactivation, and prolonging channel open time. Mitochondria may amplify InsP3-dependent Ca2+ signals by a transient permeability transition in response to Ca2+ uptake into the mitochondrion. Other evidence suggests privileged access into mitochondria for Ca2+ entering neurons by glutamatergic receptor channels. This enables specific signal modulation as the Ca2+ wave is propagated into neurons, such that mitochondria located close to glutamate channels can prolong the neuronal cytosolic response time by successive uptake and release of Ca2+. Disruption of mitochondrial function deregulates the ability of CNS-derived cells to undergo normal Ca2+ signaling and wave propagation.

Images of Ca2+ flux in astrocytes: evidence for spatially distinct sites of Ca2+ release and uptake

In this study, we have developed a mathematical method to derive the Ca2+ fluxes underlying agonist-evoked Ca2+ waves in cultured rat cortical astrocytes. Astrocytes were stimulated with norepinephrine (100 nM) to evoke Ca2+ waves, which were recorded by measuring Fluo-3 fluorescence changes with high spatial and temporal resolution. Normalized fluorescence (delta F/F) was analyzed in discrete cellular spaces in a series of successive slices along the length of the cell. From these data, Ca2+ flux was then calculated using a one dimensional reaction-diffusion equation which utilizes the temporal and spatial derivatives of the fluorescence data and the diffusion coefficient of Ca2+ in the cytosol. This method identified distinct sites of positive flux (Ca2+ release into the cytosol) and of negative flux (Ca2+ removal from cytosol) and showed that in astrocytes, sites of Ca2+ release from stores regularly alternate with sites of Ca2+ removal from the cytosol. Cross correlation analysis of the two distribution patterns gave positive correlation at 2 microns out of phase and a negative correlation in phase. Thapsigargin-induced Ca2+ waves were analyzed to determine if the negative flux was due to Ca2+ uptake via thapsigargin-sensitive Ca2+ pumps. Negative flux sites were still found under these conditions, suggesting that multiple mechanisms of Ca2+ removal from the cytosol may contribute to negative flux sites. This method of calculation of flux may serve as a means to describe the distribution of functional ion channels and pumps participating in cellular Ca2+ signalling.

Characterization of ryanodine receptors in oligodendrocytes, type 2 astrocytes, and O-2A progenitors

In this study we have investigated the expression of ryanodine receptors (RyRs), and the ability of caffeine to evoke RyR-mediated elevation of intracellular Ca2+ levels ([Ca2+]i) in glial cells of the oligodendrocyte/type 2 astrocyte lineage. Immunocytochemistry with specific antibodies identified ryanodine receptors in cultured oligodendrocytes, type 2 astrocytes, and O-2A progenitor cells, at high levels in the perinuclear region and in a variegated pattern along processes. Glia acutely isolated from rat brain and in aldehydefixed sections of cortex were similarly found to express RyRs. Caffeine (5-50 mM) caused an increase in [Ca2+]i in most cultured type 2 astrocytes and in 50% of oligodendrocytes. Responses elicited by caffeine were inhibited by pretreatment with ryanodine (10 microM) or thapsigargin (1 microM), and the peak response was unaffected by removal of [Ca2+]o. O-2A progenitor cells, in contrast, were largely unresponsive to caffeine treatment. Pretreatment with kainate (200 microM) to activate Ca2+ entry increased the magnitude of caffeine-evoked [Ca2+]i elevations in type 2 astrocytes and oligodendrocytes, and caused caffeine to activate responses in a significant proportion of previously non-responding O-2A progenitors. In both type 2 astrocytes and oligodendrocytes, caffeine evoked Ca2+ changes which propagated as wavefronts from several initiation sites. These wave amplification sites were characterized by significantly higher local Ca2+ release kinetics. Our results indicate that several glial cell types express RyRs, and that their functionality differs within different cell types of the oligodendrocyte lineage. In addition, ionotropic glutamate receptor activation fills the caffeine-sensitive Ca2+ stores in these cells.

Calcium waves in a model with a random spatially discrete distribution of Ca2+ release sites

We study the propagation of intracellular calcium waves in a model that features Ca2+ release from discrete sites in the endoplasmic reticulum membrane and random spatial distribution of these sites. The results of our simulations qualitatively reproduce the experimentally observed behavior of the waves. When the level of the channel activator inositol trisphosphate is low, the wave undergoes fragmentation and eventually vanishes at a finite distance from the region of initiation, a phenomenon we refer to as an abortive wave. With increasing activator concentration, the mean distance of propagation increases. Above a critical level of activator, the wave becomes stable. We show that the heterogeneous distribution of Ca2+ channels is the cause of this phenomenon.

High density distribution of endoplasmic reticulum proteins and mitochondria at specialized Ca2+ release sites in oligodendrocyte processes

In oligodendrocyte processes, methacholine-evoked Ca2+ waves propagate via regions of specialized Ca2+ release kinetics (wave amplification sites) at which the amplitude and rate of rise of local Ca2+ signals are markedly higher than in surrounding areas (Simpson, P. B., and Russell, J. T. (1996) J. Biol. Chem. 271, 33493-33501). In the present study we have examined the effects of other phosphoinositide-coupled agonists on Ca2+ in these cells, and the structural specializations underlying regenerative wave amplification sites. Both bradykinin and norepinephrine evoke Ca2+ waves, which initiate at the same loci and propagate through the cell body and multiple processes via identical wave amplification sites. Antibodies against type 2 inositol 1,4,5-trisphosphate receptors (InsP3R2) and calreticulin identify expression of these proteins in oligodendrocyte membranes in Western blots. Immunocytochemistry followed by high resolution fluorescence microscopy revealed that both InsP3R2 and calreticulin are expressed in high intensity patches along processes. Cross-correlation analysis of the profiles of local Ca2+ release kinetics during a Ca2+ wave and immunofluorescence for these proteins along cellular processes showed that the domains of high endoplasmic reticulum protein expression correspond closely to wave amplification sites. Staining cells with the mitochondrial dye, MitoTracker(R), showed that mitochondria are only found in intimate association with these sites possessing high density endoplasmic reticulum proteins, and they remain in the same locations over relatively long periods of time. It appears, therefore, that multiple specializations are found at domains of elevated Ca2+ release in oligodendrocyte processes, including high levels of calreticulin, InsP3R2 Ca2+ release channels, and mitochondria.

Role of sarcoplasmic/endoplasmic-reticulum Ca2+-ATPases in mediating Ca2+ waves and local Ca2+-release microdomains in cultured glia

We have characterized the sarcoplasmic-endoplasmic reticulum Ca2+-ATPase (SERCA) pumps in cultured rat cortical type-1 astrocytes, type-2 astrocytes and oligodendrocytes. Perfusion with 10 microM cyclopiazonic acid (CPA) or 1 microM thapsigargin evoked a large and persistent elevation in cytosolic [Ca2+] in normal Ca2+-containing medium and a small and transient increase in nominally Ca2+-free medium. Subtraction of the response in Ca2+-free medium from that in the control revealed a slow-onset Ca2+-entry response to SERCA inhibition, which began after most of the store depletion had occurred. Thapsigargin- and CPA-induced responses propagated as Ca2+ waves, which began in several distinct cellular sites and travelled throughout the cell and through nearby cells, in confluent cultures. Propagation was supported by specialized Ca2+-release sites where the amplitude of the response was significantly higher and the rate of rise steeper. Such higher Ca2+-release kinetics were observed at these sites during Ins(1,4,5)P3-mediated Ca2+ waves in the same cells. Fluorescently tagged thapsigargin labelled SERCA pumps throughout glial cell bodies and processes. In oligodendrocyte processes, multiple domains with elevated SERCA staining were always associated with mitochondria. Our results are consistent with a model in which only a single Ca2+ store, expressing Ins(1,4,5)P3 receptors and SERCAs sensitive to both thapsigargin and CPA, is present in rat cortical glia, and indicate that inhibition of SERCA activates both Ca2+ release as a wavefront and Ca2+ entry via store-operated channels. The spatial relationship between SERCAs and mitochondria is likely to be important for regulating microdomains of elevated Ca2+-release kinetics.

Mitochondria support inositol 1,4,5-trisphosphate-mediated Ca2+ waves in cultured oligodendrocytes

We have examined the spatial and temporal nature of Ca2+ signals activated via the phosphoinositide pathway in oligodendrocytes and the cellular specializations underlying oligodendrocyte Ca2+ response characteristics. Cultured cortical oligodendrocytes were incubated with fluo 3 or fura 2, and digital video fluorescence microscopy was used to study the effect of methacholine on [Ca2+]i. Single peaks, oscillations, and steady-state plateau [Ca2+]i elevations were evoked by increasing agonist concentration. The peaks and oscillations were found to be Ca2+ wave fronts, which propagate via distinct amplification regions in the cell where the kinetics of Ca2+ release (amplitude and rate of rise of response) are elevated. Staining with 5,5',6,6'-tetrachloro-1,1',3,3'-tetraethylbenzimidazolecarbocyanine++ + iodide (JC-1) and 3,3'-dihexyloxacarbocyanine iodide revealed that mitochondria are found in groups of three or more in oligodendrocyte processes and that the groups are distributed with considerable distance separating them. Cross-correlation analysis showed a high degree of correlation between sites where mitochondria are present and peaks in the amplitude and rate of rise of the Ca2+ response. Intramitochondrial Ca2+ concentration, measured using rhod 2, increased upon treatment with methacholine. Methacholine also evoked a rapid change in mitochondrial membrane potential as measured by the J-aggregate fluorescence of JC-1. Pretreatment with the mitochondrial inhibitors carbonyl cyanide p-(trifluoromethoxy)phenylhydrazone (1 microM, 2 min) or antimycin (2 microg/ml, 2 min) altered the methacholine-evoked Ca2+ response in most cells studied, responses being either markedly potentiated or inhibited. The results of this study demonstrate that stimulation of phosphoinositide-coupled muscarinic acetylcholinoceptors activates propagating Ca2+ wave fronts in oligodendrocytes and that the characteristics of these waves are dependent on mitochondrial location and function.