# A Computational Model to Investigate GABA-Activated Astrocyte Modulation of Neuronal Excitation

Authors: Licong Li; Jin Zhou; Hongji Sun; Jing Liu; Hongrui Wang; Xiuling Liu; Changyong Wang
Journal: Computational and Mathematical Methods in Medicine
 (2020)

---

## Abstract

Gamma-aminobutyric acid (GABA) is critical for proper neural network function and can activate astrocytes to induce neuronal excitability; however, the mechanism by which astrocytes transform inhibitory signaling to excitatory enhancement remains unclear. Computational modeling can be a powerful tool to provide further understanding of how GABA-activated astrocytes modulate neuronal excitation. In the present study, we implemented a biophysical neuronal network model to investigate the effects of astrocytes on excitatory pre- and postsynaptic terminals following exposure to increasing concentrations of external GABA. The model completely describes the effects of GABA on astrocytes and excitatory presynaptic terminals within the framework of glutamatergic gliotransmission according to neurophysiological findings. Utilizing this model, our results show that astrocytes can rapidly respond to incoming GABA by inducing Ca2+ oscillations and subsequent gliotransmitter glutamate release. Elevation in GABA concentrations not only naturally decreases neuronal spikes but also enhances astrocytic glutamate release, which leads to an increase in astrocyte-mediated presynaptic release and postsynaptic slow inward currents. Neuronal excitation induced by GABA-activated astrocytes partly counteracts the inhibitory effect of GABA. Overall, the model helps to increase knowledge regarding the involvement of astrocytes in neuronal regulation using simulated bath perfusion of GABA, which may be useful for exploring the effects of GABA-type antiepileptic drugs.


---

## Body

## 1. Introduction

The brain is an adaptive nonlinear dynamic system, in which excitatory-inhibitory (EI) balance is vital for normal brain function [1]. Regulation of neuronal excitability is one of the key factors in achieving EI balance. Among the regulatory transmitters involved in neuronal excitability, GABA, released by the GABAergic neurons or astrocytes [2, 3], is the major inhibitory transmitter, which decreases neuronal excitability and prevents the neurons from overfiring. It is well established that GABA plays a central role in memory consolidation [4], motor coordination [5], and motor learning [6]. The dysregulation of GABA is closely related to certain neurological disorders, such as epilepsy [7, 8], Parkinson’s disease [9], and anxiety and depression [10]. Experimental studies and the computational models have increased our understanding of the significant functions of GABA. Many reports suggest that elevated GABA concentrations exert both inhibitory and excitatory effects on neuronal firing [11–13]. Such contradictory scenario has also been reported in bath perfusion of a GABABR agonist [14, 15]. One potential mechanism is that elevation in ambient GABA concentration can reduce the inhibitory effect of GABAergic synapses on the connected excitatory neurons, resulting in an overall increase in excitability via the disinhibition of the neurons [12, 14].Due to the tight morphological arrangement of the astrocytes and neurons, named “tripartite synapses” [16], astrocytes actively participate in the modulation of neuronal activity and synaptic behavior (for reviews, see [17–19]). Recent studies have demonstrated that GABA-activated astrocytes in the hippocampus or cortex can induce depolarized current and enhance neuronal excitability via the release of gliotransmitter glutamate [20, 21]. These findings initiate a discussion regarding astrocytes as modulators in the conversion of inhibitory stimuli to excitatory signals, which provide novel insight into the function of astrocytes in terms of neural network excitation. However, the signaling machinery involved in the intracellular cascade of GABABR-mediated Ca2+ increase and glutamate release in astrocytes remains unclear [20]. In this case, computation modeling can be used as a viable alternative to understand the role of GABA-activated astrocytes in the regulation of neuronal excitation. Most astrocyte models are based on the Ca2+ dynamics models [22, 23], which are overwhelmingly driven by synaptically released glutamate [24–26]. According to the existing glutamate-activated astrocyte mathematical framework, Li et al. developed a GABA-activated astrocyte model coupled to the seizure-firing neurons [27]. In particular, the astrocyte model added a complex differential equation and involved dual signaling of GABA released by interneuron and glutamate released by astrocyte. Nevertheless, an accurate description of the regulation of the excitatory neurons by GABA-activated astrocytes remains elusive.In the present study, we addressed this issue by devising a novel neuron-astrocyte interaction computational model and simulating injection of exogenous GABA (“virtual” GABA,GABAex=0μM,1μM,5μM,and10μM) into the extracellular space. The goal of this paper was to describe neuronal population activity and investigate the modulation of astrocytes on the neurons in stimulus conditions by lumping the microscale tripartite synapse model and macroscale neuronal-astrocytic network. On the cell level, the GABA-activated astrocyte model was incorporated into a modeling framework of glutamatergic gliotransmission [28]. The model described the activation of astrocytes by synaptic glutamate, in addition to the crosstalk between the GABABRs and mGluRs to induce astrocytic calcium elevation [29]. By considering these two signaling pathways, we may better understand the effect of GABA on astrocytes. Furthermore, synaptic regulation by astrocyte-derived glutamate is not a straightforward process in contrast to the previous synaptic models [28, 30]. At presynaptic terminals of glutamatergic synapses, the GABABR can strongly inhibit synaptic activity [31]. Thus, the comodulation of presynaptic behavior by astrocytic glutamate and increasing concentrations of exogenous GABA was considered. Using the model, we can collect neuronal spikes and network data simultaneously and understand the role of astrocytic modulation in response to exogenous GABA by recording the changes in a set of synaptic parameters over time. Numerical simulation shows that the frequency and amplitude of GABA-evoked astrocytic Ca2+ oscillations increase with increasing GABA concentrations. Moreover, the enhanced astrocytic glutamate release promotes astrocyte-regulated synaptic release and increases astrocyte-mediated postsynaptic slow inward currents (SICs, i.e., depolarizing currents); both of which can finely tune neuronal excitation. These are efficient and adjustable regulatory mechanisms for balancing the excitation and inhibition of neuronal networks.


## 2. Model and Methods

Based on neurophysiological features and experimental observations, the complete model we developed includes the inhibitory and excitatory neurons, synapses, and astrocytes, as well as the glutamate and exogenous GABA signaling pathways. The response of astrocytes to GABA occurs through activation of GABABRs. Similar to the activation of astrocytes by glutamate released from presynaptic terminals (red arrows in Figure 1), production of the second messenger inositol 1,4,5-trisphosphate (IP3) is triggered, thereby evoking Ca2+ oscillations [20]. Ca2+ elevation in astrocytes causes glutamate releasing into the extracellular space and ultimately modulating presynaptic vesicle release probability and postsynaptic neuronal excitability (orange arrows). Another important signaling pathway involves the opposing effects of astrocytic glutamate and exogenous GABA at presynaptic terminals of excitatory synapses (blue arrows), where they comodulate synaptic release. Moreover, exogenous GABA can directly regulate the probability of presynaptic GABA release and reshape the activity of postsynaptic GABAARs [13]. Given that GABA is mainly taken up by the GABAergic neurons [32], the clearance of extracellular GABA by astrocytic uptake is ignored in our model. The model architecture is presented in Figure 1, and subsequent sections describe the mathematical formulas in detail.Figure 1
The main pathways of exogenous GABA (GABAex) action in the tripartite glutamatergic synapse model. Red arrows indicate the glutamatergic gliotransmission pathways, including astrocytic metabotropic glutamate receptor (mGluR) activation, IP3-gated Ca2+ mobilization, astrocytic glutamate release, presynaptic mGluR activation, and postsynaptic AMPAR and NMDAR-mediated depolarization currents. The orange and blue arrows represent the regulation of the astrocyte and presynaptic neuron by exogenous GABA, respectively. GABAex represents simulated GABA injection to the neuronal-astrocytic model. ER refers to the endoplasmic reticulum. It is necessary to note that locations of astrocytic and neuronal receptors in this figure do not represent their exact physiological distribution.


### 2.1. Neuron Model

In the present study, the dynamics of a single neuron receiving different stimuli were simulated using a conductance-based leaky integrate-and-fire (LIF) model. This simple spiking model greatly reduces the complexity of neuronal dynamic equation and computational expense, in particular when large groups of the neurons are coupled together into neural astrocytic networks. Both the excitatory and inhibitory neurons are given by different initial values of the model variables. The neuronal model is described by the following equation:(1)τmdvdt=vrest−v+geve−v+givi−v+Iex×1gleak,where ge, gi, and gleak refer to the total excitatory, inhibitory, and leak conductance, with a reversal potential ve=0 mV and vi=−80mV, respectively. Iex is a constant input current, denoting the background current to each neuron. The neuron’s membrane potential v resets to the default value vrest=−60 mV when it reaches the set threshold vth=−50 mV. Each spike originating from presynaptic neuron i causes a conductance change △g in postsynaptic neuron j, irrespective of whether the inputs are excitatory or inhibitory [33], i.e., gj→gj+△gij. The postsynaptic conductance can be calculated as follows:(2)dgedt=−geτe,dgidt=−giτi,(3)△g=g¯j∙rt,where τe=5ms and τi=10 ms are synaptic time constants, and g¯j is the maximal postsynaptic conductance. Excitatory synaptic action was simulated by modeling AMPA and NMDA conductance, while GABAA conductance was modeled to simulate inhibitory synaptic inputs. The item rt is postsynaptic channel open probability, depending on the neurotransmitter concentration in the synaptic cleft (see Equation (22)). The values for neuronal parameters are listed in Table 1.Table 1
Parameters used in the neuronal model [34].
ParameterDescriptionValueτmMembrane time constant20 msvthFiring threshold-50 mVvrestResting potential-60 mVveExcitatory reversal potential0 mVviInhibitory reversal potential-80 mVgleakLeak conductance10 nSIexBackground current105 pAτrefAbsolute refractory period5 msτeExcitatory conductance decay time constant5 msτiInhibitory conductance decay time constant10 ms


### 2.2. Astrocyte Model

Unlike the generation of action potentials in the neurons, the computational models describing astrocytic functions always employ astrocytic Ca2+ signaling [35]. Following GABA application, GABABR-mediated Ca2+ waves depend on the intracellular IP3 cascade and are absent in IP3R2-/- mice [20]. However, a paucity of exact experimental measurement for GABA-induced IP3 activation limits the precise parameters of the IP3 dynamic model. Thus, according to experimental results of astrocyte response to GABA [20, 29] and the existing glutamate-activated astrocyte models [24–26], we simplified the model to concentrate only on the IP3-evoked Ca2+ signaling pathway, in which the rate of IP3 production depended on the amount of synaptic glutamate and exogenous GABA in the synaptic cleft. To model the dynamics of Ca2+ oscillations, we employed the classical Li–Rinzel model [23]. The modified model of astrocytic IP3 production can be explained with the following set of equations:(4)dIP3dt=IP3∗−IP3τIP3+JGABAex+Jglutamate+JG,(5)JGABAex=vgabaIP3∙GABAexn1kgaban1+GABAexn1,(6)Jglutamate=vgluIP3∙Glun2kglun2+Glun2,where the first equation describes the process of IP3 degradation to a steady state at a degradation rate τIP3=7s and an equilibrium concentration IP3∗=0.16 μM. The second and third equations describe the production rate of cytoplasmic IP3 induced by exogenous GABA and presynaptic glutamate, respectively. GABAex and Glu are the exogenous GABA and synaptic glutamate concentration, respectively. vgabaIP3 and vgluIP3 represent the rate of IP3 production via GABA and glutamate, respectively, and kgaban1 and kglun2 are their respective dissociation constants with Hill coefficients n1 and n2. Given the amplification interaction between the GABA and glutamate signaling pathways [29], we added a GABA concentration-dependent proportional coefficient k to glutamate-induced IP3 production:(7)Jglutamate=Jglutamate∙1+k∙GABAex.In our model, the available magnitude of extracellular glutamate and GABA abruptly increases following presynaptic release and exogenous input, respectively, and then exponentially decays. The estimated glutamate (Equation (20)) and GABA concentration in the synaptic cleft can be represented mathematically as follows:(8)dGABAexdt=−gGABAc∙GABAex∙Θt−Tstim,where gGABAc is the clearance rate of exogenous GABA with value ln2/3s−1 [36]. Θ is the Heaviside function. Tstim represents the constant stimulus of GABA concentration with a duration 0.5 s [27].The above dynamic IP3 model is similar to the formulation of IP3 production in [27]; however, there are some important differences. Firstly, compared with the autoreceptor-mediated astrocyte activation model, the astrocyte model we developed was coactivated by exogenous GABA and synaptically released glutamate. Also, we considered the generation of IP3 via gap junction diffusion from neighboring astrocytes in a nonlinear coupling manner [37]. JG can be described by the following two equations:(9)JG=−Fex21+tanh∆ijIP3−IP3thrωsgn∆ijIP3,(10)∆ijIP3=IP3i−IP3j,where IP3thr is the threshold of IP3 diffusion at 0.3 μM, ω is the scaling factor of diffusion, and Fex is the IP3 permeability.Finally, the dynamic equations of astrocyticIP3t are integrated into the Li–Rinzel model to describe the kinetics and properties of Ca2+ oscillations in astrocytes:(11)dCa2+dt=Jchannel−Jpump+Jleak,where the three terms denote the calcium flux from the endoplasmic reticulum (ER) to the cytosol gated by IP3 (Jchannel), the leakage flux from the ER (Jleak), and Ca2+ reuptake from the cytosol to the ER via ATP-dependent pumps (Jpump), respectively. For the sake of simplicity, the expansion of terms is omitted and the three fluxes and parameters refer to those in [38].When IP3-induced elevation in Ca2+ concentration exceeds the threshold Cθ at time tj, a fraction rAtj of astrocytic glutamate resources can be released into the extracellular space by vesicle exocytosis [39]. The fraction of readily releasable is given as follows:(12)rAt=UA∙xAt,where UA is the resting glutamate release probability and xA is the fraction of available glutamate resources for release, according to the following:(13)dxAdt=1−xAτG−rAt∙δt−tj,where τG is the glutamate resource reintegration time constant. The astrocyte-derived glutamate concentration GA in the extrasynaptic cleft is described as follows:(14)dGAdt=−gAc∙GA+rA∙ϱe∙GT∙δt−tj∙1+k∙GABAex,where gAc is the clearance rate of glutamate, ϱe denotes the volume ratio of vesicles to periastrocytic space, and GT is the total vesicular glutamate concentration in astrocytes. The parameters are given in Table 2.Table 2
Parameters used in the astrocyte model [24, 28, 37].
ParameterDescriptionValueIP3∗Steady-state concentration of IP30.16μMτIP3IP3 degradation constant7 svgabaIP3Rate of IP3 production via GABA0.1625μM s-1vgluIP3Rate of IP3 production via glutamate0.062μM s-1n1GABA Hill coefficient0.3n2Glutamate Hill coefficient0.2kGABAnDissociation constant for GABA-stimulated IP3 production0.6μMkglunDissociation constant for glutamate-stimulated IP3 production0.78μMIP3thrThreshold for IP3 diffusion0.3μMFexIP3 permeability0.09μM s-1ωScaling factor of IP3 diffusion0.05μMCθCa2+ threshold for exocytosis0.2μMτGGlutamate reintegration time constant1.66 sUAResting glutamate release probability0.6ϱeVolume ratio of vesicles to periastrocytic space6.5×10−4gAcClearance rate of astrocytic glutamate60 s-1GTTotal vesicular glutamate concentration200 mMkGABA concentration-dependent proportionality coefficient0.3


### 2.3. Synapse Model

The effects of astrocytic glutamate and exogenous GABA on presynaptic terminals are mediated by activation of mGluRs and GABABRs, respectively. Both of these receptor types regulate the amount of Ca2+ influx into presynaptic terminals; the former promotes [40] and the latter retards [41]. Release of vesicles is dependent on elevation in intracellular Ca2+ concentration [42]. Therefore, open proportion U of Ca2+ channels modulated by the two transmitters is a critical parameter that decides the probability of presynaptic neurotransmitter release. To mimic the effects of astrocytic glutamate and exogenous GABA comodulation on synaptic release probability, we used a modified kinetic model of synaptic release regulated by astrocyte [28], which is based on the Tsodyks–Markram (TM) phenomenological model of synaptic activity [43], and added GABA activation kinetic scheme. The corresponding kinetic equation reads as follows:(15)U=U0+ξ−U0∙Γ−U0∙rGABAB,where U0 is defined as a constant in the TM model [44], denoting the open ratio of Ca2+ channels induced by action potentials, i.e., resting synaptic release probability. The second term is the effect of astrocyte-released glutamate on presynaptic glutamate receptors, in which Γ is the proportion of activated receptors and parameter ξ0≤ξ≤1 determines the type and strength of the action of astrocytes on presynaptic terminals. The third term denotes the probability of modulation of synaptic basal release by GABA, and rGABAB is the fraction of GABAB-mediated Ca2+ channels that close. The probability that channels are open or closed is determined by the following:(16)dΓdt=OG∙GA∙1−Γ−ΩG∙Γ,(17)drGABABdt=αGABAB∙GABAex∙1−rGABAB−βGABAB∙rGABAB,where GABAex and GA represent the concentrations of exogenous GABA and astrocyte-released glutamate, respectively, given by Equation (8) and Equation (14). The parameters OG and ΩG are the forward and backward rate constants for astrocytic glutamate binding to excitatory presynaptic receptors, respectively, and αGABAB and βGABAB are the rate constants of GABABR-mediated binding or unbinding, respectively.Accordingly, the modified resting synaptic release probabilityU will determine the fraction uS of available neurotransmitter to be utilized occurring at each presynaptic spike time tk:(18)duSdt=−uSτfac+U∙1−uS∙δt−tk,where τfac is the relaxation time constant of facilitation and determines the rate of uS decaying to zero in the absence of spikes. The fraction of total neurotransmitter available for release xS is described by the following:(19)dxSdt=1−xSτrec−rS∙δt−tk,where τrec is the recovery time constant and determines the reintegration rate of resources to the available pool. rS denotes the fraction of neurotransmitter effectively released into the extracellular space, i.e., probability of neurotransmitter release (Pr), defined as rS=uS∙xS. The estimated glutamate concentration GS from presynaptic release can be represented mathematically as follows:(20)dGSdt=−gSc∙GS+rS∙ϱc∙YT∙δt−tk,where gSc is the clearance rate of glutamate via neuronal and astrocytic uptake or spillover (flowing out of the synaptic cleft). ϱc is the vesicular to mixing volume ratio, and YT represents the total vesicular glutamate concentration in the presynaptic bouton. Akin to those described in glutamate release process in Equations (18)–(20), when a single action potential arrived at GABAergic synapse, GABA in synaptic vesicles are released from presynaptic terminal and then diffuse in the synaptic cleft, where they are retrieved by endocytosis [45]. Because this study focused on the regulation of astrocytes on the excitatory neurons, the open proportion U in presynaptic terminal of GABAergic synapse was set a constant for simplification.Once vesicles are released from presynaptic terminals or astrocytes, glutamate diffuses across the synaptic cleft and binds to postsynaptic AMPA and NMDA receptors, resulting in an increase in conductance. The dynamics of the synaptic conductance gating variablert in Equation (3) can be represented using a two-state model [46]. Transition between closed and open states of the channel is considered a Markov process based on the following diagram:(21)C+Tα⇌βO,where C and O indicate the unbound and bound states of the channel to neurotransmitter, respectively, T is the glutamate or GABA concentration, and α and β are the forward and backward rate constants, respectively. The two-state kinetic model describes the process of direct receptor binding to neurotransmitter, changing the ion channel. The fraction of the receptors in the open state r is described by the first-order dynamic equation:(22)drdt=α∙T∙1−r−β∙r.Finally, NMDA-, AMPA-, and GABAA-mediated postsynaptic currents are expressed as functions of rt and the membrane potential v by the following equations [46]:(23)IAMPA=g¯AMPA∙rt∙v−EAMPA,(24)INMDA=g¯NMDA∙Mgv∙rt∙v−ENMDA,(25)IGABAa=g¯GABAA∙rt∙v−EGABAA,where g¯ is the maximum conductance, g¯AMPA=0.35 nS, g¯NMDA=0.1 nS, g¯GABAA=0.62 nS, v is the postsynaptic membrane potential, E is the synaptic reversal potential, EAMPA=ENMDA=0 mV, and EGABAA=−80 mV. Notably, a unique property of NMDA currents is that NMDA receptors contain a voltage-dependent term representing magnesium (Mg2+) block at Mg2+=1mM [47]:(26)Mgv=11+exp−0.062∗vMg2+/3.57.The slow inward currentIsic elicited by astrocytic glutamate can be defined as follows [48, 49]:(27)Isict=INMDAt+IAMPAt.In this fashion, postsynaptic mechanisms of astrocytic glutamate regulation mediated by NMDA and AMPA receptors complicate the characterization of regular postsynaptic currents by glutamate released from presynaptic terminals. Thus, the two types of excitatory events (SICs and EPSCs) can be distinguished according to their origin, astrocyte-derived glutamateGA (Equation (14)) and synaptically released glutamate GS (Equation (20)), respectively. Noted worthily, although the NMDA receptors that mediate SICs and EPSCs are different [50], NMDA receptor subtypes are often not distinguished in the computational models for simplicity. The parameters for the synapse model are listed in Table 3.Table 3
Parameters used in the synapse model [28, 46].
ParameterDescriptionValueU0Resting synaptic release probability0.3ξType and strength of astrocytes acting on presynapses0.8OGActivation rate of astrocytic glutamate1.5 M-1 s-1ΩGInactivation rate of astrocytic glutamate0.2 s-1αGABABGABAB forward rate constant16μM-1 s-1βGABABGABAB backward rate constant6 s-1τfacFacilitation time constant0.3 sτrecRecovery time constant0.5 sgScPresynaptic glutamate clearance rate40 s-1ϱcVesicular to mixed volume ratio0.005YTTotal vesicular glutamate concentration500 mMαAMPAAMPA forward rate constant1.1μM-1 s-1βAMPAAMPA backward rate constant190 s-1αNMDANMDA forward rate constant0.072μM-1 s-1βNMDANMDA backward rate constant6.6 s-1αGABAAGABAA forward rate constant0.53μM-1 s-1βGABAAGABAA backward rate constant180 s-1gGABAuptakeGABA reuptake rate constant6 s-1


### 2.4. Neuronal–Astrocyte Network

Based on the network components mentioned above, we simulated a two-dimensional network model (Figure2). Each neuron and astrocyte were assigned a spatial location in a square grid with dimensions of 10×10mm planar units, providing a suitable physical arrangement for signaling transmission among them [51]. The neuronal network model comprises 400 excitatory and 100 inhibitory neurons with 0.2 connection probability as well as 400 astrocytes. Experimental results show that individual cortical astrocyte can enwrap the nearby neurons within its territory and modulate the synapses associated with hundreds of dendrites [52]. According to their spatial coordinates, therefore, we calculated the minimum distance between astrocytes and the pre- and postsynaptic neurons and then determined which astrocyte should be responsible for this synapse. Connection between astrocytes was also based on their coordinates. Considering a mean numerical density of astrocytes of 17,575 cells per mm3 in the rat visual cortex [53], average intersomal distance between astrocytes was set 675 μm in our model by enlarging in proportion. In this way, an astrocyte can be connected to about 100 nearby excitatory synapses and 4 astrocytes on average. Multiple variables were calculated and recorded during each 10 s simulation, including neuronal firing, astrocytic calcium oscillations, and excitatory presynaptic release. The model was implemented in the Brian 2.0 simulator [54, 55].Figure 2
The 2D neural astrocytic network model. Astrocytes, the inhibitory neurons, and the excitatory neurons are represented by green stars and blue and red dots, respectively. Green lines represent astrocyte network and the neuronal connections are omitted for clarity.


## 2.1. Neuron Model

In the present study, the dynamics of a single neuron receiving different stimuli were simulated using a conductance-based leaky integrate-and-fire (LIF) model. This simple spiking model greatly reduces the complexity of neuronal dynamic equation and computational expense, in particular when large groups of the neurons are coupled together into neural astrocytic networks. Both the excitatory and inhibitory neurons are given by different initial values of the model variables. The neuronal model is described by the following equation:(1)τmdvdt=vrest−v+geve−v+givi−v+Iex×1gleak,where ge, gi, and gleak refer to the total excitatory, inhibitory, and leak conductance, with a reversal potential ve=0 mV and vi=−80mV, respectively. Iex is a constant input current, denoting the background current to each neuron. The neuron’s membrane potential v resets to the default value vrest=−60 mV when it reaches the set threshold vth=−50 mV. Each spike originating from presynaptic neuron i causes a conductance change △g in postsynaptic neuron j, irrespective of whether the inputs are excitatory or inhibitory [33], i.e., gj→gj+△gij. The postsynaptic conductance can be calculated as follows:(2)dgedt=−geτe,dgidt=−giτi,(3)△g=g¯j∙rt,where τe=5ms and τi=10 ms are synaptic time constants, and g¯j is the maximal postsynaptic conductance. Excitatory synaptic action was simulated by modeling AMPA and NMDA conductance, while GABAA conductance was modeled to simulate inhibitory synaptic inputs. The item rt is postsynaptic channel open probability, depending on the neurotransmitter concentration in the synaptic cleft (see Equation (22)). The values for neuronal parameters are listed in Table 1.Table 1
Parameters used in the neuronal model [34].
ParameterDescriptionValueτmMembrane time constant20 msvthFiring threshold-50 mVvrestResting potential-60 mVveExcitatory reversal potential0 mVviInhibitory reversal potential-80 mVgleakLeak conductance10 nSIexBackground current105 pAτrefAbsolute refractory period5 msτeExcitatory conductance decay time constant5 msτiInhibitory conductance decay time constant10 ms


## 2.2. Astrocyte Model

Unlike the generation of action potentials in the neurons, the computational models describing astrocytic functions always employ astrocytic Ca2+ signaling [35]. Following GABA application, GABABR-mediated Ca2+ waves depend on the intracellular IP3 cascade and are absent in IP3R2-/- mice [20]. However, a paucity of exact experimental measurement for GABA-induced IP3 activation limits the precise parameters of the IP3 dynamic model. Thus, according to experimental results of astrocyte response to GABA [20, 29] and the existing glutamate-activated astrocyte models [24–26], we simplified the model to concentrate only on the IP3-evoked Ca2+ signaling pathway, in which the rate of IP3 production depended on the amount of synaptic glutamate and exogenous GABA in the synaptic cleft. To model the dynamics of Ca2+ oscillations, we employed the classical Li–Rinzel model [23]. The modified model of astrocytic IP3 production can be explained with the following set of equations:(4)dIP3dt=IP3∗−IP3τIP3+JGABAex+Jglutamate+JG,(5)JGABAex=vgabaIP3∙GABAexn1kgaban1+GABAexn1,(6)Jglutamate=vgluIP3∙Glun2kglun2+Glun2,where the first equation describes the process of IP3 degradation to a steady state at a degradation rate τIP3=7s and an equilibrium concentration IP3∗=0.16 μM. The second and third equations describe the production rate of cytoplasmic IP3 induced by exogenous GABA and presynaptic glutamate, respectively. GABAex and Glu are the exogenous GABA and synaptic glutamate concentration, respectively. vgabaIP3 and vgluIP3 represent the rate of IP3 production via GABA and glutamate, respectively, and kgaban1 and kglun2 are their respective dissociation constants with Hill coefficients n1 and n2. Given the amplification interaction between the GABA and glutamate signaling pathways [29], we added a GABA concentration-dependent proportional coefficient k to glutamate-induced IP3 production:(7)Jglutamate=Jglutamate∙1+k∙GABAex.In our model, the available magnitude of extracellular glutamate and GABA abruptly increases following presynaptic release and exogenous input, respectively, and then exponentially decays. The estimated glutamate (Equation (20)) and GABA concentration in the synaptic cleft can be represented mathematically as follows:(8)dGABAexdt=−gGABAc∙GABAex∙Θt−Tstim,where gGABAc is the clearance rate of exogenous GABA with value ln2/3s−1 [36]. Θ is the Heaviside function. Tstim represents the constant stimulus of GABA concentration with a duration 0.5 s [27].The above dynamic IP3 model is similar to the formulation of IP3 production in [27]; however, there are some important differences. Firstly, compared with the autoreceptor-mediated astrocyte activation model, the astrocyte model we developed was coactivated by exogenous GABA and synaptically released glutamate. Also, we considered the generation of IP3 via gap junction diffusion from neighboring astrocytes in a nonlinear coupling manner [37]. JG can be described by the following two equations:(9)JG=−Fex21+tanh∆ijIP3−IP3thrωsgn∆ijIP3,(10)∆ijIP3=IP3i−IP3j,where IP3thr is the threshold of IP3 diffusion at 0.3 μM, ω is the scaling factor of diffusion, and Fex is the IP3 permeability.Finally, the dynamic equations of astrocyticIP3t are integrated into the Li–Rinzel model to describe the kinetics and properties of Ca2+ oscillations in astrocytes:(11)dCa2+dt=Jchannel−Jpump+Jleak,where the three terms denote the calcium flux from the endoplasmic reticulum (ER) to the cytosol gated by IP3 (Jchannel), the leakage flux from the ER (Jleak), and Ca2+ reuptake from the cytosol to the ER via ATP-dependent pumps (Jpump), respectively. For the sake of simplicity, the expansion of terms is omitted and the three fluxes and parameters refer to those in [38].When IP3-induced elevation in Ca2+ concentration exceeds the threshold Cθ at time tj, a fraction rAtj of astrocytic glutamate resources can be released into the extracellular space by vesicle exocytosis [39]. The fraction of readily releasable is given as follows:(12)rAt=UA∙xAt,where UA is the resting glutamate release probability and xA is the fraction of available glutamate resources for release, according to the following:(13)dxAdt=1−xAτG−rAt∙δt−tj,where τG is the glutamate resource reintegration time constant. The astrocyte-derived glutamate concentration GA in the extrasynaptic cleft is described as follows:(14)dGAdt=−gAc∙GA+rA∙ϱe∙GT∙δt−tj∙1+k∙GABAex,where gAc is the clearance rate of glutamate, ϱe denotes the volume ratio of vesicles to periastrocytic space, and GT is the total vesicular glutamate concentration in astrocytes. The parameters are given in Table 2.Table 2
Parameters used in the astrocyte model [24, 28, 37].
ParameterDescriptionValueIP3∗Steady-state concentration of IP30.16μMτIP3IP3 degradation constant7 svgabaIP3Rate of IP3 production via GABA0.1625μM s-1vgluIP3Rate of IP3 production via glutamate0.062μM s-1n1GABA Hill coefficient0.3n2Glutamate Hill coefficient0.2kGABAnDissociation constant for GABA-stimulated IP3 production0.6μMkglunDissociation constant for glutamate-stimulated IP3 production0.78μMIP3thrThreshold for IP3 diffusion0.3μMFexIP3 permeability0.09μM s-1ωScaling factor of IP3 diffusion0.05μMCθCa2+ threshold for exocytosis0.2μMτGGlutamate reintegration time constant1.66 sUAResting glutamate release probability0.6ϱeVolume ratio of vesicles to periastrocytic space6.5×10−4gAcClearance rate of astrocytic glutamate60 s-1GTTotal vesicular glutamate concentration200 mMkGABA concentration-dependent proportionality coefficient0.3


## 2.3. Synapse Model

The effects of astrocytic glutamate and exogenous GABA on presynaptic terminals are mediated by activation of mGluRs and GABABRs, respectively. Both of these receptor types regulate the amount of Ca2+ influx into presynaptic terminals; the former promotes [40] and the latter retards [41]. Release of vesicles is dependent on elevation in intracellular Ca2+ concentration [42]. Therefore, open proportion U of Ca2+ channels modulated by the two transmitters is a critical parameter that decides the probability of presynaptic neurotransmitter release. To mimic the effects of astrocytic glutamate and exogenous GABA comodulation on synaptic release probability, we used a modified kinetic model of synaptic release regulated by astrocyte [28], which is based on the Tsodyks–Markram (TM) phenomenological model of synaptic activity [43], and added GABA activation kinetic scheme. The corresponding kinetic equation reads as follows:(15)U=U0+ξ−U0∙Γ−U0∙rGABAB,where U0 is defined as a constant in the TM model [44], denoting the open ratio of Ca2+ channels induced by action potentials, i.e., resting synaptic release probability. The second term is the effect of astrocyte-released glutamate on presynaptic glutamate receptors, in which Γ is the proportion of activated receptors and parameter ξ0≤ξ≤1 determines the type and strength of the action of astrocytes on presynaptic terminals. The third term denotes the probability of modulation of synaptic basal release by GABA, and rGABAB is the fraction of GABAB-mediated Ca2+ channels that close. The probability that channels are open or closed is determined by the following:(16)dΓdt=OG∙GA∙1−Γ−ΩG∙Γ,(17)drGABABdt=αGABAB∙GABAex∙1−rGABAB−βGABAB∙rGABAB,where GABAex and GA represent the concentrations of exogenous GABA and astrocyte-released glutamate, respectively, given by Equation (8) and Equation (14). The parameters OG and ΩG are the forward and backward rate constants for astrocytic glutamate binding to excitatory presynaptic receptors, respectively, and αGABAB and βGABAB are the rate constants of GABABR-mediated binding or unbinding, respectively.Accordingly, the modified resting synaptic release probabilityU will determine the fraction uS of available neurotransmitter to be utilized occurring at each presynaptic spike time tk:(18)duSdt=−uSτfac+U∙1−uS∙δt−tk,where τfac is the relaxation time constant of facilitation and determines the rate of uS decaying to zero in the absence of spikes. The fraction of total neurotransmitter available for release xS is described by the following:(19)dxSdt=1−xSτrec−rS∙δt−tk,where τrec is the recovery time constant and determines the reintegration rate of resources to the available pool. rS denotes the fraction of neurotransmitter effectively released into the extracellular space, i.e., probability of neurotransmitter release (Pr), defined as rS=uS∙xS. The estimated glutamate concentration GS from presynaptic release can be represented mathematically as follows:(20)dGSdt=−gSc∙GS+rS∙ϱc∙YT∙δt−tk,where gSc is the clearance rate of glutamate via neuronal and astrocytic uptake or spillover (flowing out of the synaptic cleft). ϱc is the vesicular to mixing volume ratio, and YT represents the total vesicular glutamate concentration in the presynaptic bouton. Akin to those described in glutamate release process in Equations (18)–(20), when a single action potential arrived at GABAergic synapse, GABA in synaptic vesicles are released from presynaptic terminal and then diffuse in the synaptic cleft, where they are retrieved by endocytosis [45]. Because this study focused on the regulation of astrocytes on the excitatory neurons, the open proportion U in presynaptic terminal of GABAergic synapse was set a constant for simplification.Once vesicles are released from presynaptic terminals or astrocytes, glutamate diffuses across the synaptic cleft and binds to postsynaptic AMPA and NMDA receptors, resulting in an increase in conductance. The dynamics of the synaptic conductance gating variablert in Equation (3) can be represented using a two-state model [46]. Transition between closed and open states of the channel is considered a Markov process based on the following diagram:(21)C+Tα⇌βO,where C and O indicate the unbound and bound states of the channel to neurotransmitter, respectively, T is the glutamate or GABA concentration, and α and β are the forward and backward rate constants, respectively. The two-state kinetic model describes the process of direct receptor binding to neurotransmitter, changing the ion channel. The fraction of the receptors in the open state r is described by the first-order dynamic equation:(22)drdt=α∙T∙1−r−β∙r.Finally, NMDA-, AMPA-, and GABAA-mediated postsynaptic currents are expressed as functions of rt and the membrane potential v by the following equations [46]:(23)IAMPA=g¯AMPA∙rt∙v−EAMPA,(24)INMDA=g¯NMDA∙Mgv∙rt∙v−ENMDA,(25)IGABAa=g¯GABAA∙rt∙v−EGABAA,where g¯ is the maximum conductance, g¯AMPA=0.35 nS, g¯NMDA=0.1 nS, g¯GABAA=0.62 nS, v is the postsynaptic membrane potential, E is the synaptic reversal potential, EAMPA=ENMDA=0 mV, and EGABAA=−80 mV. Notably, a unique property of NMDA currents is that NMDA receptors contain a voltage-dependent term representing magnesium (Mg2+) block at Mg2+=1mM [47]:(26)Mgv=11+exp−0.062∗vMg2+/3.57.The slow inward currentIsic elicited by astrocytic glutamate can be defined as follows [48, 49]:(27)Isict=INMDAt+IAMPAt.In this fashion, postsynaptic mechanisms of astrocytic glutamate regulation mediated by NMDA and AMPA receptors complicate the characterization of regular postsynaptic currents by glutamate released from presynaptic terminals. Thus, the two types of excitatory events (SICs and EPSCs) can be distinguished according to their origin, astrocyte-derived glutamateGA (Equation (14)) and synaptically released glutamate GS (Equation (20)), respectively. Noted worthily, although the NMDA receptors that mediate SICs and EPSCs are different [50], NMDA receptor subtypes are often not distinguished in the computational models for simplicity. The parameters for the synapse model are listed in Table 3.Table 3
Parameters used in the synapse model [28, 46].
ParameterDescriptionValueU0Resting synaptic release probability0.3ξType and strength of astrocytes acting on presynapses0.8OGActivation rate of astrocytic glutamate1.5 M-1 s-1ΩGInactivation rate of astrocytic glutamate0.2 s-1αGABABGABAB forward rate constant16μM-1 s-1βGABABGABAB backward rate constant6 s-1τfacFacilitation time constant0.3 sτrecRecovery time constant0.5 sgScPresynaptic glutamate clearance rate40 s-1ϱcVesicular to mixed volume ratio0.005YTTotal vesicular glutamate concentration500 mMαAMPAAMPA forward rate constant1.1μM-1 s-1βAMPAAMPA backward rate constant190 s-1αNMDANMDA forward rate constant0.072μM-1 s-1βNMDANMDA backward rate constant6.6 s-1αGABAAGABAA forward rate constant0.53μM-1 s-1βGABAAGABAA backward rate constant180 s-1gGABAuptakeGABA reuptake rate constant6 s-1


## 2.4. Neuronal–Astrocyte Network

Based on the network components mentioned above, we simulated a two-dimensional network model (Figure2). Each neuron and astrocyte were assigned a spatial location in a square grid with dimensions of 10×10mm planar units, providing a suitable physical arrangement for signaling transmission among them [51]. The neuronal network model comprises 400 excitatory and 100 inhibitory neurons with 0.2 connection probability as well as 400 astrocytes. Experimental results show that individual cortical astrocyte can enwrap the nearby neurons within its territory and modulate the synapses associated with hundreds of dendrites [52]. According to their spatial coordinates, therefore, we calculated the minimum distance between astrocytes and the pre- and postsynaptic neurons and then determined which astrocyte should be responsible for this synapse. Connection between astrocytes was also based on their coordinates. Considering a mean numerical density of astrocytes of 17,575 cells per mm3 in the rat visual cortex [53], average intersomal distance between astrocytes was set 675 μm in our model by enlarging in proportion. In this way, an astrocyte can be connected to about 100 nearby excitatory synapses and 4 astrocytes on average. Multiple variables were calculated and recorded during each 10 s simulation, including neuronal firing, astrocytic calcium oscillations, and excitatory presynaptic release. The model was implemented in the Brian 2.0 simulator [54, 55].Figure 2
The 2D neural astrocytic network model. Astrocytes, the inhibitory neurons, and the excitatory neurons are represented by green stars and blue and red dots, respectively. Green lines represent astrocyte network and the neuronal connections are omitted for clarity.


## 3. Results

### 3.1. Astrocytic Glutamate Increases the Average Network Firing Rate

We began our analysis by describing the firing behavior of the neurons. The neural-astrocytic network was treated with increasing GABA concentrations (0, 1, 5, and 10μM). A quarter of the neurons were chosen for clarification using raster plot (Figure 3). The network firing rate, measured as an average across the spikes of all neurons (time bin width ∆t=1ms), is shown in the bottom panel of the raster plot of neuronal firing. Furthermore, we fixed inhibitory synaptic release probability at a constant of 0.5 [56] to ignore the effect of GABABR-mediated disinhibition at inhibitory presynaptic terminals. In order to produce weak and strong inhibitory effects on excitatory inputs, maximal conductance of GABAA increases with the change of exogenous concentration depending on the experimental results [13]. Figures 3(a)–3(d) (top panel) show that the firing activity of the excitatory (red) and inhibitory (blue) neurons during treatment with different concentrations of GABA decreased, in line with the previous experimental observation that we performed acute GABA treatment on the cortical cultures [57]. The release time of gliotransmitter glutamate (green) from astrocytes was significantly affected. Astrocytic glutamate triggered synchronized responses in the neurons, which increased the average network firing rate (bottom panel). The synchronized neuronal activities are consistent with the results of electrophysiological experiments [50, 58] and computational simulation [59], which imply that astrocytes have a vital role in modulation the firing of the neurons.Figure 3
The raster plot and mean firing rate of neural activity during treatment with different concentrations of GABA: (a) 0μM, (b) 1 μM, (c) 5 μM, and (d) 10 μM. The top panel shows astrocytic glutamate release (green) and excitatory (red) and inhibitory (blue) neuronal firing. The average network firing rate is shown at the bottom.
(a)(b)(c)(d)


### 3.2. GABA-Activated Astrocytes Augment Calcium Oscillations and Glutamate Release

The key model parameters for GABA-induced astrocytic activity are shown in Figure4, including the intracellular cascade of IP3-mediated Ca2+ response and Ca2+-dependent glutamate release. The dynamic changes in intracellular Ca2+ in astrocyte were regulated by the integration of the exogenous GABA and presynaptic glutamate signaling pathways (Figure 1), akin to spatial and temporal integrators [60]. With the exogenous GABA concentration increased (GABAex=0,1,5,and10μM), the concentration of glutamate resulting from presynaptic release decreased significantly (Figure 4(a)), while the activation of glutamate receptors on astrocyte increased (Figure 4(b)). This may be attributed to the amplification of mGluRs via crosstalk with the GABABRs [29]. In these cases, compared to the condition without exogenous stimuli (Figure 4(c), black line), the time course of IP3 production differed across the three paradigms due to the increase of GABA concentrations (GABAex=1,5,and10μM) in the synaptic cleft (Figure 4(c), colored lines), which resulted in greater Ca2+ oscillatory behavior (Figure 4(d)). In the Ca2+ variation curve, the single pulse elevation may arise from direct stimulation with exogenous GABA, and the subsequent long-lasting Ca2+ oscillations could be evoked by glutamate [27]. Indeed, the astrocytic Ca2+ response closely depends on the amount of glutamate released from presynaptic terminals in the synaptic cleft [51]. Although the neurons were suppressed, an amplification interaction between the GABA and glutamate signaling pathways augmented the amplitude and frequency of Ca2+ oscillations in astrocyte. Simulated results indicate that a higher concentration of GABA leads to a stronger and prolonged Ca2+ wave, which changes the rate of glutamate release (Figure 4(e)). Nevertheless, there are no experimental data available to precisely quantitate the effect of these two receptors. Further experimental studies are required to elucidate the exact mechanism and quantitative properties of Ca2+ oscillations evoked by multiple neurotransmitters.Figure 4
The dynamic response of astrocyte to presynaptic glutamate and exogenous GABA. (a) The amount of presynaptically released glutamate (Glusyn) decreases as the concentration of exogenous GABA increases, but (b) there is greater activation of mGluRs on astrocyte due to the crosstalk with GABAB receptors, which results in (c) IP3 production and (d) IP3-gated [Ca2+] elevation. (e) Once the [Ca2+] threshold has been reached (0.2 μM, red dashed line), astrocyte releases glutamate (Gluastro) into the extracellular space. To avoid triggering repeatedly the release in all the time steps, gliotransmitter release event is triggered only once in threshold conditions, a mechanism similar to the refractory period of the neuron.


### 3.3. Enhanced Astrocytic Glutamate Promotes Presynaptic Release

The variations in presynaptic release probability in our model were regulated by astrocytic glutamate and exogenous GABA (Figure1). As a contrast, the direct response of presynaptic terminal to exogenous GABA in the absence of astrocytic modulation is shown in Figure 5(a). The release probability Pr significantly decreased with GABA inhibition, while Pr in central excitatory synapses was between 0.1 and 0.9 [61]. This indicates that increased GABA concentrations robustly inhibit excitatory presynaptic release [31]. In the presence of astrocyte-derived glutamate, the scenario could be reversed (Figure 5(b)). In particular, the regulation of presynaptic release by astrocyte-derived glutamate was significantly enhanced as the concentration of exogenous GABA increased. The rate of astrocytic glutamate potentiation was larger than the rate of GABA depression, which gave rise to an increase in Pr, illustrating that glutamate exocytosis from astrocytes promoted presynaptic transmitter release [40]. Furthermore, the variations in excitatory postsynaptic currents (EPSCs) elicited by presynaptic glutamate release are consistent with the quantal size variability in presynaptic vesicle release (Figure 6). The results suggest that exogenous GABA causes a large decrease in the EPSC amplitude [14], whereas the release of astrocytic glutamate enhances excitatory postsynaptic currents.Figure 5
Variations in presynaptic glutamate release probability (Pr) during (a) astrocyte-independent and (b) astrocyte-dependent pathways. The dots represent each presynaptic release event.Figure 6
Variations in excitatory postsynaptic currents (EPSCs) induced by glutamate released from presynaptic terminal during treatment with different concentrations of GABA: (a) 0μM, (b) 1 μM, (c) 5 μM, and (d) 10 μM.


### 3.4. Astrocyte-Mediated Slow Inward Currents Increase Postsynaptic Neuronal Excitability

Both excitatory and inhibitory inputs to the postsynaptic neurons are accompanied by a decrease in impedance and an increase in membrane conductance [33]. The inhibitory and excitatory synaptic conductance values gi and ge were consistently balanced at a certain ratio [62]; the former was several fold higher than the latter [63]. In our model, postsynaptic events were characterized by changes in current instead of changes in conductance. To investigate the impact of NMDAR- and AMPAR-mediated SICs on the postsynaptic neurons, we removed the excitatory effect of regular synaptic currents driven by synaptically released glutamate. These results could be the direct response of the postsynaptic neuron to astrocyte.GABA-evoked astrocytic Ca2+ oscillations triggered the release of glutamate (Figure 7(a)), which diffused into the extracellular space and bound to post-extrasynaptic glutamate receptors, eliciting depolarizing SICs (Figures 7(b)–7(e)). Astrocyte-mediated SICs caused membrane depolarization in the postsynaptic neurons, which was sufficient to reach the firing threshold and ultimately influenced the time course of postsynaptic neuronal firing (Figure 7). Pathological conditions, such as neuronal depolarization, termed paroxysmal depolarization shifts, may contribute to seizure generation [64]. Indeed, intracellular Ca2+ elevation in astrocytes is a crucial factor in the regulation of SICs, which has been previously reported in in vitro [49] and numerical studies [48, 65]. Moreover, the SIC amplitude is determined by the extracellular glutamate concentration. Both an enhanced astrocytic glutamate release and a downregulated glutamate clearance rate can dramatically modulate neuronal excitability, therefore increasing seizure susceptibility [48, 66, 67]. The results of this simulation show that the higher concentration of exogenous GABA is injected to neuronal-astrocytic network, the larger amplitude of astrocyte-mediated SICs is evoked due to the enhanced release of astrocytic glutamate. The increase in excitability could partly counteract the resulting decrease in EPSCs induced by presynaptic glutamate, thereby fine-tuning neuronal network excitation.Figure 7
Astrocyte-mediated SICs and the change in postsynaptic membrane potential. (a) Enhanced glutamate release by astrocytes due to the elevated concentrations of GABA ((b–e)GABAex=0μM,1μM,5μM,and10μM) results in a stronger Isic and subsequently makes the neurons more excitable to resist the increased GABA inhibition.


## 3.1. Astrocytic Glutamate Increases the Average Network Firing Rate

We began our analysis by describing the firing behavior of the neurons. The neural-astrocytic network was treated with increasing GABA concentrations (0, 1, 5, and 10μM). A quarter of the neurons were chosen for clarification using raster plot (Figure 3). The network firing rate, measured as an average across the spikes of all neurons (time bin width ∆t=1ms), is shown in the bottom panel of the raster plot of neuronal firing. Furthermore, we fixed inhibitory synaptic release probability at a constant of 0.5 [56] to ignore the effect of GABABR-mediated disinhibition at inhibitory presynaptic terminals. In order to produce weak and strong inhibitory effects on excitatory inputs, maximal conductance of GABAA increases with the change of exogenous concentration depending on the experimental results [13]. Figures 3(a)–3(d) (top panel) show that the firing activity of the excitatory (red) and inhibitory (blue) neurons during treatment with different concentrations of GABA decreased, in line with the previous experimental observation that we performed acute GABA treatment on the cortical cultures [57]. The release time of gliotransmitter glutamate (green) from astrocytes was significantly affected. Astrocytic glutamate triggered synchronized responses in the neurons, which increased the average network firing rate (bottom panel). The synchronized neuronal activities are consistent with the results of electrophysiological experiments [50, 58] and computational simulation [59], which imply that astrocytes have a vital role in modulation the firing of the neurons.Figure 3
The raster plot and mean firing rate of neural activity during treatment with different concentrations of GABA: (a) 0μM, (b) 1 μM, (c) 5 μM, and (d) 10 μM. The top panel shows astrocytic glutamate release (green) and excitatory (red) and inhibitory (blue) neuronal firing. The average network firing rate is shown at the bottom.
(a)(b)(c)(d)


## 3.2. GABA-Activated Astrocytes Augment Calcium Oscillations and Glutamate Release

The key model parameters for GABA-induced astrocytic activity are shown in Figure4, including the intracellular cascade of IP3-mediated Ca2+ response and Ca2+-dependent glutamate release. The dynamic changes in intracellular Ca2+ in astrocyte were regulated by the integration of the exogenous GABA and presynaptic glutamate signaling pathways (Figure 1), akin to spatial and temporal integrators [60]. With the exogenous GABA concentration increased (GABAex=0,1,5,and10μM), the concentration of glutamate resulting from presynaptic release decreased significantly (Figure 4(a)), while the activation of glutamate receptors on astrocyte increased (Figure 4(b)). This may be attributed to the amplification of mGluRs via crosstalk with the GABABRs [29]. In these cases, compared to the condition without exogenous stimuli (Figure 4(c), black line), the time course of IP3 production differed across the three paradigms due to the increase of GABA concentrations (GABAex=1,5,and10μM) in the synaptic cleft (Figure 4(c), colored lines), which resulted in greater Ca2+ oscillatory behavior (Figure 4(d)). In the Ca2+ variation curve, the single pulse elevation may arise from direct stimulation with exogenous GABA, and the subsequent long-lasting Ca2+ oscillations could be evoked by glutamate [27]. Indeed, the astrocytic Ca2+ response closely depends on the amount of glutamate released from presynaptic terminals in the synaptic cleft [51]. Although the neurons were suppressed, an amplification interaction between the GABA and glutamate signaling pathways augmented the amplitude and frequency of Ca2+ oscillations in astrocyte. Simulated results indicate that a higher concentration of GABA leads to a stronger and prolonged Ca2+ wave, which changes the rate of glutamate release (Figure 4(e)). Nevertheless, there are no experimental data available to precisely quantitate the effect of these two receptors. Further experimental studies are required to elucidate the exact mechanism and quantitative properties of Ca2+ oscillations evoked by multiple neurotransmitters.Figure 4
The dynamic response of astrocyte to presynaptic glutamate and exogenous GABA. (a) The amount of presynaptically released glutamate (Glusyn) decreases as the concentration of exogenous GABA increases, but (b) there is greater activation of mGluRs on astrocyte due to the crosstalk with GABAB receptors, which results in (c) IP3 production and (d) IP3-gated [Ca2+] elevation. (e) Once the [Ca2+] threshold has been reached (0.2 μM, red dashed line), astrocyte releases glutamate (Gluastro) into the extracellular space. To avoid triggering repeatedly the release in all the time steps, gliotransmitter release event is triggered only once in threshold conditions, a mechanism similar to the refractory period of the neuron.


## 3.3. Enhanced Astrocytic Glutamate Promotes Presynaptic Release

The variations in presynaptic release probability in our model were regulated by astrocytic glutamate and exogenous GABA (Figure1). As a contrast, the direct response of presynaptic terminal to exogenous GABA in the absence of astrocytic modulation is shown in Figure 5(a). The release probability Pr significantly decreased with GABA inhibition, while Pr in central excitatory synapses was between 0.1 and 0.9 [61]. This indicates that increased GABA concentrations robustly inhibit excitatory presynaptic release [31]. In the presence of astrocyte-derived glutamate, the scenario could be reversed (Figure 5(b)). In particular, the regulation of presynaptic release by astrocyte-derived glutamate was significantly enhanced as the concentration of exogenous GABA increased. The rate of astrocytic glutamate potentiation was larger than the rate of GABA depression, which gave rise to an increase in Pr, illustrating that glutamate exocytosis from astrocytes promoted presynaptic transmitter release [40]. Furthermore, the variations in excitatory postsynaptic currents (EPSCs) elicited by presynaptic glutamate release are consistent with the quantal size variability in presynaptic vesicle release (Figure 6). The results suggest that exogenous GABA causes a large decrease in the EPSC amplitude [14], whereas the release of astrocytic glutamate enhances excitatory postsynaptic currents.Figure 5
Variations in presynaptic glutamate release probability (Pr) during (a) astrocyte-independent and (b) astrocyte-dependent pathways. The dots represent each presynaptic release event.Figure 6
Variations in excitatory postsynaptic currents (EPSCs) induced by glutamate released from presynaptic terminal during treatment with different concentrations of GABA: (a) 0μM, (b) 1 μM, (c) 5 μM, and (d) 10 μM.


## 3.4. Astrocyte-Mediated Slow Inward Currents Increase Postsynaptic Neuronal Excitability

Both excitatory and inhibitory inputs to the postsynaptic neurons are accompanied by a decrease in impedance and an increase in membrane conductance [33]. The inhibitory and excitatory synaptic conductance values gi and ge were consistently balanced at a certain ratio [62]; the former was several fold higher than the latter [63]. In our model, postsynaptic events were characterized by changes in current instead of changes in conductance. To investigate the impact of NMDAR- and AMPAR-mediated SICs on the postsynaptic neurons, we removed the excitatory effect of regular synaptic currents driven by synaptically released glutamate. These results could be the direct response of the postsynaptic neuron to astrocyte.GABA-evoked astrocytic Ca2+ oscillations triggered the release of glutamate (Figure 7(a)), which diffused into the extracellular space and bound to post-extrasynaptic glutamate receptors, eliciting depolarizing SICs (Figures 7(b)–7(e)). Astrocyte-mediated SICs caused membrane depolarization in the postsynaptic neurons, which was sufficient to reach the firing threshold and ultimately influenced the time course of postsynaptic neuronal firing (Figure 7). Pathological conditions, such as neuronal depolarization, termed paroxysmal depolarization shifts, may contribute to seizure generation [64]. Indeed, intracellular Ca2+ elevation in astrocytes is a crucial factor in the regulation of SICs, which has been previously reported in in vitro [49] and numerical studies [48, 65]. Moreover, the SIC amplitude is determined by the extracellular glutamate concentration. Both an enhanced astrocytic glutamate release and a downregulated glutamate clearance rate can dramatically modulate neuronal excitability, therefore increasing seizure susceptibility [48, 66, 67]. The results of this simulation show that the higher concentration of exogenous GABA is injected to neuronal-astrocytic network, the larger amplitude of astrocyte-mediated SICs is evoked due to the enhanced release of astrocytic glutamate. The increase in excitability could partly counteract the resulting decrease in EPSCs induced by presynaptic glutamate, thereby fine-tuning neuronal network excitation.Figure 7
Astrocyte-mediated SICs and the change in postsynaptic membrane potential. (a) Enhanced glutamate release by astrocytes due to the elevated concentrations of GABA ((b–e)GABAex=0μM,1μM,5μM,and10μM) results in a stronger Isic and subsequently makes the neurons more excitable to resist the increased GABA inhibition.


## 4. Discussion and Conclusions

Numerousin vitro and in vivo studies indicate that astrocytes play a vital role in neuronal excitability and synaptic transmission. Based on neurophysiological findings [20], we developed a biophysical neuronal-astrocytic network model to quantitatively analyze the impact of astrocytes on the modulation of neuronal excitability at different concentrations of exogenous GABA. This is an important issue in view of recent studies regarding the mechanism of the GABAergic neuron to astrocyte signaling, and the simulation results give important clues to the involved GABA signaling mechanism.Our research yielded several predictions. Firstly, the period and amplitude of Ca2+ signaling were dramatically elevated by increasing GABA dose strength. Without exogenous stimulus, i.e., GABAex=0μM, neuronal activity can also trigger relatively slower astrocytic Ca2+ oscillations. In fact, astrocyte Ca2+ events occur in a slow and prolonged manner under physiological condition [68]. This is probably due to the full intracellular cascade of glutamate-activated astrocyte requiring successive events including astrocytic mGluR activation, second messenger IP3 production, astrocyte Ca2+ elevation, and astrocytic gliotransmitter release. Consistent with [20], our results show that exogenous GABA produces a concentration-dependent increase in the elevation rate of astrocyte Ca2+.It is known that astrocytes release glutamate through multiple pathways [69], in which Ca2+-dependent exocytosis may be the most widely examined. Considering the spontaneous astrocytic Ca2+ oscillations [70], neuronal activity, or external stimulus-induced Ca2+ increase [71], there seems to be a consensus that increased Ca2+ concentration results in glutamate release and subsequent neuronal regulation. Yet, Ca2+-dependent mechanism remains a controversial topic. Recent study suggests that Ca2+-dependent gliotransmission is a pharmacological phenomenon rather than a physiological process [72]. Given the complexity of astrocyte Ca2+ activity [73], it is not surprising that the literatures report some discrepant results regarding the mechanisms of glutamate release and gliotransmission. Although the mechanism is not fully settled, a number of published astrocyte computational models have been based on astrocytic Ca2+ signaling [35]. Our model also considers Ca2+-dependent mechanism as the method of glutamate release from astrocytes. Numerical results in this paper verify astrocyte function as the conversion of GABA inhibition into glutamatergic excitation, in which astrocytic Ca2+ activity is a crucial factor contributing to network excitability. Although the release of synaptic glutamate was inhibited by GABA, astrocyte Ca2+ response in this model was enhanced, which may be the amplifying interaction between the GABABR and mGluR-induced Ca2+ signaling [29]. If these are further verified by experimentation, irrespective of which neurotransmitter excites astrocytes, an intervention in Ca2+ signaling may be a potential mechanism for some drug targets. For instance, anticonvulsants reduce neuronal excitability by blocking astrocytic Ca2+ signaling [64].Another significant finding of our model is that GABA and astrocytic glutamate play opposing roles at excitatory presynaptic terminals. We used neurotransmitter depletion and reintegration to describe the modulation of basal release probability by these two factors. The net effect of GABA and astrocytic glutamate on excitatory presynaptic terminals is due to the coexistence of mixed scenarios. The results suggest that the enhanced release of astrocytic glutamate results in stronger vesicle recruitment and release as compared with the inhibitory effect of GABA. From the perspective of information transmission, the probability of presynaptic glutamate release depends on the pattern of presynaptic activity that integrates the contradictory effects of GABA and astrocytic glutamate. Thus, the synapse can be thought of as a filter [74], determining whether synaptic information is transmitted to the postsynaptic neurons.Our model illustrates that astrocyte-elicited SIC by activation of extrasynaptically located postsynaptic NMDA receptors is the other mechanism that affects neuronal activity. Unlike the results of postsynaptic high frequency activity produced by SICs in pathological conditions [48, 66, 75], SIC-mediated increase of neuronal excitability in our study seems to be slightly in modulation of amplitude and period. This scenario may result from integrated inhibitory inputs to the postsynaptic neurons, which are connected to about 20 inhibitory neurons in our model. In fact, SICs can be amplified by upregulating mGluRs on astrocytes, but this was commonly observed in pathological tissues such as epilepsy [76]. Although it is hard to fight against the extra amount of GABA, these results evidence the involvement of astrocytes in modulation of neuronal activity, suggesting that astrocytic glutamate mediates the increase of excitation in the neural-astrocyte network.Overall, the model described in the present study attempts to mimic the bidirectional neuron-astrocyte interaction under exogenous stimuli. These results demonstrate the involvement of astrocytes in synaptic transmission and suggest the excitatory drive of GABA-activated astrocytes in the neuronal network. Indeed, there exists a fine balance between excitation and inhibition in the central nervous system. Astrocytes have the potential ability to modulate synaptic activity and maintain this excitatory-inhibitory balance in particular brain states, which is far-reaching within brain function.



---

*Source: 8750167-2020-09-15.xml*