Abstract

Structural rearrangement of the dentate gyrus has been described as the underlying cause of many types of epilepsies, particularly temporal lobe epilepsy. It is said to occur when aberrant connections are established in the damaged hippocampus, as described in human epilepsy and experimental models. Computer modelling of the dentate gyrus circuitry and the corresponding structural changes has been used to understand how abnormal mossy fibre sprouting can subserve seizure generation observed in experimental models when epileptogenesis is induced by status epilepticus. The model follows the McCulloch-Pitts formalism including the representation of the nonsynaptic mechanisms. The neuronal network comprised granule cells, mossy cells, and interneurons. The compensation theory and the Hebbian and anti-Hebbian rules were used to describe the structural rearrangement including the effects of the nonsynaptic mechanisms on the neuronal activity. The simulations were based on neuroanatomic data and on the connectivity pattern between the cells represented. The results suggest that there is a joint action of the compensation theory and Hebbian rules during the inflammatory process that accompanies the status epilepticus. The structural rearrangement simulated for the dentate gyrus circuitry promotes speculation about the formation of the abnormal mossy fiber sprouting and its role in epileptic seizures.

1. Introduction

Epileptic syndromes are a group of neurological disorders with different etiology and clinics characterized by recurrent seizures with excessive and hypersynchronous neuronal activity [1]. Nowadays, approximately 50 million people are affected with epilepsy and 40% of which are classified as temporal lobe epilepsy (TLE) [2].

The human TLE is frequently associated with neuronal loss in the hippocampus and dentate gyrus (DG) [3]. Throughout this degeneration, called hippocampal sclerosis, the granule cells are less affected. The neurotransmitter gamma-aminobutyric acid and proteins expressed in these cells are allegedly responsible for neuroprotection, increasing resistance to excitatory injury [4]. The DG has been proposed to be a gate for entry of intense neuronal activity into the hippocampus. However, in TLE, the synaptic reorganization involving the DG granule cells (abnormal sprouting of hippocampal mossy fibers) may act to reduce the DG filtration properties, increasing the propensity to intensify the neuronal activity and, therefore, the seizure propagation through the hippocampus and other limbic system structures [5].

Mossy fiber sprouting (MFS) has been observed in several animal models of TLE. The recurrent synaptic transmission between granule cells is one of the supposed inductors of a recurrent excitation between these cells and therefore is responsible for the generation of spontaneous seizures with probable origin in the hippocampus [8, 9]. Some investigators believe that the MFS makes the hippocampal circuitry susceptible to prolonged bursts of action potentials, contributing to the seizure onset [9, 10]. Thus, the abnormal MFS could play an important role in sustaining seizures in TLE with mesial sclerosis, contributing to the provision of the epileptogenic substrate of the hippocampus. According to Lytton et al. [11], the abnormal organization of the brain circuits can be considered crucial for several types of epilepsy.

On the other hand, the occurrence or lack of occurrence of spontaneous recurrent seizures in an animal model is not directly related to the presence of mossy fiber sprouting. Therefore, the complexity involved in the investigation of the main mechanisms responsible for seizure eruption makes the computational simulation of DG an indispensable approach in providing new insights. Of course any model is an oversimplification of the reality. In the particular case of a mathematical model, the reductionism must be a rule to make the computational implementation reliable. However, it is a powerful tool to help evaluate if the minimum conditions represented are enough to reproduce the phenomenon investigated.

In 1943, McCulloch and Pitts created the first formalized model of an artificial neural network. Since then, computer modeling has enabled important advances in the understanding of complex phenomena such as epilepsy [1214].

The aim of this study was to investigate the possible role of abnormal MFS on the DG excitability, involved in experimental models of status epilepticus. The computational model was based on the McCulloch-Pitts formalism and neuroanatomical data as well as on patterns of connectivity of the main cell types present in the hippocampus [12, 1519]. The model also represents sustained depolarizations that are typical of nonsynaptic mechanisms during the ictal periods and changes in GABAa transmission reported during inflammatory response observed after pilocarpine-induced seizures. The simulations create the possibility for discussion concerning the complex role of the synaptic reorganization associated with experimental models of status epilepticus.

2. Methods

2.1. Neuronal Network Model

The neural network implemented for the simulation study is based on the McCulloch-Pitts formalism [15, 20], including the Compensation Algorithm for synaptogenesis, proposed by Dammasch et al. [19, 21], and Hebbian and anti-Hebbian rules. The cell population represents the main types described in the dentate gyrus (DG) as well as its connectivity (Figure 1). Since the cytoarchitecture of DG is favourable to the non-synaptic connections, they were also represented in the model.

The neuronal network, with neurons, is composed of NE excitatory and inhibitory neurons. The excitatory neurons were divided into granule cells and hilar mossy cells. The connectivity is represented by a connection matrix , where is the synaptic weight of the connection of the neuron on the neuron .

The collection of states of each neuron of the network is represented by a vector , where . Each neuron may have two possible states in the instant of time : active () or inactive (). The probability of neuron turning active in time could be estimated by the threshold function [20]: where represents the stimulation input, is the membrane potential of the neuron , is the threshold, and represents the noise of the threshold function.

The membrane potential in the instant of time of the neuron , , is calculated by where is the effect of the non-synaptic mechanisms on the transmembrane potential. The connection weights are (connections from granule cells), ( and ; connections from mossy cells), and (inhibition from interneurons).

According to Almeida et al. [14, 22], the Na+ influx trough ionic channels and the Na+ efflux trough Na/K pump define the transition points between the ictal and interictal periods during non-synaptic epileptiform activity. In the present model, the ictal period of the epileptiform activities was associated with the spontaneous increase of cell activities (SICA). Representating these two Na+ currents in the generation of a SICA, was calculated based on two components: NSE, representing the excitatory effect of the Na+ influx, and NSI, representing the inhibitory effect of the Na+ efflux. Therefore, where “” and “” are parameters used to adjust the balance between excitation and inhibition, respectively.

In the proposed model, to simulate the ictal and interictal periods, two states of the neuronal activity were identified: sustained depolarization state and polarized state. The former state comprises the ictal period and the last the interictal period. In order to represent the transition between both states, the sustained depolarization state was divided in two: the sustained depolarization itself and the in repolarization state. Defining as the average activity of the cell , as described in (9), and as the threshold for the action of the non-synaptic mechanisms, the sustained depolarization happens when . During the sustained depolarization, the neuron is firing and the sustained depolarization is done due to Na+ influx trough ionic channels greater than its efflux trough pump, which can be represented by . Consequently, when in the in repolarization state is required, the essential condition for the interictal period (polarized state) is . The mathematical description of the three states is basically defined by the same functional composition. The changing rates of NSE and NSI follow the equations where and are dependent on the state and and are time constants governing the recovering rate of NSE and NSI. To Summarize, consider the following(i)Sustained depolarization state ( and ): with constant and dependent on , according to where is constant. In this state, represents an increment of the Na+ influx through the channels, based on the increment of the neuronal discharge. On the other side, the increment of intracellular Na+ increases the Na/K pump action, promoting the inhibitory effect of the electrogenic current, which is represented in the present model by an effect proportional to NSE.(ii)In repolarization state (): where and are constants and and are functions that contribute to cell repolarization: where and are constants. The aim of the above equations was to represent the moment when the electrogenic current of the Na/K pump is able to interfere with the membrane potential promoting the repolarization. In these circumstances, the voltage-dependent sodium channels start to reduce conductance and the repolarization takes place. Equation (9) is a phenomenological representation of this behavior.(iii)Polarized state () with .

2.2. Neuronal Network Connectivity
2.2.1. Compensation Theory

Experimental observations indicate that neurons react to anatomical changes induced by long-term metabolic disequilibrium. These reactions depend on whether the state is high or low [19, 20, 23]. Therefore, each neuron follows the local rules aiming at compensatory changes to the input activity, caused by structural adaptations. Consequently, for a specific neuron , its deviation from the average activity expected, in the present work 0.2, is called the morphogenetic state of the neuron, , which is estimated by means of the following relationship: where the average activity is calculated from to and is the morphogenetic time step, corresponding to the period of time needed to perform the connectivity changes of the network. , used as a measure of the neuron activity, is the average firing probability.

The morphogenetic state of neuron represents the neuronal capacity to form, stabilize, or degrade pre- or postsynaptic elements or synapses [20, 21]. The possible interactions between these three states follow the synaptogenesis compensation theory proposed by Wolff and Wagner [17]. According to this theory, the afferent synapse spectrum influences—via transmission—the functional state of the neuron, estimated as an average of activity along the determinate period of time . Depending on its level of activity, the functional state leads to a morphogenetic state that tries to compensate its deviations. This is performed by changing the neuronal elements, bound and free, pre- and postsynaptic, after each morphogenetic time step. The degeneration of bound elements and recombination of free synaptic elements lead to a new structural state. The new afferent spectrum leads the neuron to a new functional state.

To simulate the network morphogenesis, the synapses are distinguished in pre- (pr) and postsynapses (po), which are then divided into bound (b) and free (f). The postsynaptic elements have additional divisions: excitatory (e) and inhibitory (i) [20]. This classification leads to the following representation: (i)bound (bpr) and free (fpr) presynaptic elements;(ii)bound (bepo) and free (fepo) excitatory postsynaptic elements;(iii)bound (bipo) and free (fipo) inhibitory postsynaptic elements.

The sum of the synaptic elements of the neuron , according to this classification, represented by , is defined as follows [19, 20]: (i)sum of the bound excitatory postsynaptic elements: ,(ii)sum of the bound inhibitory postsynaptic elements: ,(iii)sum of the free excitatory postsynaptic elements: ,(iv)sum of the free inhibitory postsynaptic elements: ,(v)sum of the bound presynaptic elements: ,(vi)sum of free presynaptic elements: .

The morphogenetic rules depend on the morphogenetic state of the neuron , on the actual structural state of the neuron, and also on the cell sensitivity to structural variations [19, 20]. These rules are represented in Table 1.

The current structural state of the neuron is characterized by the number of afferent and efferent synapses. The synaptogenesis and the modification of the synaptic contact can change the structural state of the network [19]. According to the formalization of the compensation theory of synaptogenesis [17], two kinetic parameters are defined: and , where and refer to the level of excitability of the neuron: high () and low (). Changes in the amount of synaptic elements of the neuron , , are calculated based on the dependence of the level of excitability ( and are proportionality constants), on the neuron’s deviation from a desired medium average activity () and on the current amount of synaptic elements (Table 1). In the present model, the changes in the connectivity of the network followed the rules proposed by Dammasch et al. [19] and Dammasch et al. [21].

The network connectivity variation also happens with dependence on the reduction of the bound synaptic elements. The presynaptic elements’ decay is proportional to the strength of the actual connections [19]:

The loss of postsynaptic elements is computed separately for inhibitory and excitatory synapses, according to the following equation [20]: for all .

The presynaptic elements of altered synapses are redistributed postsynaptically while the postsynaptic elements are removed. Therefore, the amount of presynaptic elements previously bound (equal to the amount of degraded postsynaptic elements ) is transferred to the amount of free presynaptic elements [20]:

New synaptic contacts and the reinforcement of the current synaptic contacts take place from the recombination of free elements. The recombination is given by Butz et al. [20]:

for all , where .

The amount of free synaptic elements remained during the recombination is also renewed subsequently.

2.2.2. Hebbian and Anti-Hebbian Rules

The synaptic modification during learning is made by a mechanism that depends on the simultaneous activity of the presynaptic button and the postsynaptic cell. According to the Hebbian learning rule [24], a synapse is potentiated when the activity in the synaptic button increases at the same time that the postsynaptic neuron is depolarized. The synaptogenesis, according to the Hebbian rules depend on the activities of connected cells and on the contribution of each cell to the firing of the other. Therefore, the change of the connection weight between two neurons, and , is determined by where is constant and is the number of times that the discharge of the presynaptic neuron is followed by discharge of the postsynaptic neuron during one morphogenetic time step. Assuming that the formation of new synapses is limited as well as the strengthening of the synaptic weight, in the present model a sigmoid function was used to calculate where is a proportionality constant.

When the postsynaptic terminals are active and the presynaptic neuron remains inactive, the corresponding synapses are depressed. This process is known as anti-Hebbian plasticity, because of its antagonism in respect to the Hebb learning. Thus, the changes in the connectivity were performed according to where is constant and is the number of times that the discharge of the postsynaptic neuron is not preceded by discharge of the pre-synaptic neuron during one morphogenetic time step.

2.3. Connection Matrix and Parameters of the Model

The concept of the connection matrix representing the DG was based on the current literature. According to Dyhrfjeld-Johnsen et al. [16], the hippocampus comprises 1.000.000 granule cells, 30.000 mossy cells, 10.000 basket cells, 2.000 axo-axonic cells, 4.000 MOPP cells (molecular layer perforant path-associated cell), 12.000 HIPP cells (hilar perforant path-associated cell), 3.000 HICAP cells (hilar commissural-associational pathway related cells), and 3.000 IS cells (interneuron selective cells). The model network was performed in a reduced scale (1 : 10.000), that is, with 99 granule cells, 3 mossy cells, and 3 interneurons (basket, axoaxonic, MOPP, HIPP, HICAP, and IS). Combining these three types of cells, nine groups of connections were represented and the connections’ arrangement for the initial conditions of the simulations was set according to experimental investigations [12, 16]: (a)group I (granule cells → granule cells): no connections; (b)group II (granule cells → mossy cells): each mossy cell receives connection from 33 granule cells;(c)group III (granule cells → interneurons): each interneuron receives connection from 33 granule cells;(d)group IV (mossy cells → granule cells): each mossy cell sends a connection to all granule cells;(e)group V (mossy cells → mossy cells): each mossy cell sends a connection to all mossy cells;(f)group VI (mossy cells → interneurons): each mossy cell sends a connection to all interneurons;(g)group VII (interneurons → granule cells): each interneuron sends a connection to all granule cells;(h)group VIII (interneurons → mossy cells): each interneuron sends a connection to all mossy cells;(i)group IX (interneurons → interneurons): each interneuron sends a connection to all interneurons.

The initial values of the elements of the connectivity matrix followed a normal distribution with mean and standard deviation set to .

According to Dammasch et al. [19], Dammasch et al. [21], and Cromme and Dammasch [23], it is hypothetically assumed that degeneration of the free synaptic elements is faster than the degeneration of the bound synaptic elements. Dammasch et al. [21] found in simulations that while the relation between decrease and increase of bound elements is critical for stability, the decrease of free elements is of less importance. Therefore, the following rations of kinetic parameters were chosen:

To maintain the oscillation properties of the network during the morphogenesis, as proposed by Cromme and Dammasch [23], it was assumed that

Assuming no bias decrease of the synaptic elements, then

The relationships of the kinetic parameters described above allow the description of one parameter in terms of the other. Choosing as the adjusting parameter, where can be interpreted as a constant that defines the speed of the morphogenetic changes. During the process of network connectivity change, the kinetic parameters must guarantee smooth convergence in the change of synaptic elements [21] and thus must be small. In this model, studying the effect of the compensation theory on the network activity, was set in the interval . Taking into account the investigations by Dammasch and Wagner [18], Dammasch et al. [19], Dammasch et al. [21], Cromme and Dammasch [23], and Butz et al. [20], it was assumed that , , and .

For the excitatory pre-synaptic cells the synaptic weights were , , when , and , when . These parameters were estimated by running simulations of the network activity without stimulation and allowing it to reach stable activity [19]. The network was considered stable when the average activity reaches values in the range 0.15–0.25.

The parameters related to the increase of the neuronal activity due to the nonsynaptic mechanisms were adjusted to induce depolarization and repolarization during SICA mimicking the effects of the Na+ influx through the channels and its efflux through the pump, during a spontaneous nonsynaptic epileptiform activity, as described by Almeida et al. [14]. With this aim, the adjusted values were , , , , , , , , , , and .

Two timescales were adopted in the simulations, one for the cellular activity changes, corresponding to the duration of each interaction , and a timescale for the connectivity changes (duration of one morphogenetic time step ). Thus, one morphogenetic time step of the proposed model () corresponds to the spiking remodelling and, similarly, to the axonal ramification between different postsynaptic targets [25, 26]. In the model, each morphogenetic time step corresponds to 150 interactions of and each interaction to 0.36 min.

The effect of the Hebbian rule on the connectivity changes and on the network activity was investigated to guarantee smooth connectivity changes. For this conditions, the possible range for was determined.

3. Results

By investigating the effect of different strengths of stimulation (; ) on the neuronal network activity, an increase of induced activity can be seen (Figure 2). It can be observed that the increase of the network connectivity also causes a progressive increase of the network activity with afterdischarges. Mimicking the nonsynaptic mechanisms, during afterdischarges, the depolarization function NSE increases quickly, also raising the value. The inhibition function, NSI, increases more slowly and, at the end of afterdischarges, overcomes the excitation function. From that moment on, the values of NSE and NSI are restored to the levels at which the neuronal activity remains low.

To investigate how different situations of combinations of actuation of the compensation theory and the Hebbian rule act on the network, three simulations were made. The results are shown in Figure 2: (i) without stimulation; (ii) with stimulation; (iii) with stimulation and after the stimulation, assuming the transient change of the interneuron action from inhibitory to excitatory. For the stimulations performed in the simulations (ii) and (iii), the stimulus was applied at the day 1.5, with , during 70 min. The stimulus intensity was adjusted to be less intense but sufficient to induce maximal activation of the network. This procedure took into account the objective of representing the experimental protocols of drug applications in animals to induce epileptiform seizures and, consequently, the status epilepticus [2729]. The transient change of the inhibitory effect of interneurons to excitatory effect aimed at representing the excitatory effect of the GABAa receptors, typical of the inflammatory states in the brain, which can persist from the first 24 h for up to two weeks after the status epilepticus [30]. Experimental observations show that after status epilepticus, at the beginning of the latent period, changes in the transmembrane Cl-gradients may occur, due to anionic accumulation in the neuronal cytoplasm. This accumulation is enough to change the GABAa receptors activation effect from inhibitory to excitatory. The cause of this change was demonstrated to be related to the diminished expression of the cotransporter KCC2, responsible for the Cl-extrusion [30]. This mechanism is reported by Morimoto et al. [7] as the GABAa-mediated excitation hypothesis. This transient change was performed with the synaptic weight of two of the three interneurons altered according to the following equation: where days indicates the onset of the inflammatory state and days the end. Only two interneurons were submitted to the transient changes, since in the DG not all inhibitory synapses are GABAa [31]. The weight of the two interneuron connections was altered for the connection to 22 granule cells. This group of cells corresponds to the 22% of excitatory GABAa observed in epileptic human brain tissue [32, 33]. Similar findings were also encountered by Pathak et al. [30] studying DG granule cells in rat hippocampus.

By comparing the simulation with experimental data the model can be evaluated according to its ability to represent mechanisms subserving the synaptogenesis and epileptogenesis involved in the status epilepticus (Figure 3). The three main periods that characterize the seizure progression after status epilepticus can be distinguished in the simulations: status epilepticus, latent period, and spontaneous seizures. According to Morimoto et al. [7], in the model of status epilepticus the administration of chemical agents such as pilocarpine or kainate is followed by an emergence of continual recurrent seizures sustained for several hours. The same process is represented in the simulations with the network stimulation. Following the induction of status epilepticus, the absence of seizures characterizes the latent (or silent) period. When the seizure returns, the period of spontaneous seizures starts. The simulations involving the stimulation period are in close correspondence with the occurrence of these periods. In Figure 3, the event rates registered experimentally by van Puttena et al. [6] and a diagrammatic representation of the periods, according to Morimoto et al. [7], are correlated with the simulations.

To analyze the SICA induction in the neuronal network, the occurrence of SICA along the 100 days, for all simulated situations, is shown in Figure 4. The SICAs are shown to have dependence on the combination of the parameters related to the compensation theory and the Hebbian rule. To perform the stimulation (Figure 4(b)), the stimulus was also applied at day 1.5, with , during 70 min. The demarked SICAs are the ones whose activity exhibited variations greater than 30%, assuming that the less intense activities are not necessarily related to epileptiform activities. The average activity of the cells and the excitation and inhibition functions showed similar behaviour during the SICA to those observed during the induced increase by stimulation (Figure 2). What distinguishes the two events is the manner of triggering the activity. The SICA deflagration was due to the increased efficacy of the synaptic connections.

In Figure 4(a), when the network is stimulated, it can be observed that the combinations of the compensation theory with Hebbian rule result in SICA for a limited range of values for and . In the case where the network was stimulated (Figure 4(b)), it is evident the increase of region where the SICA emerges, meaning larger range for and . Furthermore, the number of SICA observed for each combination of and was higher. When the GABAa excitation was represented (Figure 4(c)), the conditions favorable for SICA occurred over a wider range of values of and .

In the situation where the network has not been stimulated, it is possible to investigate the average activity of the different cell types (Figure 5(a)) as well as for each group of connections (Figure 5(b)). After 41 days, applying the theory of compensation and the Hebbian rule caused variation in the average activity of granule cells that ranged between 0.0034 and 0.0096% (<1%). This variation is associated with small changes in the connectivity between the granular cells (group I), . For most combinations of and (Figure 5(b), yellow region), there are changes in the connections that the mossy cells and the interneurons send to the granule cells (groups IV and VII, resp.). The major variations of the connectivity of the groups IV and VII occurred for higher values of and lower values of (Figure 5(b), light blue region). Reductions were observed for both groups of connection in the order of . The changes in the average activity of the mossy cells were in the range −0.005 to 0.0068%, being associated with the following changes in the connectivity: group (II), ; group (V), ; and group (VIII), . Similarly for groups (I) and (IV), the major reductions of connectivity occurred for higher values of and lower for (Figure 5(b), light blue region). The change in the connectivity of the interneurons was in the range −0.0196 to 0.0079, and the reduction of the interneuron activity (Figure 5(a), IN, light blue region) can be attributed to the increased role of the compensation theory. The changes in the average activity were associated with the following connectivity changes: group (III), ; group (VI), ; and group (IX), .

When the simulation includes the stimulation (, during 70 min), the percentage variation of the average activity in the 41st day, calculated in respect to the 1st day, was more intense for higher values and lower values (Figure 6(a), dark red regions, for GC and MC, and dark blue for IN). The increased average activity of the GC, in the order of 3% causing the most pronounced number of SICA for and , is associated with the connectivity increase of the groups (I) and (VII), respectively, and , as well as reduction of the group (IV), with (Figure 6(b)). For the same combination of and values, the average activity of the MC also increased almost 3% and the following connectivity changes were observed: group (II), ; group (V), ; and group (VIII), . On the other hand, the average activity of the interneurons decreased 3% and with the following connectivity changes with interneurons as the postsynaptic cells: group (III), ; group (VI), ; group (IX), . When , that is, considering only the compensation theory, stimulation caused, especially for larger values of , reduction in the level of activity and connectivity of the neurons.

Including GABAa changing from inhibitory to excitatory for two interneurons and stimulating the neural network, the granule and mossy cells activities increase more intensely for a larger number of combinations of and values: and (Figure 7(a)). The increase of the average granule cell activities (~5%) increased the SICA occurrence for and and is associated with connectivity increase for the groups (I) and (VII) and connectivity decrease for the group (IV) (Figure 7(b)). For the same and values, the average activity of the mossy cells increased 4.5% and the following connectivity changes occurred: group (II), ; group (V), ; group (VIII), . The average activity of the interneurons reduced ~4% and the following changes in the connectivity occurred for postsynaptic interneurons: group (III), ; group (VI), ; group (IX), . When , therefore considering only the compensation theory, the stimulation followed by the transient excitatory GABAa caused the most intense activity and connectivity reduction (dark blue regions).

The higher incidence of SICA, reproducing the induction of epileptiform activity in experimental models, occurred for and . For these combinations of the compensation theory with the Hebbian rule, the occurrence of SICA is depicted in Figure 8(a). Concomitantly, the inhibitory synapses formation rate increases throughout the network (Figure 8(c)). On the other hand, the connectivity of the excitatory synapses decreases in response to the activity increase (Figure 8(b)) due to the compensatory theory (Figure 8(c)) and the anti-Hebbian rule (Figure 8(c)). Following the afterdischarges induced by the network stimulation, between the second and twelfth days, the average activity of the network reduces (Figure 8(a)); it is when a progressive increase of the excitatory synapses takes place (Figure 8(b)—green and blue curves). These changes were correlated with the synaptic changes, with progressive increase of the excitatory synapses (Figure 8(b)), promoted by the Hebbian rule (Figure 8(c)) and the compensation theory (Figure 8(c)), and decrease of the inhibitory synapses (Figure 8(b)), due to compensation. Later, synaptic activity increases until the emergence of SICA. Changes in activity and connectivity during and subsequent to the SICA were similar to those observed for the increased activity induced by stimulation. Furthermore, it can be seen that the connectivity of all groups of connections, including group (I) representing the interconnections between the granule cells, increases to a maximum level.

When the simulation includes the transient excitatory GABAa, an activity increase (Figure 9(a)) after the activity induced by stimulation can be observed. During this period, it can be observed that (Figures 9(b) and 9(c)) (i) the connectivity between the granule cells increased simultaneously with the increase of the Hebbian rule rate (); (ii) connections sent by granule cells to other cell types were less intense; (iii) connectivity of the MC to the other cell types reduced due to the compensation theory (); (iv) inhibitory connectivity increased for all cell types, due to compensation. After this period, a progressive increase of the connectivity takes place between the GC and also the excitatory connectivity of the MC. In respect to the inhibitory synapses, the main effect is reduction of the inhibitory input on the MC (Figure 9(b)).

4. Discussion

The aim of the work presented was to investigate, by means of mathematical modelling and computational simulation, the interplay between changes in synaptic connectivity and the induction of SICA on neuronal networks that represents the circuitry of the DG of the rat hippocampus. The neuronal networks used described the cellular activity with dependence on the synaptic connectivity, according to McCulloch-Pitts formalism [15], and on the non-synaptic mechanisms typical of the DG region. The correspondence between the simulations and the experimental findings allows for the proposal of possible mechanisms involved. Assuming that the synaptic connectivity can be altered following the compensation theory [17, 18] and the Hebbian rules [24], the simulations performed suggest the following:(i)the combined action of these synaptogenesis rules contributes to the occurrence of abnormal mossy fiber sprouting and SICA sustaining; (ii)besides the abnormal mossy fiber sprouting, other modifications of connectivity, such as those that result in reduction of inhibitory cell activity, contribute to the occurrence of SICA; (iii)the transient excitatory GABAa, during the inflammatory period subsequent to the status epilepticus, acts, intensifying the Hebbian rules action and, consequently, strengthening the abnormal mossy fiber sprouting.

4.1. Network in Absence of Stimulation

According to Schaefers et al. [34], the DG exhibits high synaptic turnover rate and has a highly dynamic structure not restricted to the prenatal ontogenesis. These properties are also maintained after birth, during development, and can also take place in an adult brain. Mimicking these features in the simulations, the fluctuation of the activities induces the synaptogenesis processes. Without stimulation, the network remains in equilibrium for almost all combinations of the compensatory theory with the Hebbian rules. Only in a narrow range of lower values and higher do the SICA increase (Figure 4). These simulations suggest that if alterations on the brain tissue reinforce the role of compensation theory and depress the Hebbian rules’ action, the network will evolve in time favouring the emergence of epileptiform events, even in the absence of stimulation or another action to increase excitation. In this situation, two alterations in the connectivity can be qualified as the main reason for the SICA (Figure 5): (i) the decrease of the inhibition action on the mossy cells and (ii) the increase of the inhibition action on the interneurons.

The first alteration in the connectivity contributes to the emergence of the SICA, which in turn enhances the activation of the granule cells. This mechanism indicated by the simulations is consistent with the hypothesis formulated by Santhakumar et al. [35]. This hypothesis, named “irritable mossy cells hypothesis,” proposes that the increase of the mossy cell excitability is responsible for the generation and sustaining of epileptic seizures.

The increase of the inhibitory action on the interneurons, therefore the second connectivity alteration, results in granule and mossy cells disinhibition. This action is also favourable to the SICA emergence. This mechanism can be associated with the recurrent inhibition hypothesis, suggested by Sharma et al. [3], and also to explain the generation and sustaining of the epileptic seizures. In this hypothesis, the interneurons’ death or activation reduction is the main cause of the seizure disruption. The simulations present additional mechanisms that can produce the same effect. The increase of the interconnectivity between interneurons, aiming at compensating the network activity changes, can also contribute to reduce the inhibition efficacy, enhancing the granule cell excitation and therefore the seizures.

In the simulations, without stimulation, the SICAs were not due to the abnormal mossy fibers sprouting, since there was no significant increase in the interconnectivity of the granule cells (Figure 4—group (I)). This indicates that the abnormal mossy fiber sprouting, although often observed in experimental epilepsy [36], is not necessary for the seizure disruption. This result is consistent with the experimental findings from Buckmaster and Lew [37]. The epileptic rats treated with rapamycin, a suppressant of the axon sprouting, had no effect on seizure occurrence. It has been shown that amygdala, subiculum, and other temporal limbic structures could replace the mossy fiber effects. In fact, other brain areas show increases in recurrent excitation in models of TLE. The simulations offer an additional insight into this field, the network disconnected from other structures, and without mossy fiber sprouting was able to support synaptic reorganization with generation of spontaneous seizure by means of the Hebbian rules and compensation theory.

4.2. Network with Stimulation

When stimulated, the network simulation shows an increase in the neuronal activity enough to overcome the threshold for actuation of nonsynaptic mechanisms with subsequent triggering of SICA (Figure 2). The SICA duration was about 1.8 h and, associated with the experimental conditions, represents the status epilepticus. According to the simulations, the status epilepticus may be modulated by nonsynaptic mechanisms. In fact, experimental observations offer support to this speculation, since it was confirmed that nonsynaptic mechanisms are able to sustain epileptiform seizures in DG of rats’ hippocampus [14, 22, 3840].

During the period corresponding to the status epilepticus, the compensation theory causes excitatory reduction and inhibitory increase in the network connectivity. In consequence, immediately after the end of the period, the activity of the neural network was reduced, characterized as the latency period for the occurrence of SICA, which may be associated with the latency period observed for epilepsy induction with experimental models [7, 41]. During this period, following the compensation theory, new excitatory synapses are formed with the goal of reversing the activity reduction. Among these excitatory synapses included the interconnectivity between granule cells and the increase which characterizes the abnormal mossy fiber sprouting. The simulations suggest that the abnormal mossy fiber sprouting starts as a protective mechanism of the neural network, in order to compensate the increase of inhibitory synapses that occur during status epilepticus. This hypothesis reinforces the view that the axonal sprouting and the reactive synaptogenesis of the dentate gyrus serve as a repair mechanism, in order to restore normal function of the tissue [36]. However, an increase in the interconnectivity of the granule cells improves the performance of the Hebbian rules and the firing of a pre-synaptic cell is able to influence even more the firing of the postsynaptic cell connected. Consequently, the simulations suggest that the Hebbian rules, reinforced by the compensation theory, contribute to intensify the abnormal mossy fiber sprouting. This may happen along several days, during the latency period and even after SICA disruption. The intense interconnectivity between granule cells (Figures 5(b) and 7(b)) and the resulting increase of SICA occurrence (Figure 4(b)) take place for higher values (compensation theory) and lower, however not null, value (Hebbian rule), supporting the hypothesis that both theories (compensation and Hebbian) may be realized in the same network [19]. This dependence on and suggests that the interconnectivity increases between granule cells (abnormal mossy fiber sprouting) and the SICA happens only if the compensation theory forms inhibitory synapses and degenerates excitatory synapses more efficiently than the Hebbian rules increaseing the synaptic weight of the excitatory connections of the network. Additionally, the simulations show that the Hebbian rules are indispensable. When , the stimulation leads the network to a degeneration due to the action of the compensation theory.

4.3. Network with Stimulation and Transient Excitatory GABAa

The simulation of the network including the stimulation and also the excitatory effect of the GABAa synapses increased the occurrence of SICA and the range of combining values for and able to sustain them. The main changes observed in these simulations were the intense increase of the interconnectivity between granule cells (Figures 7 and 9). The simulations suggest that the excitatory GABAa synapses counteract the activity reduction caused by the increase of inhibitory connectivity along of the status epilepticus. This process enhanced granule cell firing, which, after the status epilepticus, was more interconnected due to the compensation theory. Thus, the granule cell promotes even more firing of action potentials between them. Therefore, the Hebbin rules were reinforced and the abnormal mossy fiber sprouting was intensified. These simulations support the conjecture that slow plasticity changes (Hebbian plasticity) after brain injury or status epilepticus may act on the eventual transition to the epilepsy state [42].

5. Conclusion

The simulations of the present paper allow speculation about the conjoint action of the compensation theory and the Hebbian rules during status epilepticus and the latency period promoting synaptic circuitry changes favourable to spontaneous seizure (Figure 10). According to the simulations, the network stimulation (electrical, chemical, injury, etc.) increases the neuronal activity leading to seizure mimicking the status epilepticus. During status epilepticus, seizures are strongly modulated by nonsynaptic mechanisms that sustain the depolarization even with the synaptic circuitry changes promoted by the compensation theory. The Hebbian rules take place, strengthening the current excitatory connections. When the nonsynaptic mechanisms of inhibition, in the case the Na/K pump, act to overcome the excitation mechanisms (e.g., current Na+ channels), the status epilepticus ends. From that moment, the excitability shows a dramatic reduction characterizing the beginning of the latency period. At the beginning of this period the compensation theory acts, reducing inhibition and increasing excitation. Contributing particularly to the excitation increase the simulations show the emergence of excitatory connections between granule cells, which indicates the onset of abnormal mossy fiber sprouting. During this period, the non-synaptic mechanisms do not contribute significantly to the excitability. In this sense, the abnormal mossy fiber sprouting arises as a protective mechanism, which tends to reverse the reduction of granule cell activities induced by the status epilepticus. The mossy fiber sprouting increases the connectivity between granule cells and reinforces the action of the Hebbian rules. In positive feedback, the Hebbian rule enhances the mossy fiber sprouting and helps to increase granule cell activity. Additionally, at the beginning of the latency period, the effect of the transient excitatory GABAa, promoted by the inflammatory process, contributes to reverse the neuronal activity reduction. This effect further enhances the performance of the Hebbian rules on the mossy fiber sprouting. As long as the abnormal mossy fiber sprouting intensifies, the excitation increases and, following the compensation theory, the network reacts to reverse this situation. When the abnormal mossy fiber sprouting is able to increase the granule cell activities overcoming the threshold of the nonsynaptic mechanisms actuation, the transition to the epilepsy state takes place.

The predictions proposed in this paper could be used to experimentally test the conjoint action of the compensation theory and Hebbian rule on the epileptogenesis induced by models of status epilepticus.

Acknowledgments

This work was supported by the Brazilian agencies Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG), Fundação CAPES, FAPESP/FAPEMIG, and INNT.