Characterizing Deep Brain Stimulation effects in computationally efficient neural network models
© Latteri et al; licensee BioMed Central Ltd. 2011
Received: 8 February 2011
Accepted: 15 April 2011
Published: 15 April 2011
Recent studies on the medical treatment of Parkinson's disease (PD) led to the introduction of the so called Deep Brain Stimulation (DBS) technique. This particular therapy allows to contrast actively the pathological activity of various Deep Brain structures, responsible for the well known PD symptoms. This technique, frequently joined to dopaminergic drugs administration, replaces the surgical interventions implemented to contrast the activity of specific brain nuclei, called Basal Ganglia (BG). This clinical protocol gave the possibility to analyse and inspect signals measured from the electrodes implanted into the deep brain regions. The analysis of these signals led to the possibility to study the PD as a specific case of dynamical synchronization in biological neural networks, with the advantage to apply the theoretical analysis developed in such scientific field to find efficient treatments to face with this important disease. Experimental results in fact show that the PD neurological diseases are characterized by a pathological signal synchronization in BG. Parkinsonian tremor, for example, is ascribed to be caused by neuron populations of the Thalamic and Striatal structures that undergo an abnormal synchronization. On the contrary, in normal conditions, the activity of the same neuron populations do not appear to be correlated and synchronized.
To study in details the effect of the stimulation signal on a pathological neural medium, efficient models of these neural structures were built, which are able to show, without any external input, the intrinsic properties of a pathological neural tissue, mimicking the BG synchronized dynamics.
We start considering a model already introduced in the literature to investigate the effects of electrical stimulation on pathologically synchronized clusters of neurons. This model used Morris Lecar type neurons. This neuron model, although having a high level of biological plausibility, requires a large computational effort to simulate large scale networks. For this reason we considered a reduced order model, the Izhikevich one, which is computationally much lighter. The comparison between neural lattices built using both neuron models provided comparable results, both without traditional stimulation and in presence of all the stimulation protocols. This was a first result toward the study and simulation of the large scale neural networks involved in pathological dynamics.
Using the reduced order model an inspection on the activity of two neural lattices was also carried out at the aim to analyze how the stimulation in one area could affect the dynamics in another area, like the usual medical treatment protocols require.
The study of population dynamics that was carried out allowed us to investigate, through simulations, the positive effects of the stimulation signals in terms of desynchronization of the neural dynamics.
The results obtained constitute a significant added value to the analysis of synchronization and desynchronization effects due to neural stimulation. This work gives the opportunity to more efficiently study the effect of stimulation in large scale yet computationally efficient neural networks. Results were compared both with the other mathematical models, using Morris Lecar and Izhikevich neurons, and with simulated Local Field Potentials (LFP).
PD is a degenerative disorder of the central nervous system that often impairs motor skills, speech, and other functions. It is characterized by muscle rigidity, tremor, a slowing of physical movements (bradykinesia) and, in extreme cases, a complete loss of physical movement (akinesia). The primary symptoms are the results of decreased stimulation of the motor cortex by the basal ganglia and other brain stem structures, traditionally considered as a consequence of the insufficient formation and action of dopamine, which is produced in the dopaminergic neurons of the Substantia Nigra reticulata (SNr). Other symptoms may include high level cognitive dysfunction and subtle language problems, postural instability and gait disturbances. In some cases, it would be inaccurate to say that the cause is "unknown", because a small proportion is caused by genetic mutations. It is possible for a patient to be initially diagnosed with PD but then to develop additional features, requiring revision of the diagnosis.
At present, there is no cure for PD, but medications or surgery can provide relief from the symptoms. The most widely used form of treatment is L-dopa in various forms. However, only 1-5% of L-DOPA enters the dopaminergic neurons. The remaining L-DOPA is often metabolized to dopamine elsewhere, causing a wide variety of side effects. Due to feedback inhibition, L-dopa results in a reduction in the endogenous formation of L-dopa, and so eventually becomes counterproductive. In the 1990s the surgical ablation has been used to treat PD. Ablative brain surgery is the surgical lesion by burning or freezing or with chemical substances of brain tissue to treat neurological or psychological disorders. The thalamus was a potential target for treating tremor, especially in the Ventral Intermediate Nucleus (VIM) and Centrum Medianum-Parafascicular (CM-Pf) nuclei. The lesions caused by this type of surgery are irreversible, so generally DBS surgery is considered preferable to lesion because it has the same effect and is adjustable and reversible [4–6]. Treating PD with surgery was once a common practice, but after the discovery of levodopa, surgery was restricted to only a few cases. DBS is a surgical treatment involving the implantation of a medical device called a brain pacemaker, which sends electrical impulses to specific parts of the brain. DBS in selected brain regions has provided remarkable therapeutic benefits for otherwise treatment-resistant movements and affective disorders such as chronic pain, PD or essential tremor and dystonia. Despite the long history of DBS, , its underlying principles and mechanisms are still unclear. DBS directly changes brain activity in a controlled manner, its effects are reversible (unlike those of lesioning techniques). The Food and Drug Administration (FDA) approved DBS as a treatment for essential tremor in 1997, for PD in 2002 , and dystonia in 2003 . DBS leads are placed in the brain according to the type of symptoms to be addressed. For dystonia and symptoms associated with PD (rigidity, bradykinesia/akinesia and tremor), the lead may be placed in either the Globus Pallidus or Subthalamic Nucleus . The right side of the brain is stimulated to address symptoms on the left side of the body and vice versa. DBS does not cure PD, but it can help to manage some of its symptoms and subsequently to improve the patient's quality of life . Presently, the procedure is used only for patients whose symptoms cannot be adequately controlled with medications, or whose medications have severe side effects . Its direct effect on the physiology of brain cells and/or neurotransmitters is currently debated, but it is apparent that sending high frequency electrical impulses into specific areas of the brain can mitigate symptoms  and/or directly decrease the side effects induced by Parkinsonian medications, , since it allows a sometimes huge decrease in medications, making the medication regime more tolerable. More recent neurophysiological data suggest that the DBS can modify also the connections among the cellular networks, giving rise to a holistic interpretation of DBS action these aspects that should be also considered in the future .
There are a few sites in the brain that can be targeted to achieve different results, so each patient must be assessed individually, and the particular site (or concurrent sites) to be stimulated is chosen based on the particular needs. Traditionally, the two most common sites are the Subthalamic Nucleus (STN) and the Globus Pallidus internus (GPi), but other sites, such as the CM-Pf  and Nucleus Tegmenti Peduncolopontini (PPTg) have recently shown important benefits for PD treatment . A tailored DBS with the targeting of multiple nuclei was proposed to obtain the high clinical results in complex parkinsonian syndromes . The main objective of stimulation is to re-establish desynchronization via a pulse train, whose parameters are selected by the Neurosurgeon at the main aim to decrease the disease symptoms. The patient with PD has a high synchronization into band β (10÷30 Hz) of the activities STN and GPi (brakinesia., akinesia., etc.), whereas the patient with L-Dopa into same band has high desynchronization, with an improvement of movements . These results confirm what is expected from the Gate Control Theory , which states that the synchronization of neuronal activity obstructs information flow in brain structures, whereas, the desynchronization allows information flow.
A increase of the neuronal synchronization causes a increase of power spectral density (PSD) of the Local Field Potential (LFP) taken from STN and GPi, whereas, the neuronal desynchronization causes a decrease of the PSD.
Such a desynchronization in the β band takes place also when the patient is subjected to DBS, sending a pulse train with frequency larger than 70 Hz . The DBS of the STN and GPi, involves a synchronization of neurons in the high frequencies and a desychronization in the β band. These stimuli have the common goal of reducing the synchronized activity of specific clusters of neurons, in order to desynchronize the target population, re-establishing a normal physiological activity in a highly synchronized population.
The first step toward the understanding of the dynamic effect of the DBS on the brain nuclei, is to design a model of these areas and simulate the effects caused by the stimulation. At this aim, after a deep scanning of the state of the art, in , a model of the neuronal BG was derived using a neural network made up of Morris-Lecar neurons, arranged in a specific population and showing, without any external input, basic oscillations mimicking the Parkinsonian tremor. As in , the positive effect of an electrical stimulation induced by implanted electrodes was demonstrated in such a population model. Therefore this model was initially taken into consideration as a reference. Then a fundamental issue, when the need of simulating large scale populations arises, lies in the fact that the number of floating point operations needed to simulate a time unit for a single neuron is relevant for the time needed to appreciate the results of the whole population. This led us to analyze other neuron models that could reveal the same characteristics (basically autonomous bursting oscillations) as the Morris Lecar model, but being at the same time much less demanding as regards the computational burden. Therefore, after analysing, in the first part of this paper in the method section, the effect produced by a population of Morris Lecar units, the same conditions were reproduced using a population consisting of Izhikevich neurons , where parameters were selected so as to show a resonant-like bursting characteristics, similar to the Morris Lecar model, but with a much less computational power demand. Of course, all the other characteristics, like the synaptic dynamics and all the other relevant parameters, were left unchanged.
In the second part of this paper, this new reduced order model gave us the possibility to increase, where needed, the population size using roughly the same computation time. The positioning of the electrodes within the population so built revealed the expected effects i.e. a desynchronization within the neurons, which led to a spreading of the power spectral density over a wide frequency range. This effect even increases in the presence of multiple stimulating sites within the population.
In this paper the population targets simulated represented the (STN) and Globus Pallidus Externus (GPe) but recent researches have a new target like the PPn .
where v j denoted the membrane potential of the j th neuron and w j was an auxiliary variable Eq. 4 describes the evolution of the parameter (the index s denotes synaptic), relevant for the synaptic coupling and modulating the dynamics of the current (see Eq. 10). The is calculated as a function of the membrane potential v j . The other parameters denote:
C: membrane capacitance;
g ca , g k, g l : leak, Ca++, and K+ conductances through membranes channel;
v ca , v k, v l : equilibrium potential of relevant ion channels;
v 1 , v 2, v 3, v 4, : tuning parameters for steady state and time constant.
The dynamics of the neuron was controlled by external currents: a slowly varying current , whose dynamics is described in Eq. 3; this current had been proposed by Rinzel and Ermentrout  as a source of bursting, which reflects the inhibitory feedback from the GPe.
where ξ(t) is Gaussian distributed, zero mean noise with unit standard deviation . The neuronal dynamics displayed by the simplified model STN neurons are characterized by a bursting activity, wherein the simulations the bursting frequency is chosen to be close to frequencies usually observed in patients suffering from parkinsonian tremor. The bursts are formed by a small number of spikes, six to ten per burst, which are controlled by the feedback from the GPe.
The bursting pattern displayed by the network elements is qualitatively similar to the bursting pattern observed in wide areas of the nervous system [35, 36, 26, 27]. For illustration, we represent the neurons in the populations as arranged in square lattices. Two stimulation protocols were used: the first one using one electrode positioned at the lattice center, and a second one consisting of four stimulation electrodes, equally spaced within the population.
For the numerical integration of the system we used the routine by Runge-Kutta method of order four.
Simplifying the complex organization of the neuronal population of the STN, we mimic excitatory couplings between different neurons of the population. In literature, the origin of the synchronized activity established between neurons in pathological conditions is not known, but could be induced by the modification of the bursting activity of the neurons.
Also, in our model, single spikes are not able to cause a significant change of the spiking times of the other neurons, only synchronized bursts are able to induce such a change in the firing pattern.
To design the connection within STN neurons we used a simple excitatory coupling ( ) that does not take into account any input from other BG nuclei. The synaptic interaction is modelled as suggested by Terman [37, 38].
The stimulation was applied via one or four electrodes located within the network. The strength of the stimulation typically decays with distance from the stimulation electrode. Which compartments of a neuron are activated by an electrical extracellular stimulation is still not clear [44, 45]. The final result of the stimulation is most probably a combination of excitatory action directly at the soma of the neuron together with an activation of afferent fibres [46, 31]. Therefore in our mathematical model we elaborated excitatory actions, which take place at the soma and within the population mediated by excitatory synapses connecting the STN neurons [47, 40, 30, 41, 48], and inhibitory actions coming from an activation of inhibitory afferent fibres [30–33].
Hence, the absolute value of the stimulation current was used to control the slowly varying current. Therefore, this type of stimulation mimicked the activation of afferent inhibitory fibres. The exact interplay between these different types of action remains unknown, to date, and is essential in determining the shape of the electrical pulses . To cope with this challenge we first investigated excitatory effects. The corresponding stimulation term was given in Eq. (1) as . Here s1 = 1 is used to weight the local effectiveness of the stimulation and the step function X1 defines the onset and offset of the excitatory stimulation acting directly at the soma of the stimulated neuron.
where ||y j - y k || is the distance between the k th electrode and j th neuron, and cs is the parameter that controls the strength of the stimulation. x1(t) (Eq. 1) determines the onset and offset of the excitatory stimulation. W(t) is the continuous pulse train consisting of a high-frequency repetition of the biphasic pulses presented through the electrode (Figure 2). c n is a normalization factor, which is dependent on the number of stimulation electrodes used and which guarantees that the total amount of stimulation remains independent on the number of electrodes used.
Here R is the measure of synchronization and Θ is mean phase. It results that 0 ≤ R(t) ≤ 1 for all times t; R = 1 corresponds to perfect in-phase synchronization, whereas R = 0 means complete desynchronization .
Stimulation consists in introducing a stimulation current that influences Eq. 1. In this equation we remind that the parameter X1 controls the onset and offset of the excitatory stimulation acting directly at the soma of the stimulated neuron. Moreover W(t) is a pulse train with a fixed high-frequency equal to 130 Hz, consisting of a repetition of the biphasic pulses delivered through the electrode. This approach was applied in both cases studied, i.e. when only one electrode is implanted at the centre of the network and when four electrodes are positioned symmetrically with respect to the lattice. Multi target stimulation reflects the common neurosurgical protocol. However, both the techniques (single and multi-electrodes) have the common goal of reducing the synchronized activity of the target population, either by re-establishing a normal physiological activity in a highly synchronized population of neurons or by reversibly mimicking a tissue lesioning (DBS).
Stimulation methods via one central electrode
The figure shows that, in absence of stimulation, the neurons bursting was strongly synchronized in-phase. At t = 4.5 s we introduced the stimulating signal: this caused a desynchronization for t = 4.5 s , leading a large decrease in R(t) (average value R(t) = 0.4). The pulse input, simulating DBS, acted so as to mask the synchronizing effects of the excitatory interconnections within the STN.
Stimulation methods via four electrodes
in the multi electrodes case, the stimulation strength (in terms of current in Eq. 11) decayed with distance from the stimulation site and each neuron received a mixture of the four stimulation signals.
The model of the neural network that we used to simulate the biological behavior, in the previous section, was made up of a limited number of neurons, imposed by the limitations of computational resources.
Till now we modelled the neural network using Morris Lecar neurons consisting of six differential equations: in these conditions, considering a reasonable computation burden, we could analyze a number of neurons not larger than 100. For this reason we decided to study a reduced order neuron model, with less differential equations but preserving the same characteristics.
We used the Izhikevich model equation, in particular the chattering model for bursting neurons . This model offers a reduced complexity, if compared with the Morris Lecar equations, retaining, at the same time, the most important features. In particular the neurons can fire with stereotypical bursts of closely spaced spikes .
All the other characteristics, like the synaptic dynamics (I jsyn ), the noise component (I jnoise ), the stimulation current (I jstim ) and other relevant parameters for the network, were left unchanged. The Izhikevich parameters were suitably modified in order to allow a (time) comparison with the Morris Lecar model.
Here u j is an auxiliary variable. The term was obtained by fitting the spike initiation dynamics of cortical neurons so that units of v j correspond to mV and units of time correspond to ms. We modulated the parameters so as to obtain realistic bursting patterns. The parameter a determines the rate of recovery, whereas b determines the sensitivity of the recovery variable u j to the membrane potential v j .
The parameter c determines the after-spike reset value, which depends on fast high threshold conductance. Similarly, the parameter d determines the after-spike reset of the recovery variable, which depends on slow high-threshold conductance.
The Eq. 16, was the same as Eq.4 related the system studied in the previous section. The I jnoise was calculated offline by using a second order algorithm with decay time τ noise . We used a numerical integration method based on the Eulero algorithm.
Here ||y j - y k || is the distance between the k th and the j th neuron and G r is a random variable with Gaussian distribution, zero mean and 0.1 standard deviation; it is added to the distance among neurons.
After the definition of the reduced neuron model and the corresponding network, we use the same stimulation protocol as that applied with the Morris Lecar system but with a different number of neurons.
Using numerical methods we investigated the dynamical properties of a neural population consisting of 250 neurons. In the simulations we used the same periodic train pulse discussed above. In the first case, as seen before, we stimulated with an electrode in the middle of the network and then with four electrodes positioned in the centre of the population. A further step was to connect two identical populations to appreciate the effect of the stimulation in both populations, resembling the effect of multiple nuclei, typical of DBS.
This shows that the neuron bursting was strongly synchronized in-phase and R(t) = 1 for t ∈ 0s. The mean synchronization measure was R(t) = 0.6 for t ∈ [3,6.5]s. In this case the network desynchronized after applying the stimulation signal.
These results, compared with those ones reported in the previous section, allow us to appreciate the same dynamical effect of the stimulation signals also in this reduced order network. In summary we used:
• Different, reduced order model consisting of Izhikevich- Chattering;
• Number of neurons equal to 225;
• Random network.
With this different system we obtained good values of index synchronization. For this reason we decided to continue the study with two coupled populations to study the mutual interaction effects.
To study the challenging situation of strong driving, we assume that the mean coupling within the driving population is equal to the mean coupling between the two populations. The excitatory coupling between the pacemaker and the driven system is given by . Within the driven population a weak excitatory synaptic coupling exists, which by itself does not induce synchronization. For illustration, we represent the neurons in the populations as arranged in square lattices and stimulated either via four stimulation electrodes equally spaced within the population or with one electrode situated in the central part of the network.
The stimulation in population 1 causes an instantaneous desynchronization of population 2, which shows a decrease of the synchronization measure R.
We can see that, when stimulated, the first population acts as peacemaker for population 2 (driven population), Figure 13, resulting in a decrease of the synchronization index R. For the first population R(t) = 0.7 (where t ∈ [3,6.5]s). For the second population, Figure 14, the value of R changes from an average of R(t) = 0.7 ∈ (when t 0s), to R(t) = 0.5 (when t ∈ [3,6.5]s), as a consequence of the coupling with population 1.
We can see in this simulation how the increasing the number of electrodes causes an increasing in the desynchronization effect. Stimulating the first population (Figure 15) acts as peacemaker for the driven population 2 (Figure 16).
For the first population we have R(t) = 0.2 where t ∈ [3,6.5]s, moreover for the second population the value of R change from R(t) = 0.5 with t ∈ 0s, to R(t) = 0.3 with t ∈ [3,6.5]s, given by coupling with population 1.
where r j is the distance between neuron j and the recording electrode . All the current components, except for the stimulation current, contributed to the ionic current I j (t)(Eq. 1). Therefore I j (t) is composed of , , , and the three ionic currents responsible for the basic spiking features of the j th STN model neurons (i.e., the right hand side of Eq. 1). Re = 1 is the extracellular resistivity that was assumed to be homogeneous.
In the following the effect of stimulation is considered. We apply the traditional stimulation via four electrodes in the same 100 neuron population.
Referring to one population, in order to compare the dynamics of the un-stimulated and the stimulated network of the Izhijevich and Morris Lecar network, we studied how the application of a stimulation signal caused a change in the LFP of a population. As seen before, we will appreciate a desynchonization of the LFP through the spectrogram in dB of Izhikevich network. From Eq.18 we calculated the LFP of the Izhikevich population.
From these plots we can see that there is a high synchronization at low frequency and this is a characteristic, mimicking the pathological network. This behavior is clearly visible also in the Izhikevich neuron population: between the [5-20] Hz range there is a high activity.
When we apply the pulse train at 130 Hz, the network desynchronizes at the low frequencies and we notice a spreading of the power density. The same effect is seen in the literature, when the patient is treated with L-Dopa , and in Figure 18, referring to the Morris Lecar network. If we compare the Figure 17 with Figure 19 and 18 with Figure 20, we can see that we obtained a same effects both in the stimulated and in the not stimulated case.
In this work some models for studying the effect of stimulation in neural populations mimicking the dynamics met in the relevant brain neural tissues of PD patients are analysed. A new structure, built using a reduced order neuron model (the Izhikevich one), computationally much lighter, was demonstrated to lead to the same results as other neural networks already introduced in the literature. This gives the opportunity to more efficiently study the effect of stimulation in large scale neural networks. Results were compared with the other mathematical models, i.e. Morris Lecar neurons, analyzing both parameters related to the neural level (the synchronization effect) and indexes related to the macroscopic level (LFP). The results agree in terms of the positive effect of the stimulation in terms of power spectral density. We consider these results very promising: they open the way to the possibility of simulating large scale networks in pathological conditions at the aim to design new control strategies for the PD effect mitigation.
List of abbreviations
- PD :
- DBS :
Deep Brain Stimulation
- BG :
- FDA :
Food and Drug Administration
- STN :
- GPi :
Globus Pallidus internus
- PPTg :
Nucleus Tegmenti Peduncolopontini
- PSD :
Power Spectral Density
- LFP :
Local Field Potential
- GPe :
Globus Pallidus externus
- VIM :
Ventral Intermediate Nucleus
- CM-Pf :
The authors would like to thank prof. Luigi Fortuna, Dean of the Engineering Faculty of the University of Catania, for his continuous support in this research activity.
- Jankovic J: Parkinson's disease: clinical features and diagnosis. J Neurol Neurosurg Psychiatr. 2008, 79 (4): 368-76. 10.1136/jnnp.2007.131045.View ArticleGoogle Scholar
- National Institute for Health and Clinical Excellence. Clinical guideline 35: Parkinson's disease. London. 2006, [http://www.nice.org.uk/]Google Scholar
- Cilia R: Long-term Efficacy of Entacapone in Patients with Parkinson's Disease and Motor Fluctuations - A Six-Year Clinical Follow-Up Study. 2006Google Scholar
- Samii , Ali : Stereotactic Surgery. American Parkinson Disease Association. 2002Google Scholar
- Tass PA: Pioneering book on the principles of phase resetting with a focus on potential applications in medicine and biology. Phase Resetting in Medicine and Biology. 1999, Springer Verlag, NY, USAView ArticleGoogle Scholar
- Tass PA, Hauptmann C, Popovych O: Development of therapeutic brain stimulation techniques with methods from nonlinear dynamics and statistical physics. Int J Bifurcat Chaos. 2006, 16 (7): 1889-1911). 10.1142/S0218127406015787.View ArticleMathSciNetGoogle Scholar
- Kringelbach ML, Jenkinson N, Owen SLF, Aziz TZ: Translational principles of deep brain stimulation. Nature Reviews Neuroscience. 2007, 8: 623-635. 10.1038/nrn2196.View ArticleGoogle Scholar
- Gildenberg PL: Evolution of neuromodulation. Stereotact Funct Neurosurg. 2005, 83 (2-3): 71-79. 10.1159/000086865.View ArticleGoogle Scholar
- U.S. Department of Health and Human Services: FDA approves implanted brain stimulator to control tremors. Retrieved October 18, 2006, [http://www.fda.com/]
- "Brain pacemaker" treats dystonia. 2003, KNBC TV, Retrieved October 18, 2006, [http://knbctvhealth.ip2m.com/]
- Deep brain stimulation. Surgery Encyclopedia. Retrieved January 25, 2007, [http://www.surgeryencyclopedia.com/Ce-Fi/Deep-Brain-Stimulation.html]
- Kleiner-Fisman G, Herzog J, Fisman DN: Subthalamic nucleus deep brain stimulation: summary and meta-analysis of outcomes. Mov Disord. 2006, 21 (Suppl 14): S290-304. 10.1002/mds.20962.View ArticleGoogle Scholar
- National Institute of Neurological Disorders and Stroke. Deep brain stimulation for Parkinson's Disease information page, Retrieved November 23, 2006, [http://www.ninds.nih.gov/]
- Moro E, Lang AE: Criteria for deep-brain stimulation in Parkinson's disease: review and analysis. Expert Review of Neurotherapeutics. 2006, 6 (11): 1695-705. 10.1586/14737126.96.36.1995. PMID 17144783View ArticleGoogle Scholar
- Apetauerova D, Ryan RK, Ro SI, Arle J: End of day dyskinesia in advanced Parkinson's disease can be eliminated by bilateral subthalamic nucleus or globus pallidus deep brain stimulation. Movement Disorders. 2006, 21 (8): 1277-9. 10.1002/mds.20896. PMID 16637040View ArticleGoogle Scholar
- Mazzone P, Scarnati E, Garcia-Rill E: Commentary: The pedunculopontine nucleus: clinical experience, basic questions and future directions. J Neural Transm. 2010.Google Scholar
- Mazzone P, Scarnati E: Deep brain stimulation of the medial thalamus for movement disorders: the role of centromedianparafascicular complex. Edited by: Krames ES, Peckham PH, Rezai AR. 2009, Neuromodulation. Academic Press, New York, 599-615.Google Scholar
- Mazzone P, Sposato S, Insola A, Scarnati E: The deep brain stimulation of the pedunculopontine tegmental nucleus: towards a new stereotactic neurosurgery. J Neural Trans. 2011Google Scholar
- Riederer P, Reichmann H, Youdim MBH, Gerlach M: Parkinson's Disease and Related Disorders. Journal of Neural Transmission. 2006, Supplementa 70:Google Scholar
- Gerfen CR: D1 and D2 dopamine receptor-regulated gene expression of striatonigral and striatopallidal neurons. Science. 1990, 250: 1429-1432. 10.1126/science.2147780.View ArticleADSGoogle Scholar
- Michael Cassidy, Paolo Mazzone, Antonio Oliviero, Angelo Insola, Pietro Tonali, Di Lazzaro Vincenzo, Peter Brown: Movement-related changes in synchronization in the human basal ganglia. Brain: a journal of neurology. 2002, 125 (Pt 6): 1235-46.Google Scholar
- Hauptmann C, Popovych O, Tass PA: Effectively desynchronizing deep brain stimulation based on a coordinated delayed feedback stimulation via several sites: a computational study. Biological Cybernetics. 2005, 93 (6): 463-470. 10.1007/s00422-005-0020-1.View ArticleMathSciNetGoogle Scholar
- Eugene M: Izhikevich: Simple Model of Spiking Neurons. IEEE Transactions on Neural Networks. 2003, 14 (6):
- Nini A, Feingold A, Slovin H, Bergmann H: Neurons in the globus pallidus do not show correlated activity in the normal monkey, but phase-locked oscillations appear in the MPTP model of parkinsonism. J Neurophysiol. 1995, 74: 1800-1805.Google Scholar
- Morris C, Lecar H: Voltage oscillations in the barnacle giant muscle fiber. Biophys J. 1981, 35: 193-213. 10.1016/S0006-3495(81)84782-0.View ArticleGoogle Scholar
- Luhmann HJ, Mudrick-Donnon LA, Mittmann T, Heinemann U: Ischaemiainduced long-term hyperexcitability in rat neocortex. Eur J Neurosci. 1995, 7: 180-191. 10.1111/j.1460-9568.1995.tb01054.x.View ArticleGoogle Scholar
- Zirh TA, Lenz FA, Reich SG, Dougherty PM: Patterns of bursting occurring in thalamic cells during parkinsonian tremor. Neuroscience. 1998, 83: 107-121. 10.1016/S0306-4522(97)00295-9.View ArticleGoogle Scholar
- Rinzel J, Ermentrout GB: Analysis of Neural Excitability and Oscillations. Methods in Neuronal Modelling From Synapses to Networks. Edited by: Koch CH, Segev I. 1989, MIT Press, MA, USA, 135-169.Google Scholar
- Hauptmann C, Mackey MC: Stimulus dependent onset latency of the inhibitory recurrent activity. Biol Cybern. 2003, 88: 459-467.Google Scholar
- Kitai ST, Kita H: Anatomy and physiology of the subthalamic nucleus: a driving force of the basal ganglia. 1987, The Basal Ganglia II: Structure and Function; Current Concepts. Plenum, NY, USA, 357-373.Google Scholar
- Wichmann T, Bergman H, Starr PA, Subramanian T, Watts RL, DeLong MR: Comparison of MPTP-induced changes in spontaneous neuronal discharge in the internal pallidal segment and in the substantia nigra pars reticulata in primates. Exp Brain Res. 1999, 125: 397-409. 10.1007/s002210050696.View ArticleGoogle Scholar
- Bolam JP, Hanley JJ, Booth PAC, Bevan MD: Synaptic organisation of the basal ganglia. J Anat. 2000, 196: 527-542. 10.1046/j.1469-7580.2000.19640527.x.View ArticleGoogle Scholar
- Hamani C, Saint-Cyr JA, Fraser J, Kaplitt M, Lozano AM: The subthalamic nucleus in the context of movement disorders. Brain. 2004, 127: 4-20. 10.1093/brain/awh029.View ArticleGoogle Scholar
- Dolan K, Witt A, Spano ML, Neiman A, Moss F: Surrogates for finding unstable periodic orbits in noisy data sets. Phys Rev E. 1999, 59: 5235-5241. 10.1103/PhysRevE.59.5235.View ArticleADSGoogle Scholar
- Magill PJ, Bolam JP, Bevan MD: Dopamine regulates the impact of the cerebral cortex on the subthalamic nucleus-globus pallidus network. Neuroscience. 2001, 106 (2): 313-320. 10.1016/S0306-4522(01)00281-0.View ArticleGoogle Scholar
- Beurrier C, Garcia L, Bioulac B, Hammond C: Subthalamic nucleus: a clock inside basal ganglia?. Thalamus Relat Syst. 2002, 2: 1-8.View ArticleGoogle Scholar
- Rubin D, Terman J: High frequency stimulation of the subthalamic nucleus eliminates pathological thalamic rhythmicity in a computational model. J Comput Neurosci. 2004, 16: 211-235. 10.1023/B:JCNS.0000025686.47117.67.View ArticleGoogle Scholar
- Terman D, Rubin JE, Yew AC, Wilson CJ: Activity patterns in a model for the subthalamo pallidal network of the basal ganglia. J Neurosci. 2002, 22: 2963-2976.Google Scholar
- Beurrier C, Congar P, Bioulac B, Hammond C: Subthalamic nucleus neurons switch from single-spike activity to burstfiring mode. J Neurosci. 1999, 19 (2): 599-609.Google Scholar
- Levesque J-C, Parent A: GABAergic interneurons in human subthamic nucleus. Mov Disord. 2005, 20: 574-584. 10.1002/mds.20374.View ArticleGoogle Scholar
- Humphries MD, Gurney KN: A pulsed neural network model of bursting in the basal ganglia. Neural Netw. 2001, 14: 845-863. 10.1016/S0893-6080(01)00060-0.View ArticleGoogle Scholar
- Traub RD, Miles R: Neural Networks of the Hippocampus. 1991, Cambridge University Press, Cambridge, UKView ArticleGoogle Scholar
- Hellwig B: A quantitative analysis of the local connectivity between pyramidal neurons in layers 2/3 of the rat visual cortex. Biol Cybern. 2000, 82: 111-121. 10.1007/PL00007964.View ArticleGoogle Scholar
- Benabid AL, Benazzous A, Pollak P: Mechanisms of deep brain stimulation. Mov Disord. 2002, 17: 73-74. 10.1002/mds.10145.View ArticleGoogle Scholar
- Nowak LG, Bullier J: Axons, but not cell bodies, are activated by electrical stimulation in cortical gray matter. II. Evidence from selective inactivation of cell bodies and axon initial segments. Exp Brain Res. 1998, 118: 489-500. 10.1007/s002210050305.View ArticleGoogle Scholar
- Grill WM, McIntyre CC: Extracellular excitation of central neurons: implications for the mechanisms of deep brain stimulation. Thalamus Relat Syst. 2001, 1: 269-277.View ArticleGoogle Scholar
- Shen K, Zhu Z, Munhall A, Johnson SW: Synaptic plasticity in rat subthalamic nucleus induced by high-frequency stimulation. Synapse. 2003, 50: 314-319. 10.1002/syn.10274.View ArticleGoogle Scholar
- Gillies A, Willshaw D: Models of the subthalamic nucleus - the importance of intranuclear connectivity. Med Eng Phys. 2004, 26: 723-732. 10.1016/j.medengphy.2004.06.003.View ArticleGoogle Scholar
- Deep Brain Stimulation. Department of Neurological Surgery, University of Pittsburgh, Retrieved May 13, 2008, [http://www.neurosurgery.pitt.edu/imageguided/movement/stimulation.html]
- Ranck JB: Which elements are excited in electrical stimulation of mammalian central nervous system: a review. Brain Res. 1975, 98: 417-440. 10.1016/0006-8993(75)90364-9.View ArticleGoogle Scholar
- Izhikevich Eugene: Which Model to Use for Cortical Spiking Neurons?. IEEE Transactions on Neural Networks. 2004, 15 (5): 10.1109/TNN.2004.832719.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.