The Roles of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slow-wave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking (RS) cells with frequency adaptation and do not exhibit bursts in current-clamp experiments (in vitro). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings (gsyn) and external current inputs (I). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to the bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold (IL) and low-threshold (IT) ion channels play a role in increasing and decreasing the parameter conditions (gsyn and I) in which bistable dynamics occur, respectively. For high values of IL conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance IT increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.


Introduction
The classical Hodgkin-Huxley model, proposed in 1952, gave perspectives on understanding the electrical properties of neuronal dynamics considering the biophysical mechanics of the exchange of charges by sodium and potassium ion channels [1,2].Nowadays, much more is known about the variety of channels involved in the neuronal machinery [3,4] and their relation to healthy brain function [5,6] and diseases [7,8].Despite that, the complete understanding of the role of such ion channels is still inciting many research questions in this matter [9,10].Currently, an extensive perspective of research considering the neuroinformatics approach provides new insights into understanding and treating brain diseases [11,12].Neuroinformatics has been used to fit and explore the dynamics of realistic neuron activities [13], including the influence of ion channels [14] and different currents to which neurons can be exposed [15,16].In the present work, we are particularly interested in the ion channel influence as an adaptive mechanism that affects the firing pattern and synchronization transitions since adaptation plays an important role in synchronization [17,18].
The transition of neuronal single-cell firing patterns can emerge as a sum of different sources.In this framework, topology properties, such as the number [18,19], the type [20], the intensity of neuronal interactions [21][22][23][24][25], and the time delay of chemical transmission [26,27], can contribute to the firing patterns transitions.Besides that, noise input currents, as well as external perturbations, can also play a role in the spike and burst transitions [28][29][30].The relationship between how single spikes and bursts are maintained in relation to theta and gamma oscillations was recently discovered [31].However, beyond the topology factors, intrinsic neuronal properties have been pointed out as the main factor in firing pattern transitions [32,33].For example, the mechanism of the adaptation of excitability appears to be a key factor in the transition between spike and burst patterns [34,35].Typically, spike adaptation corresponds to the capability of neurons to reduce their spike frequency due to recent sub and over-threshold neuronal activities [36][37][38].It is also known that potassium and calcium currents develop a role in neuronal adaptation mechanisms [39][40][41].
The emergence of different network firing patterns is mainly associated with the intrinsic properties of the neurons, for example, ion channel density [33,42] and type [43,44], as well as the synaptic input currents [45,46], noise [47,48], and other kinds of couplings and external inputs that the neurons of the network can submit [49].A particular type of populational neuronal firing pattern in network dynamics is the synchronized state.In this context, the influence of the ion channels in synchronization was demonstrated by Boaretto et al.It was shown that small changes in ionic conductance can affect the capacity of the networks to exhibit phase synchronization [50].Landenbauer et al. studied the impact of adaptation on synchronization by comparing coupled neuronal models, and highlighting the importance of the comprehension of such a mechanism on the emergence of the network dynamics [51].Synchronized activities are essential for the correct functioning of the brain; however, highly synchronized activities with burst firing patterns are associated with epileptic seizures [52].Besides the high synchronization levels, epileptic activities can also be related to burst firing patterns [53,54] and low synaptic inhibition [55].The experimental results of cultured neuronal networks, based on multi-electrode arrays, indicated that burst activities present stronger synchrony capabilities with a sufficient level of excitation [56].These currents are also known to develop frequency-dependent resonant mechanisms, an example being the hyperpolarization-activated I h current, which is related to K + and connected to theta resonance and the T-type voltage-dependent Ca 2+ channel that acts as an amplifier [57].In addition, T-type channels are related to the increase in burst activities in seizure generation [58] and have been identified to play a central role in epileptogenesis in the pilocarpine model of epilepsy [59].Besides that, after status epilepticus, the density of T-type Ca 2+ channels was upregulated in pyramidal neurons [60].However, a deeper understanding of the complementary roles of channels in the complex firing patterns of neurons and networks has yet to be uncovered.
Bistable patterns are dynamic behaviors that can be exhibited from the levels of a single neuron [61] to a network [62,63].Bistability in synchronization emerges as a collective dynamic interplay between excitatory and inhibitory interactions between neurons [64].Recently, Akcay et al. reported that bistability with different phases can emerge in a pair of Type I neurons connected by chemical synapses [65], a phenomenon related to the history dependency of the system and known as hysteresis.Bistability and hysteresis are two mechanisms that can be associated with the emergence of burst oscillations [66].In cortical neurons, the state-dependent coexistence of tonic and burst firings gathers the conditions for the emergence of bistability and hysteresis [67].Simplified neuronal and oscillator models have reported the emergence of bistable firing patterns [34,68].Bistable and multi-stable dynamics are said to play a role in both healthy [69] and abnormal brain activities [70].For this reason, studying a detailed description of ion channels in biophysical neuronal models can bring new light to understanding such firing patterns that emerge in brain activities.
The main purpose of this work is to investigate the impact of potassium and calcium currents as an adaptive mechanism that enables the emergence of burst synchronization associated with a bistable regime [34].In our simulations, burst activities emerge with highly synchronized firing patterns for a range of synaptic couplings (g syn ) and external current input (I).Slow potassium currents promote the emergence of high-synchronous activities and spike-to-burst firing pattern transitions.This transition is bistable where physiological asynchronous states coexist with pathological burst synchronization.The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by synchronized firing states; however, asynchronous states result in an asynchronous spike firing pattern.Furthermore, we notice that for high values of highthreshold (I L ) conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input.On the other hand, when the conductance I T increases, the opposite is observed.As our main conclusions, we show how the dynamic and biophysical characteristics of the slow potassium and calcium currents in networks promote bistable transitions.We believe that this work has the potential to uncover pharmacological targets in a manner in which high synchronization can be efficiently hindered.
Our paper is organized as follows.In Section 2, we describe the fundamental equations that govern the dynamics of the model and the diagnostics considered for the analysis of neuronal dynamics.In Section 3, we present our results.We depict the neuronal dynamics by first presenting the perspective of a single network and, following that, presenting the network configuration.We gradually introduce the role of potassium and calcium currents by observing how the emergence of the bistable regime takes place.Finally, in Section 4, we present our discussion to expose our future perspectives and conclusions of this work.

Neural Model
We consider a conductance-based model in which the membrane potential V [71] is given by the following equation: where C m = 1 µF/cm 2 is the specific capacitance of the membrane, g leak is the resting membrane conductance, E leak the reversal potential, and I ionic is the sum of partial ionic currents I j .The voltage-dependent ionic currents have the same general equation: where the j-th ionic current I j is expressed as the product of the maximum conductance of each ion j with the conductance g j of the respective ion.The variables m and n are the activation variables of sodium and potassium, respectively, and h is the inactivation variable of the sodium ion channel [1].The difference between the membrane potential V and the reversal potential for a specific ion E j is (V − E j ) [71].

Description of Ionic Currents
The total ionic current I ionic described in Equation ( 1) is given by: where I Na and I K are the basic sodium and potassium currents of the Hodgkin-Huxley model [72], I M is a slow voltage-dependent potassium current responsible for spike frequency adaptation [73], I L is a high-threshold calcium current, and I T is a low-threshold calcium current [74,75].

Sodium and Potassium Currents
The voltage-dependent sodium and potassium currents are described by the Hodgkin-Huxley equations and were adapted for central neurons by Traub and Miles [72].The sodium currents are described in the following equations: , and the potassium currents are described by The conductance of sodium and potassium are g Na = 50 mS/cm 2 , g K = 5 mS/cm 2 , and the reversal potentials are E Na = 50 mV, E K = −100 mV, respectively.The V T variable is used to adjust the peak threshold in our simulations: V T = −55 mV.

Slow Potassium Current
The non-inactivating slow current of potassium ions is described by the equations: where g M = 0.03 mS/cm² and τ max = 1000 ms [40,73].

Calcium Currents
The first calcium current used to produce bursting is due to a high-threshold Ca 2+ current (I L , L-type, and long-lasting current) and is described as: where the maximum conductance of I L is g L = 0.3 mS/cm 2 [75].
The equations for the second type of calcium current, the low-threshold Ca 2+ current (I T , T-type, transient current), are: where g T = 0.4 mS/cm² is the maximal conductance of the low-threshold calcium current, and V x = 2 mV is a uniform shift of the voltage dependence [76,77].Calcium currents change intracellular calcium concentration [Ca 2+ ] i , and as a consequence, the potential reversal of calcium ions (E Ca ), which is given by where R = 8.31 J K −1 mol −1 is the universal gas constant, T = 309.15K is the temperature, and F = 96485 C mol −1 is the Faraday constant.The dynamics of the calcium concentration is given by where ] ∞ is the maximum concentration of calcium inside the cell, d = 1 µm and τ r = 5 ms.

Neuronal Network
We consider a randomly connected neuronal network composed of Hodgkin-Huxley neurons described in the previous sections, where the neurons are 80% excitatory and 20% inhibitory [78].We consider an Erdös-Rényi network with a connection probability p = 0.1, and there are no auto-connections [79].For the representation of each neuron i in a network, Equation (1) with the addition of synaptic connections is represented by where V i represents the membrane potential, I i ionic corresponds to the ion currents, I is a constant current that is the same for all neurons, and A = 0.2895 ×10 −3 cm 2 represents the membrane area.The neuron area taken as a cylinder is given by A = πdL, with d = L = 0.0096 cm [40].The specific capacitance of the membrane is C m = 1 µF/cm 2 ; therefore, the total capacitance is 0.2895 ×10 −3 µF.All the terms in Equation ( 12) have units of current (µA).The parameters g leak , g Na , g K , g M , g L , g T , and g syn are density conductances with mS/cm 2 as the unit.
The chemical synaptic current that arrives in the neuron i is represented by where V i rev represents the synaptic reversal potential, M is the adjacency matrix of the connections, and g k is the synaptic conductance from the neuron k.V i rev assumes a value equal to V i rev = V exc rev = 0 mV for excitatory and V i rev = V inh rev = −80 mV for inhibitory synaptic connection from the neuron k.The adjacency matrix is composed of element M ij equal to 1 to represent connections from neuron k to neuron i, and equal to 0 to represent the absence of such a connection.g k is updated in the time of neuron k overpass V = 0 with a positive potential derivative (dV/dt > 0).The time of spikes of a certain neuron k, denoted as t k , is also defined by these two conditions in the membrane potential.In this way, the update in synaptic conductance is represented by g k → g k + g syn , where g syn is the chemical intensity of the synaptic updates, the same for excitatory and inhibitory connections.In addition to the update rule due to spiking, each synaptic conductance g k evolves by an exponential decay described by dg k /dt = −g k /τ syn with τ syn = 5.0 ms.

Measures 2.3.1. Firing Rate
We calculate the mean firing rate (in Hz) of all neurons in the network by where t i is the time of the i-th spike, N = 1000 is the number of neurons, and T = 5 s is the time window considered for the analyses.

Coefficient of Variation
We use the inter-spike intervals (ISIs), where the i-th interval is defined as the difference between two consecutive spike times t i+1 and t i , namely ISI i = t i+1 − t i > 0. From the ISI series, the first interval is referred to as ISI 1 followed by the subsequent intervals, namely ISI 2 , ISI 3 , . . ., and ISI n .The ratio between the standard deviation and the mean (indicated by • ) gives rise to the coefficient of variation for the i-th neuron.Finally, the average of CV i over all neurons is given by

Synchronization
To measure the level of synchronization exhibited by the neuronal network, we consider the complex order parameter of Kuramoto [80] given by where the phase of each neuron j is represented by where t j,m represents the m-th spikes of the neuron j.The time t of the parameter is defined in the interval t j,m < t < t j,m+1 .The time-average order parameter is calculated by where t ini and t fin represent the initial and final times for the analyses of synchronization in the neuronal network [81].The implementation of the numerical simulations was performed using self-developed C and NetPyNE codes [82] and can be freely accessed in https://github.com/FernandoSBorges/BistabilityHH(accessed on 1 September 2023).

Neuron Dynamics
We begin by presenting some essential dynamical characteristics of the neurons and networks.Figure 1 displays the biophysical properties of neurons and their dependence on the ionic conductance of the slow potassium M-current, and Figure 2 displays the additional presence of calcium currents.In Figure 1, we present extended characteristics for the regular spiking model by changing the input amplitudes.This model includes I Na and I K for the generation of spikes, and slow K + current (I M ) for the adaptation of the spike frequency.The frequency of the action potentials increases with the input amplitude as shown by the three exemplar voltage traces (Figure 1A).For input amplitudes of I = 110 pA, the neuron exhibits a single spike after ≈ 350 ms from the start of stimulus (Figure 1B).The second spike occurs for I = 130 pA, where the first frequency F is obtained by 1/ISI 1 .Regular spiking behavior is observed for I > 135 pA (Figure 1C).As the amplitude increases, no bursting appears but rather a linear increase in frequency with progressively lower frequencies for each spike in the train as a consequence of the adaptation mechanism (Figure 1C).
In Figure 1D, we show F (in color) as a function of g M and I, where F (Equation ( 14)) is calculated considering spikes in a time window of 5.0 s.The g M has a great influence on the way the neuron fires by changing the minimum value of I, where F > 0 (white line), and in addition, higher values of g M have lower values of F for the same value of I.This shows that the slow K + current (I M ), related to the adaptation of the spike frequency, has a strong influence on the response of the neuronal firing to external stimuli.On the other hand, low-threshold calcium currents (I T ) have low alterations in F (yellow line), and high-threshold calcium currents (I L ) have an almost null effect (brown line).In order to understand the role of calcium currents in the neuron model, we show some diagrams that vary the low-threshold (g T ) and high-threshold (g L ) conductances.We delimited the range of the parameter based on the firing rate values lower than 20 Hz. Figure 2A-D display colored (g T , g L ) and (g M , g L )-diagrams for F and CV. Figure 2E exhibits some examples of selected simulations.The initial burst extends into sustained bursting due to the influence of high-threshold Ca 2+ ; this phenomenon becomes evident for g L > 0.026 mS/cm 2 (Figure 2B).This observed frequency remains below 12 Hz for g M > 0.03 mS/cm 2 and g L < 0.025 mS/cm 2 (Figure 2C), which is typically found in neurons in the rat somatosensory cortex [83].The coefficient of variation (CV) undergoes a transition from 0.3 to >1.0 depending on the value of g L for g M > 0.022 mS/cm 2 .This CV transition is abrupt and occurs at approximately g L ≈ 0.25 mS/cm 2 (Figure 2D).The two main effects of increasing g T are a reduction in the time to the first spike and to generate the initial bursting pattern (Figures 2E and S1).For some values of the g L and g M combination, regular bursts appear.Furthermore, the amplitude of the input current is also responsible for the change in the area in which regular burst activity (CV > 1.0) is observed in the colored diagram.Higher values of I increase the minimum g M values to have bursts, for I = 200 pA the minimum g M = 0.022 mS/cm 2 , while for I = 250 pA, g M > 0.038 mS/cm 2 is necessary (Figure S2).
The panels in Figure 2 provide clear evidence of the impact of adaptation and burst development resulting from the influence of slow potassium and calcium currents.In the next section, we will move forward with the network effects.

Neuron Networks
Next, we study the behavior of networks of neurons, such as the ones discussed above, connected via chemical synapses.As such, the type of behavior varies depending on the input current I and the strength of the synapses g syn .A highly interesting phenomenon that arises in the network spiking pattern is the transition from asynchronous spiking to burst synchronization.Moreover, in this transition, bistable firing patterns are observed, where both asynchronous and synchronous states coexist.
In Figure 3, we present a systematic study of the parameter combinations of constant applied current and the chemical synaptic conductance (I and g syn ) in the measures F, CV, and R (see Section 2.3).These measures are displayed, respectively, in Figure 3A-C.Notice that here, we purposely remove the effect of g M,L,T to show that without the slow potassium and calcium currents, there is no burst no matter the input amplitude and coupling (Figure 3B).Moreover, the firing rate increases with I and g syn (Figure 3A).Interestingly, however, as seen in Figure 3C,D, even without the effect of g M,L,T , the raster plot shows a synchronized behavior with activity around 4 Hz due to population dynamics.This is noticed for combined low values of input current I and g syn , whereas for higher values, only asynchronous behavior is observed (Figure 3E).The results demonstrate that the lack of slow potassium and calcium currents hinders the possibility of observing single-cell bursting, even in a network.Transitions from asynchronous to synchronous activity are still present as seen in how the circle and square are in areas of different values of R in Figure 3C, but the CV changes are slight and nearly absent, indicating only spikes.

Bistable Regime
The results are different when the effect of the slow potassium current is added.As shown in Figure 4, increasing g syn results in transitions from asynchronous activity to a burst synchronization.Noticeably, the voltage traces exhibit individual bursting for neurons with a regular spike firing pattern without coupling.These network bursts are observed for fixed values of I < 200.0 pA and g syn > 1.0 µS/cm 2 .(F-H) Curves of F, CV, and R, respectively, for increasing values of g syn (blue triangles) and decreasing values of g syn (red triangles).We use the model with slow potassium and without calcium currents, i.e., g M,L,T = [0.03,0, 0] mS/cm 2 .
The emergence of the bursting synchronization by increasing values of g syn has a marked effect on the measures F, CV, and R. When I = 170.4pA, this sudden increase that leads to the transition starts approximately at g syn = 1.09 µS/cm 2 for initial conditions with asynchronous spikes, and g syn = 1.18 µS/cm 2 for initial conditions with burst synchronization.This hysteresis is indicated in the white area in Figure 4A-C and is an observation that is a strong indication of bistability.Figure 4D shows a raster plot and voltage traces for the asynchronous spike pattern, and Figure 4E shows the same for burst synchronization.A comparison of the blue and red curves in Figure 4F-H shows the exact transition area from asynchronous activity to bursting synchronization.
The results demonstrate that the slow potassium current (here considered as g M = 0.03 mS/cm 2 ) can promote bistable dynamics in neuronal networks.This type of observation related to a biophysically grounded parameter is a key factor in understanding how bistability relates to other brain phenomena, such as decision making or pathologies such as epilepsy.The transition from spike to network burst was observed in [18] for the adaptative exponential integrate-and-fire neuron model, and, later, bistability in neuronal networks using this simple neuron model was related to epileptic seizures elsewhere [34,84].

Calcium Effects in the Bistability
We showed that a slow potassium current is necessary to observe bistability.A natural question that arises from our work is how the firing transition and the bistable dynamics depend on the single-cell calcium ionic channels.In this section, we discuss how quantitative changes emerge for different combinations of calcium currents.Figure 5 shows how CV values depend on constant current and chemical synaptic conductance (I,g syn ) with respect to different values of g M,L,T .The bistable parameter region is identified in white and separate black and reddish/blueish regions for asynchronous spikes (low CV) to synchronous bursts (high CV), respectively.Figure S3 demonstrates that these qualitative patterns are maintained for the firing frequency F and the mean order parameter R. Indeed, the transition area changes according to the presence of ionic currents.The parameters g syn and I are sensitive to changes with respect to g M,L,T .In particular, highthreshold calcium, controlled by g L , is the one that causes a greater change by allowing transitions at lower values of the g syn (Figure 5C).Furthermore, the synchronous area (reddish/blueish) is larger for the high-threshold calcium change, allowing for burst synchronization for I > 200.0 pA.In contrast, low-threshold calcium promotes the opposite effect by slightly increasing the value of g syn necessary to observe a transition (compare Figure 5A,B and Figure 5C,D).Figure S4 shows the dynamic behavior of currents generated by each active conductance in a neuron during the bistable regime.The analysis includes the total synaptic current (I syn ), the total ionic current (I ionic ), the fraction of outward ionic currents, and the fraction of inward ionic currents.See [85] for a detailed analysis of currents in neural models with different conductance densities.
Our results show that calcium has a double effect in promoting the bistability of asynchronous/synchronous activity that we observed: it can either facilitate it by allowing lower values of g syn to act on the transition when the high-threshold calcium channel is used or make it harder via an opposite effect mediated by the low-threshold calcium channel.

Discussion and Conclusions
In this work, we extend the analysis of bistable firing transitions in neurons and networks [34,[61][62][63] considering the effect of potassium and calcium currents.To quantify these firing patterns, we employ the usual measures, including the firing rate, CV, and synchronization level.This provides an important link to the role of these ionic currents in controlling the network behavior.We notice that in the absence of the joint effect of potassium and calcium, there is only a spike synchronization pattern, and bursts at the single-cell level are not observed.In the presence of slow potassium and calcium currents, bistability is identified by hysteresis among synchronized bursting and asynchronous activity.It is worth noting a significant distinction from the study conducted by [34], which utilized the adaptative exponential integrate-and-fire neuron model.In our current model, we not only observe the effects but also establish a direct link to the underlying neurobiology, enabling us to propose experimental interventions that can be tested empirically through the manipulation of ionic currents.
Different voltage-dependent ion channels are responsible for the control of excitability and can play an intrinsic role in pathologies, such as epilepsy [86][87][88], Parkinson's [89], and Alzheimer's disease [90].In particular, potassium and calcium ions have been found to play a role in many of these diseases [91][92][93].In addition, some research has evidenced that blocking such ion channels can play a role in avoiding the mechanism associated with the emergence of epileptic activities [94].Detailed investigations have been devoted to identifying the variations and mutations of calcium and potassium ion channels that exert the most influence on epileptic activity [95].
The opening and closing of such ion channels can have different effects on neuronal excitability.The majority of potassium ion channels open when the membrane depolarizes and close when it hyperpolarizes.The M current, in particular, acts in the subthreshold domain and limits the ability of the neuron to fire repetitively [96].Their involvement in epileptic activities is still being investigated [97,98], and some recent research studying the lost-and gain-of-function of such ions has started to clarify these issues [2,99].Besides that, mutations in potassium channels have also been identified as an important factor in the parthenogenesis of human epilepsy [100].Such mutations and anomalies can be generated by external factors, e.g., the use of drug substances that have evidenced potential for changing the expression of potassium channels [101].Recently, precise studies and therapies have focused on the mutation of potassium channel genes [102].It is clearly important to further investigate how channel-blocking mechanisms can serve to control network activity [103,104].
Regarding calcium, we have shown that neurons embedded with these channels can exhibit a spike-to-burst transition [105].In particular, two types of calcium currents are highlighted: L-type and T-type.The former, where "L" stands for large or long-lasting, is a high-voltage activated channel.The latter, where "T" stands for transient, is a low-voltage activated channel.Calcium currents are highly involved in network high synchronous patterns that are observed in epileptic seizures [106].Not surprisingly, calcium channel blockers may act as anti-seizure drugs for prevention and treatment [107,108], i.e., the blocking of such ion channels can attenuate the burst firing pattern [109] and hinders epileptic depolarization [110].Nonetheless, it is challenging to develop drugs that are subtype-specific, and in most cases, they have not been developed yet [111].Thus, understanding the spike-to-burst pattern through computer simulations and neuroinformatics has the potential to elucidate the emergence of synchronous patterns related to epileptic activities [112].A combined treatment using calcium channels blockers was demonstrated to be beneficial in anticonvulsant and antinociceptive effects [113] (for the reader interested in channel blockers, please see the list in [111]).According to our results, whereas L-type channels make transitions easier by lowering the value of synaptic strength (g syn ) necessary for the transition, T-type channels have the opposite effect.
Admittedly, there is yet a joint effect of potassium and calcium currents that we have not addressed in our simulations, which is the one provided by Ca 2+ -dependent K + currents (I K(Ca) ) [114][115][116].These potassium currents, influenced by calcium concentration, are known to induce prolonged hyperpolarization pauses, thereby significantly impacting the firing properties of neurons.The exploration of their intricate behavior is a key aspect of our future research plans.In future work, we will also analyze the effect of the hyperpolarization-activated cation current (I h ) [117] and the persistent Na + current [118] on the bistability.
Our results also extend the current knowledge of the joint effect of calcium and potassium ion channels in the context of other important firing pattern transitions.Based on our analysis, we can predict the ionic blockers required to avoid and treat high synchronous activity.These blockers would change the single-cell bursts that may result in epileptic seizures to asynchronous activity.In particular, we have evidence of the role of such ion channels not only in the spike-to-burst firing pattern transition but also in the relation of such transition with the synchronous patterns.Furthermore, we can gain insights by examining situations where the channels are blocked selectively, rather than all simultaneously.

Figure 1 .
Figure 1.Model of regular spiking neuron, with I Na and I K for spike generation, and slow K+ current (I M ) for spike frequency adaptation.(A) (Top) Voltage traces with different amplitudes of the depolarizing pulses (bottom), V = −85 mV before the depolarizing pulses.(B) Time to first spike (t 1 ) as a function of the injected current (amplitude of the pulse I).(C) Frequency-current curves (F/I), where the instantaneous firing rate (inverse of the inter-spike interval) is represented as a function of I.The curves indicated by different colors correspond to the 1st, the 2nd, the 3rd, and the 10th spikes in the train.(D) Spike frequency F (in color) as function of g M and I, considering 5.0 s time window.The white line represents the transition where F > 0 for g T = 0 and g L = 0. Additionally, this transition lines for g T = 0 and g L = 0.1 mS/cm 2 (brown line), and g T = 0.4 mS/cm 2 and g L = 0 (yellow line) are shown.Other parameters are L = d = 96.0µm, g leak = 0.01 mS/cm 2 , E leak = −85.0mV, g Na = 50 mS/cm 2 , V T = −55.0mV, g K = 5 mS/cm 2 , τ max = 1000 ms, and g M = 0.03 mS/cm 2 .The absolute values of the membrane potential of (A) are available at https://github.com/FernandoSBorges/BistabilityHH(accessed on 1 September 2023).

Figure 2 .
Figure 2. Firing pattern for different I M , I T , and I L conductances.(A) Firing rate in colored (g T , g L )diagram for g M = 0.03 mS/cm 2 .(B) The same as (A) for the CV.(C) Firing rate in colored (g M , g L )diagram for g T = 0.4 mS/cm 2 .(D) The same as (C) for the CV.(E) Exemplar voltage traces considering different values of g M , g L , and g T , where each parameter combination is shown atop and V = −85 mV before the depolarizing pulses.Other parameters are the same as Figure 1 with I = 200 pA.

Figure 3 .
Figure 3. Spike synchronization in space parameters as a function of I and synaptic conductance (g syn ).(A) Mean firing frequency.(B) Coefficient of variation, and (C) mean order parameter.(D,E) Raster plot and the neuronal membrane potentials for the parameter indicated in (A-C).(D) Synchronized spikes for the parameters I = 88.3 pA and g syn = 0.2 µS/cm 2 .(E) Desynchronized spikes for the parameters I = 98.9 pA and g syn = 0.5 µS/cm 2 .We use the model without slow potassium and calcium currents, i.e., g M,L,T = [0, 0, 0].

Figure 4 .
Figure 4. Burst synchronization in space parameters as a function of I and chemical synaptic conductance (g syn ).(A) Mean firing frequency.(B) Coefficient of variation.(C) Mean order parameter.(D,E) Raster plot and the neuronal membrane potentials for the parameter indicated in (A-C).(D) Desynchronized spikes for the parameters I = 170.4pA and g syn = 1 µS/cm 2 are observed and synchronized bursts are shown for the parameters I = 170.4pA and g syn = 1.3 µS/cm 2 in (E).(F-H)Curves of F, CV, and R, respectively, for increasing values of g syn (blue triangles) and decreasing values of g syn (red triangles).We use the model with slow potassium and without calcium currents, i.e., g M,L,T = [0.03,0, 0] mS/cm 2 .

Figure 5 .
Figure 5.Effect of calcium currents in the burst synchronization in space parameters as a function of input current (I) and chemical synaptic conductance (g syn ).(A-D) Coefficient of variation (CV) for different combinations of g M,L,T (see values atop).(A) Model with slow potassium and without calcium currents, the same parameters of Figure 4. (B) Model with slow potassium and only lowthreshold calcium current (I T ).(C) Model with slow potassium and only high-threshold calcium current (I L ). (D) Model with slow potassium and both low-threshold and high-threshold calcium currents (I T and I L ).