Skip to main content

Mechanisms of Intermittent State Transitions in a Coupled Heterogeneous Oscillator Model of Epilepsy

Abstract

We investigate the dynamic mechanisms underlying intermittent state transitions in a recently proposed neural mass model of epilepsy. A low dimensional model is constructed, which preserves two key features of the neural mass model, namely (i) coupling between oscillators and (ii) heterogeneous proximity of these oscillators to a bifurcation between distinct limit cycles. We demonstrate that state transitions due to intermittency occur in the abstract model. This suggests that there is a general bifurcation mechanism responsible for this behaviour and that this is independent of the precise form of the evolution equations. Such abstractions of neural mass models allow a deeper insight into underlying dynamic and physiological mechanisms, and also allow the more efficient exploration of large scale brain dynamics in disease.

1 Introduction

Epilepsy is a prevalent neurological disorder characterised by the repeated occurrence of pathological brain states known as seizures. Seizures are often accompanied by marked changes in electroencephalogram (EEG) dynamics. A pertinent example is the case of absence epilepsy in which short epochs of high amplitude, slow spike-wave rhythms spontaneously arise from a low amplitude background EEG (see, e.g. [1]). Although the reasons for these dynamic state changes are unknown, mathematical modelling can be used to help our understanding of the epileptic brain by highlighting possible underlying mechanisms. In addition, potential differences between the epileptic and non-epileptic brain can be revealed and seizure abatement or prevention strategies can be explored (see, e.g. [2]).

Physiologically inspired macroscopic brain models (for example neural mass and neural field models) have often been employed to investigate the causes of epilepsy and seizures [38]. In characterising the epileptic brain using these models, one can consider ways in which seizures can sporadically reoccur. In general, three methods have been proposed, namely bifurcations, bistability, and intermittency, with each of these potential dynamic mechanisms receiving both experimental and modelling support [4, 5, 7, 911]. In particular, intermittency has been postulated to underlie the generation of absence seizures in experimental animal models [10, 12]. A mechanistic neural mass model was recently shown to display intermittent dynamic transitions reminiscent of absence seizures [7, 13]. In the model of [7], intermittency arose due to interactions between heterogeneous populations of neurons. However, it remains unclear exactly which properties of this model imbue it with these important dynamics, and whether these properties are specific to the neural mass formalism. A greater understanding of these dynamics will highlight which features of interconnected regions of the brain are potentially responsible for seizure transitions.

Here, we approach this problem by extracting what we believe to be the fundamental dynamic features of the neural mass model, i.e. those features which, independent of the precise form of the equations, lead to the intermittent behaviour. We then construct a simple, low dimensional system preserving these features and demonstrate that this reduced system displays dynamics similar to those of the neural mass model. This makes a deeper understanding of the mechanisms of intermittent transitions possible, whilst it also shows that the proposed mechanism is dynamically robust. Furthermore, we exploit the reduced dimensionality of our new model to explore the effects of coupling in larger, heterogeneous systems. We propose this model abstraction strategy as a complement to the more detailed neural mass and neural network models for exploring the dynamic mechanisms of epilepsy.

2 Methods

2.1 Model

Our starting point is the work of [7], in which a model composed of connected neural masses was formulated as a representation of interacting populations of neurons in the cortex (e.g. interacting cortical columns). It was shown that this model displayed intermittency, autonomously switching between low amplitude, comparatively fast (i.e. “alpha band”) oscillations and high amplitude, slow spike-wave oscillations. These dynamics are reminiscent of electrographic recordings from patients with absence epilepsy. In the model, we identify the high amplitude oscillation (ghost of the limit cycle, i.e. “laminar” phase) with seizure EEG and the background EEG state is represented by the global re-injection of the model (i.e. the “turbulent” phase). Both intrinsic parameter values of the neural masses (the nodes of the network) and the connectivity between nodes were important for the observed intermittent dynamics in the model, thus pathological rhythm generation was placed into the context of interacting, heterogeneous regions of tissue.

Here, we explore the hypothesis that the intermittent dynamics of the neural mass model are due to general dynamic properties, namely (i) the intrinsic dynamic repertoire (bifurcations) of isolated nodes and (ii) the interaction between nodes when they are coupled, and when nodes are heterogeneous. We therefore proceed by engineering a simple deterministic system of ordinary differential equations, which incorporates the features we believe to be important for the observed neural mass model intermittency. We first consider the dynamics of individual nodes, which in the neural mass model were close to a bifurcation between qualitatively different oscillations, with a region of bistability between the two. We therefore begin by constructing a simple two-dimensional non-linear system with these features (see details below). We then introduce global coupling between these nodes and examine the dynamics of the coupled system to test for the emergence of intermittency.

Our simplified model of the neural mass is constructed so that the amplitude of oscillations (r) is controlled by a bifurcation parameter, μ. Choosing a cubic form for the dependence of R= r 2 on μ and then rotating the system using an angle variable, θ, gives rise to the bifurcation structure required, as shown in Fig. 1. The slowing of oscillations with increasing amplitude is obtained by making θ ˙ a decreasing function of R. In polar coordinates, the model for a single compartment is thus:

r ˙ = r ( μ a r 2 + b r 4 c r 6 ) θ ˙ = ω d r 2
(1)
Fig. 1
figure 1

Model bifurcations. Bifurcations over μ in the single compartment, 2 variable model. Top: limit cycle amplitude. Bottom: angular velocity

To couple N different compartments (specified by subscript i), we use simple additive coupling, so that in Cartesian coordinates the equations are:

x i ˙ = y i ( ω d r i 2 ) + x i ( μ a r i 2 + b r i 4 c r i 6 ) + β A x y ˙ i = x i ( ω d r i 2 ) + y i ( μ a r i 2 + b r i 4 c r i 6 )
(2)
i=1,,N

, where

r i 2 = x i 2 + y i 2

The term A x provides the additive coupling, where A is the adjacency matrix for a generic network and x=( x 1 ,, x N ). As in [7], self-coupling has been included in the intrinsic node equations and, therefore, does not appear explicitly in A (i.e. A i i =0 i).

For convenience, we fix the parameters a=2, b=3/2, and c=1/3. These determine the position of the region of bistability between different oscillation types of the uncoupled system. To see this (and to make sense of future diagrams), we give a brief description of the dynamics of (1). Since r ˙ is independent of θ, stationary points and periodic orbits lie on contours of r ˙ =0, i.e. r=0, the fixed point, which is stable if μ<0 and unstable if μ>0, and the solutions of

0=μa r 2 +b r 4 c r 6

which will correspond to a degenerate circle of fixed points if r 2 =ω/d is a solution (a special case we ignore), and a periodic orbit otherwise. It is easier to analyze solutions by setting r 2 =R and looking for positive solutions to

μ=aRb R 2 +c R 3

This explains why we have chosen the parameters above: let F(R)=aRb R 2 +c R 3 , then

F (R)=a2bR+3c R 2

and so turning points of this curve (seen as a curve of solutions in the (μ,R)-plane) occur at F (R)=0 or, for the parameters a=2, b=3/2, and c=1/3,

0=23R+ R 2 =(R2)(R1)

In other words, turning points, which correspond to saddle-node bifurcations of periodic orbits occur at R=2 and R=1 (r= 2 , r=1). This structure is shown in Fig. 1. If R=2, then μ= 2 3 , and if R=1, then μ= 5 6 . A Hopf bifurcation at μ=0 creates a stable low amplitude limit cycle. As μ increases through μ= 2 3 , a second (large amplitude) limit cycle is created together with a large amplitude unstable periodic orbit in a saddle-node bifurcation. Following this bifurcation, there is a region of bistability until μ= 5 6 when the stable small amplitude limit cycle is destroyed in the second saddle-node bifurcation. It is this region of bistability that we aim to use in the coupled compartment models to generate intermittency by choosing heterogeneous parameters a little below μ= 2 3 , so that although each individual compartment lies in a region with a unique stable small amplitude limit cycle, their coupling can cause temporary (but repeated) excitation into the high amplitude state. This is essentially the mechanism proposed in [7].

2.2 Connectivity

In this study, we explore three different connectivity schemes. In order to demonstrate the equivalence of our model dynamics with those of [7], we initially study an all-to-all coupled system of three compartments (see [13]). For simulations of spatially extended systems, we work in two dimensions with periodic boundaries and symmetric nearest neighbour or distance dependent connectivity. Connectivity weights are scaled by a connectivity parameter, β. For systems with distance dependent connectivity, an exponential fall off of connectivity strength is used as follows:

a i j = e α r i r j max i j ( e α r i r j )
(3)

where a i j is the entry of A connecting nodes i and j and r i and r j are the 2-d coordinates of the location of compartments i and j, respectively. The denominator scales the connectivity such that the maximum value of A is 1, and this is subsequently scaled in Eq. (2) by β.

3 Results

3.1 Type 1 Intermittency in 3 Coupled Compartments

It has been suggested that the sporadic nature of epileptic episodes in absence seizures is due to dynamic intermittency [10, 12]. Our recent work demonstrated that intermittency in heterogeneous neural mass models can arise due to a type 1 route, i.e. a proximal tangent bifurcation of a limit cycle [13, 14]. These dynamics arose in a system composed of 3 all-to-all coupled compartments. We therefore investigate whether these dynamics can be recreated in our abstract model, solely based on the generic dynamic features extracted, as described in Sect. 2.

In Fig. 2, we demonstrate that intermittency does emerge from our abstract model formulation. Figure 2 shows a long time series with parameter heterogeneity fixed to μ 1 =0.2, μ 2 =0.3, and μ 3 =0.6 and all-to-all connectivity scaled by β=3/2. Each compartment makes intermittent deviations from the low amplitude oscillatory state into a high amplitude, slow oscillation. A close up of these model dynamics can be seen in Fig. 3.

Fig. 2
figure 2

Intermittency in the abstract model. Intermittent switching between low amplitude and high amplitude oscillations in variable x 3 of the model (Eq. (2)) with μ 1 =0.2, μ 2 =0.3 and μ 3 =0.6 and β=3/2

Fig. 3
figure 3

Intermittent dynamics. Left: close up of a state transition in the x-variable of each of 3 compartments of the model demonstrated in Fig. 1. Right: phase portrait showing the trajectory of compartment 3 in xy space

For fixed μ i we explored the effect of changing coupling strength, β, which is shown in the bifurcation diagram of Fig. 4. For high β, the system evolves with synchronous, stable limit cycle oscillations. As coupling strength is decreased these oscillations give way to an intermittent regime, with the system exploring both the high and low amplitude oscillations of the single compartment model. We note that the way in which these model dynamics vary in relation to coupling strength is very similar to the full neural mass model (see [13]). In order to explore the dynamic mechanisms leading to this intermittent window, a long simulation was performed close to the onset of intermittency. Analysis of this system revealed a U-shaped distribution of high amplitude state durations as well as a proximate tangent bifurcation of the phase locked limit cycle (Fig. 5). Thus, our abstract model also follows a type 1 route into intermittency (Chap. 5, pp. 68–69 in [15]), [13].

Fig. 4
figure 4

Effect of coupling strength. Scan of the dynamics of the 3 compartment model over changes in the coupling strength, β. The distinct maxima of x 3 are shown, with black and red dots denoting forward and backward scans, respectively. All other parameters are as in Fig. 2

Fig. 5
figure 5

Type 1 intermittency. The main figure shows a first return map of variable y 1 of the system close to the transition into intermittency (β=1.897) and demonstrates the presence of a tangent bifurcation. The inset shows the distribution of lengths of the high amplitude states

3.2 Dynamics of Larger Systems

An advantage of the reduced dimensionality of our new model is the added computational efficiency for simulating large systems. Thus, the role of connectivity and heterogeneity in intermittent state changes can be explored more easily. In this study, we explore the dynamics of a 9×9 sheet of compartments under different connectivity schemes. We begin with nearest neighbour coupling, which is a simplistic approximation to the predominance of connectivity between proximal tissue in, for example human cortex.

The model can be used to investigate the interplay between network connectivity and the susceptibility of nodes to abnormal states (which is referred to as “epileptogenicity”) [8]. Specifically, higher values of the parameter μ can be thought of as being more epileptogenic, since in a single compartment this would render a node closer to the higher amplitude rhythm. As a demonstration of this line of enquiry, a square region of compartments with μ=0.6 is placed at the centre of a nearest neighbour coupled system. This central square is surrounded by “normal” compartments, which are further from the bifurcation and have μ drawn from a normal distribution with mean 0.3 and variance 0.05. Figure 6 shows the distribution of μ in the system and a visualisation of the connectivity matrix, A.

Fig. 6
figure 6

Schematic of a system with nearest neighbour coupling. Layout of the 9×9 system of coupled compartments. Left: colour coded image of the distribution of μ. The central 9 square s have μ=0.6, whereas the rest of the system has μ drawn from a normal distribution with mean 0.3 and variance 0.05. Numbers indicate compartment labels for reference. Right: demonstration of connectivity for the central compartment. Black indicates the presence of connections

As in the smaller system, A was scaled by β in order to explore changes in connectivity strength. We found that β can be tuned such that compartments of the system undergo intermittent deviations into the high amplitude oscillatory state. An example of the dynamics of the central compartments is given in Fig. 7. It can be seen that different combinations of these compartments undergo periods of high amplitude oscillations. Compartments 1 and 9, at the periphery of the central square are less easily perturbed into the high amplitude state than the other compartments. In particular, compartment 9 is coupled to two compartments with very low values of μ (white squares in Fig. 6) and deviates less often into high amplitude rhythms. Thus, it is demonstrated in this system that both the epileptogenicity of nodes and their connected networks is important for state transitions.

Fig. 7
figure 7

Dynamics of a system with nearest neighbour coupling. Time series of the dynamics of the 9 central compartments of the system described in Fig. 6. Compartment numbering is given in Fig. 6. β=2.75

In addition to nearest neighbour coupling, a system with distance dependent connectivity was also explored (see Sect. 2), which more closely approximates short range interactions in the brain. In this system, a graded epileptogenic region was also incorporated by allowing μ to decay exponentially from its value of 0.6 in the central compartment to 0.3 at the periphery. The connectivity matrix (βA) and the distribution of μ for this system can be seen in Fig. 8.

Fig. 8
figure 8

Schematic of a system with distance dependent coupling. Layout of the 9×9 system of coupled compartments with distance dependent coupling and exponential decay in μ. Left: colour coded image of the distribution of μ. The central compartment has μ=0.6, and the value of μ decays exponentially to 0.3 at the periphery of the system. Numbers indicate compartment labels for reference. Right: colour coded connectivity values for the central compartment

This system also displayed intermittent high amplitude bursts, as can be seen in Fig. 9. In this case, several bursts were more coherent across the central compartments, even though in this system μ varied between compartments with distance from the centre.

Fig. 9
figure 9

Dynamics of a system with distance dependent coupling. Time series of the dynamics of the 9 central compartments of the system described in Fig. 8. Compartment numbering is also provided in Fig. 8. β=5

4 Discussion

In this paper, we have shown that a network of coupled two-dimensional equations (compartments) can display intermittent behaviour with transitions between low amplitude and high amplitude oscillations. The equations were chosen so that in the absence of coupling the parameters of each compartment were close to a saddle-node bifurcation of periodic orbits creating the large amplitude oscillations, but the different compartments had different parameter values (heterogeneity). This behaviour has also been observed in a more detailed, higher dimensional model of interacting neural masses used as a model of absence epilepsy, so our results suggest that there is an underlying bifurcation mechanism leading to this type of dynamics. This, in turn, shows that the intermittency of the detailed model is robust—it does not depend on the precise terms and parameters of the model. This is important since the model is, by its nature, an approximation. Thus, we have provided strong evidence that these dynamics are due to coupling between heterogeneous systems that posses a region of bistability flanked by saddle-node bifurcations of qualitatively different limit cycles. Furthermore, we confirmed that the bifurcation into the intermittent regime in our reduced model follows the type 1 route [14], which is equivalent to the dynamics of the neural mass model [13]. In addition, we demonstrated the potential of the reduced model to simulate large networks, and hence explore the interplay between network connectivity and spatial distribution of heterogeneities in epileptogenic networks [7, 8, 1619].

Large scale dynamic models are important for understanding both function and dysfunction in the brain [1922]. However, as models become larger and more complex, they become difficult to systematically explore, and extraction of underlying principles governing their dynamic repertoire is not trivial. Epilepsy is a pertinent example, where we are challenged to relate multi-scale mechanisms and spatially distributed abnormalities (epileptogenic zones [23] and networks [24]) to the production of pathological phenotypes, which include abnormal electrographic dynamics. Here, we have demonstrated that principles underlying such dynamics can be uncovered by pursuing sequential levels of abstraction. In [7], we formulated a large dimensional neural mass model with the emergent phenomenon of intermittency resembling epileptic electrographic activity. Subsequently in [13], a reduced dimensional neural mass model (i.e. composed of fewer compartments) was shown to preserve this feature and allowed for the categorisation of the route into intermittency at some parameter values as type 1 [14]. In the current study, we have added a further level of abstraction, retaining only certain dynamic features of the original model.

In terms of brain dynamics, the theoretical prediction of this work is that oscillations in different parts of epileptic brain networks have different degrees of “excitability” due to their proximity to a bifurcation into a different oscillatory dynamic regime. When these regions communicate, e.g. via synaptic connectivity, an intermittent dynamic regime can occur, with spontaneous episodes of qualitatively different dynamics. We have therefore demonstrated a high level theoretical link between brain oscillations [25], differential excitability in brain networks (e.g. [26]) and abnormal rhythm generation, which are crucial concepts for epilepsy. An important next step will be to study models such as the one presented here in order to explore more precisely the relationship between network topology and distributed heterogeneities, which will provide a better understanding of epileptogenic networks and allow to predict the nature of inter-connected normal and abnormal regions of tissue in the epileptic brain [8, 16, 19]. An advantage of the current framework is that such questions can be explored with improved computational efficiency as compared to neural mass models.

To confirm the use of our model in this direction, we demonstrated the preservation of intermittency in systems incorporating spatially structured connectivity. Interestingly, in these systems, we did not immediately uncover a dynamic regime in which the majority of nodes displayed concomitant and coordinated switching into abnormal dynamics, as would be expected for absence epilepsy [1]. Our experience simulating spatially structured neural mass models, together with considerations of the anatomy of the brain, suggests that such a regime could require the addition of long range connections. The exploration of this hypothesis will further uncover the link between the topology of large scale brain networks and the propagation of epileptiform activity and will be an interesting avenue for future study.

An alternative dynamic regime proposed to underly switching between seizure and non-seizure states in absence seizures is noise-driven bistability [4, 2730]. Similarly to our current work, this dynamic regime has been investigated both in mechanistic neural mass models [4] and in reduced forms preserving bistability between a stable steady state and a stable limit cycle [29, 30]. In the current study, we employ similar equations, though the resulting bifurcation structure is different as we focus upon transitions between qualitatively different oscillations (following [7]), rather than between a fixed point and a limit cycle. Furthermore, [30] and [29] explored the effect of qualitative changes in network structure on the resulting dynamics. In addition to this work in bistable systems, future such explorations in our low dimensional framework for intermittency can provide insight into the role of network connectivity on seizure occurrence. In addition, recent work has explored the possibility that critical transitions are responsible for seizure termination in focal onset seizures with secondary generalisation. [11] showed that certain features of these transitions, such as slowing and flickering [31] can be seen in recordings from patients with these seizure types, and slowing of absence seizure rhythms toward seizure termination does occur. In our intermittent system, the dynamics of seizure termination are governed by the nature of the vector field where the system leaves the previously stable limit cycle. Whether such hallmarks can be seen in intermittent models is yet to be determined. Although in the current study we focused on robustness of the mechanisms of intermittent transitions, in future studies, it will be interesting to explore other details of the waveforms of seizure EEG.

A similar model approximation approach was taken recently by [32], who sought to uncover the dynamic principles underlying multi-stability and scale invariant fluctuations in the human alpha rhythm. Similarly to the current study, their starting point was a biophysically inspired neural mass model which offered a novel explanation for the observed data [33]. Studying a normal form for the dynamics of interest, the requirements were shown to be a sub-critical Hopf bifurcation and the incorporation of multiplicative noise. Other abstract modelling studies have demonstrated minimal neural mass models underlying the generation of different epileptic electroencephalographic waveforms [34], as well as the complex role of networks underlying focal and generalised epileptic seizures [17]. The approach of modelling at several levels of abstraction will be a key strategy in advancing our understanding of the brain.

In summary, we have uncovered a dynamic mechanism responsible for spontaneous state transitions in a neural mass model of epilepsy. Future explorations of this model and further developments and applications of the general methods employed herein will help to advance our understanding of large scale brain dynamics in health and disease. Since our model is framework independent, it might also be used to explore rhythm generation in other systems, for example oscillating biochemical networks [35].

References

  1. Sadleir LG, Farrell K, Smith S, Connolly MB, Scheffer IE: Electroclinical features of absence seizures in childhood absence epilepsy. Neurology 2006, 67(3):413–418. 10.1212/01.wnl.0000228257.60184.82

    Article  Google Scholar 

  2. Lopes da Silva F, Blanes W, Kalitzin SN, Parra J, Suffczynski P, Velis DN: Epilepsies as dynamical diseases of brain systems: basic models of the transition between normal and epileptic activity. Epilepsia 2003, 44: 72–83. (Suppl 12) (Suppl 12) 10.1111/j.0013-9580.2003.12005.x

    Article  Google Scholar 

  3. Wendling F, Bellanger JJ, Bartolomei F, Chauvel P: Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals. Biol Cybern 2000, 83(4):367–378. 10.1007/s004220000160

    Article  Google Scholar 

  4. Suffczynski P, Kalitzin S, Lopes da Silva F: Dynamics of non-convulsive epileptic phenomena modeled by a bistable neuronal network. Neuroscience 2004, 126(2):467–484. 10.1016/j.neuroscience.2004.03.014

    Article  Google Scholar 

  5. Breakspear M, Roberts JA, Terry JR, Rodrigues S, Mahant N, Robinson PA: A unifying explanation of primary generalized seizures through nonlinear brain modeling and bifurcation analysis. Cereb Cortex 2006, 16: 1296–1313.

    Article  Google Scholar 

  6. Marten F, Rodrigues S, Benjamin O, Richardson MP, Terry J: Onset of polyspike complexes in a mean-field model of human electroencephalography and its application to absence epilepsy. Philos Trans R Soc A, Math Phys Eng Sci 2009, 367(1891):1145–1161. 10.1098/rsta.2008.0255

    Article  MathSciNet  MATH  Google Scholar 

  7. Goodfellow M, Schindler K, Baier G: Intermittent spike-wave dynamics in a heterogeneous, spatially extended neural mass model. NeuroImage 2011, 55(3):920–932. 10.1016/j.neuroimage.2010.12.074

    Article  Google Scholar 

  8. Goodfellow M, Schindler K, Baier G: Self-organised transients in a neural mass model of epileptogenic tissue dynamics. NeuroImage 2012, 59(3):2644–2660. 10.1016/j.neuroimage.2011.08.060

    Article  Google Scholar 

  9. Suffczynski P, Lopes da Silva FH, Parra J, Velis DN, Bouwman BM, van Rijn CM, van Hese P, Boon P, Khosravani H, Derchansky M, Carlen P, Kalitzin S: Dynamics of epileptic phenomena determined from statistics of ictal transitions. IEEE Trans Biomed Eng 2006, 53(3):524–532. 10.1109/TBME.2005.869800

    Article  Google Scholar 

  10. Hramov A, Koronovskii AA, Midzyanovskaya IS, Sitnikova E, van Rijn CM: On-off intermittency in time series of spontaneous paroxysmal activity in rats with genetic absence epilepsy. Chaos 2006., 16(4): Article ID 043111 Article ID 043111

  11. Kramer M, Truccolo W, Eden UT, Lepage KQ, Hochberg LR, Eskandar EN, Madsen JR, Lee JW, Maheshwari A, Halgren E, Chu CJ, Cash SS: Human seizures self-terminate across spatial scales via a critical transition. Proc Natl Acad Sci USA 2012, 109(51):21116–21121. 10.1073/pnas.1210047110

    Article  Google Scholar 

  12. Sitnikova E, Hramov AE, Grubov VV, Ovchinnkov AA, Koronovsky AA: On-off intermittency of thalamo-cortical oscillations in the electroencephalogram of rats with genetic predisposition to absence epilepsy. Brain Res 2012, 1436: 147–156.

    Article  Google Scholar 

  13. Goodfellow M, Rummel C, Garry D, Baier G, Schindler K, Glendinning P: State transitions in a model of intermittent seizure dynamics. Preprint 2012.74, MIMS; 2012. Goodfellow M, Rummel C, Garry D, Baier G, Schindler K, Glendinning P: State transitions in a model of intermittent seizure dynamics. Preprint 2012.74, MIMS; 2012.

  14. Pomeau Y, Manneville P: Intermittent transition to turbulence in dissipative dynamical systems. Commun Math Phys 1980, 74: 189–197. 10.1007/BF01197757

    Article  MathSciNet  Google Scholar 

  15. Schuster HG, Just W: Deterministic Chaos: An introduction. 4th edition. New York, Wiley; 2005.

    Book  Google Scholar 

  16. Goodfellow M, Taylor PN, Wang Y, Garry DJ, Baier G: Modelling the role of tissue heterogeneity in epileptic rhythms. Eur J Neurosci 2012, 36(2):2178–2187. 10.1111/j.1460-9568.2012.08093.x

    Article  Google Scholar 

  17. Terry JR, Benjamin O, Richardson MP: Seizure generation: the role of nodes and networks. Epilepsia 2012, 53(9):e166-e169. 10.1111/j.1528-1167.2012.03560.x

    Article  Google Scholar 

  18. Kramer MA, Cash SS: Epilepsy as a disorder of cortical network organization. Neuroscientist 2012, 18(4):360–372. 10.1177/1073858411422754

    Article  Google Scholar 

  19. Richardson MP: Large scale brain models of epilepsy: dynamics meets connectomics. J Neurol Neurosurg Psychiatry 2012, 83(12):1238–1248. 10.1136/jnnp-2011-301944

    Article  Google Scholar 

  20. Deco G, Jirsa VK, Robinson PA, Breakspear M, Friston K: The dynamic brain: from spiking neurons to neural masses and cortical fields. PLoS Comput Biol 2008., 4(8): Article ID e1000092 Article ID e1000092

  21. Coombes S: Large-scale neural dynamics: simple and complex. NeuroImage 2010, 52(3):731–739. 10.1016/j.neuroimage.2010.01.045

    Article  Google Scholar 

  22. Jirsa VK, Sporns O, Breakspear M, Deco G, McIntosh AR: Towards the virtual brain: network modeling of the intact and the damaged brain. Arch Ital Biol 2010, 148(3):189–205.

    Google Scholar 

  23. Rosenow F, Lüders H: Presurgical evaluation of epilepsy. Brain 2001, 124: 1683–1700. 10.1093/brain/124.9.1683

    Article  Google Scholar 

  24. Spencer SS: Neural networks in human epilepsy: evidence of and implications for treatment. Epilepsia 2002, 43(3):219–227. 10.1046/j.1528-1157.2002.26901.x

    Article  Google Scholar 

  25. Thut G, Miniussi C, Gross J: The functional importance of rhythmic activity in the brain. Curr Biol 2012, 22(16):R658-R663. 10.1016/j.cub.2012.06.061

    Article  Google Scholar 

  26. Valentín A, Alarcón G, García-Seoane JJ, Lacruz ME, Nayak SD, Honavar M, Selway RP, Binnie CD, Polkey CE: Single-pulse electrical stimulation identifies epileptogenic frontal cortex in the human brain. Neurology 2005, 65(3):426–435. 10.1212/01.wnl.0000171340.73078.c1

    Article  Google Scholar 

  27. Kalitzin SN, Velis DN, da Silva FHL: Stimulation-based anticipation and control of state transitions in the epileptic brain. Epilepsy Behav 2010, 17(3):310–323. 10.1016/j.yebeh.2009.12.023

    Article  Google Scholar 

  28. Kalitzin S, Koppert M, Petkov G, Velis D, da Silva FL: Computational model prospective on the observation of proictal states in epileptic neuronal systems. Epilepsy Behav 2011, 22: S102-S109. (Suppl 1) (Suppl 1)

    Article  Google Scholar 

  29. Benjamin O, Fitzgerald TH, Ashwin P, Tsaneva-Atanasova K, Chowdhury F, Richardson MP, Terry JR: A phenomenological model of seizure initiation suggests network structure may explain seizure frequency in idiopathic generalised epilepsy. J Math Neurosci 2012. 10.1186/2190-8567-2-1

    Google Scholar 

  30. Koppert MMJ, Kalitzin S, Silva FL, Viergever MA, Model A: Connectivity and phase coherence in neural network models of interconnected Z 4-bistable units. Engineering in Medicine and Biology Society 2012 Annual International Conference of the IEEE (EMBC) 2012, 5458–5461.

    Chapter  Google Scholar 

  31. Scheffer M, Bascompte J, Brock W, Brovkin V, Carpenter SR, Dakos V, Held H, van Nes EH, Rietkerk M, Sugihara G: Early-warning signals for critical transitions. Nature 2009, 461(7260):53–59. 10.1038/nature08227

    Article  Google Scholar 

  32. Freyer F, Roberts J, Ritter P, Breakspear M: A canonical model of multistability and scale-invariance in biological systems. PLoS Comput Biol 2012., 8(8): Article ID e1002634

  33. Freyer F, Roberts J, Becker R, Robinson P, Ritter P, Breakspear M: Biophysical mechanisms of multistability in resting-state cortical rhythms. J Neurosci 2011, 31(17):6353–6361. 10.1523/JNEUROSCI.6693-10.2011

    Article  Google Scholar 

  34. Wang Y, Goodfellow M, Taylor PN, Baier G: Phase space approach for modeling of epileptic dynamics. Phys Rev E 2012, 85(6):1–11.

    Google Scholar 

  35. Novák B, Tyson JJ: Design principles of biochemical oscillators. Nat Rev Mol Cell Biol 2008, 9(12):981–991. 10.1038/nrm2530

    Article  Google Scholar 

Download references

Acknowledgements

MG acknowledges funding from the EPSRC through a postdoctoral prize fellowship.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Marc Goodfellow.

Additional information

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

MG and PG designed the study, performed analysis and wrote the paper. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Goodfellow, M., Glendinning, P. Mechanisms of Intermittent State Transitions in a Coupled Heterogeneous Oscillator Model of Epilepsy. J. Math. Neurosc. 3, 17 (2013). https://doi.org/10.1186/2190-8567-3-17

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2190-8567-3-17

Keywords