Phase-dependent modulation as a novel approach for therapeutic brain stimulation

Azodi-Avval, Ramin; Gharabaghi, Alireza

doi:10.3389/fncom.2015.00026

ORIGINAL RESEARCH article

Front. Comput. Neurosci., 26 February 2015
Volume 9 - 2015 | https://doi.org/10.3389/fncom.2015.00026

Phase-dependent modulation as a novel approach for therapeutic brain stimulation

Ramin Azodi-Avval^1,2

Alireza Gharabaghi^1,2^*

¹Division of Functional and Restorative Neurosurgery and Division of Translational Neurosurgery, Department of Neurosurgery, Eberhard Karls University, Tuebingen, Germany
²Neuroprosthetics Research Group, Werner Reichardt Centre for Integrative Neuroscience, Eberhard Karls University, Tuebingen, Germany

Closed-loop paradigms provide us with the opportunity to optimize stimulation protocols for perturbation of pathological oscillatory activity in brain-related disorders. In this vein, spiking activity of motor cortex neurons and beta activity of local field potentials in the subthalamic nucleus have both been used independently of each other as neuronal signals to trigger deep brain stimulation for alleviating Parkinsonism. These approaches were superior to the standard continuous high-frequency stimulation protocols used in daily practice. However, they achieved their effects by bursts of stimulation that were applied at high-frequency as well, i.e., independent of the phase information in the stimulated region. In this context, we propose that, by timing stimulation pulses relative to the ongoing oscillation, an alternative approach, namely the targeted perturbation of pathological rhythms, could be obtained. In this modeling study, we first captured the underlying dynamics of neuronal oscillations in the human subthalamic nucleus by phased coupled neuronal oscillators. We then quantified the nature of the interaction between these coupled oscillators by obtaining a physiologically informed phase response curve from local field potentials. Reconstruction of the phase response curve predicted the sensitivity of the phase oscillator to external stimuli, revealing phase intervals that optimally maximized the degree of perturbation. We conclude that our specifically timed intervention based on the coupled oscillator concept will enable us to identify personalized ways of delivering stimulation pulses in closed-loop paradigms triggered by the phase of pathological oscillations. This will pave the way for novel physiological insights and substantial clinical benefits. In addition, this precisely phased modulation may be capable of modifying the effective interactions between oscillators in an entirely new manner.

Introduction

Brain neuromodulation by deep brain stimulation (DBS) is meanwhile a recognized form of treatment for several neurological and neuropsychiatric disorders such as severe Parkinson's disease (PD) (Schuepbach et al., 2013), dystonia (Vidailhet et al., 2005), and essential tremor (Deuschl et al., 2011). However, general application of this therapeutic modality remains limited. This might be due to stimulation-induced side effects and/or partial efficacy of the intervention which is probably related to a misalignment between stimulation parameters and the current disease state (Moro et al., 2006; Mure et al., 2011). While patients often display variable clinical symptoms, the stimulation parameters, such as continuous high-frequency stimulation, are predefined and remain unchanged until manual modifications are performed by the physician in charge.

Closed-loop paradigms modulating the stimulation parameters on the basis of online recorded physiological markers provide us with the opportunity to adjust stimulation protocols and improve therapeutic efficacy. In this regard, the first studies in both non-human primates (Rosin et al., 2011) and Parkinsonian patients (Little et al., 2013) addressed the current limitations by applying stimulation in an adaptive manner only when specific physiological markers were detected. More specifically, adaptive DBS of the globus pallidus internus controlled by spiking activity of motor cortex neurons was more effective than continuous high-frequency stimulation in a non-human primate model of PD (Rosin et al., 2011). In a recent clinical study, adaptive DBS of the subthalamic nucleus (STN)—triggered by beta-band activity (a physiological marker of motor impairment in PD) and recorded in the immediate vicinity of the stimulating electrode in the STN—was shown to be more energy-efficient than, and clinically superior to continuous high-frequency DBS (Little et al., 2013). Despite being superior to the standard stimulation protocols, these closed-loop approaches nonetheless achieved their effects by bursts of stimulation applied at high frequency independent of the phase information in the stimulated region (Rosin et al., 2011; Little et al., 2013).

In this context, we propose that, by timing stimulation pulses relative to the ongoing oscillation in the respective area, an alternative approach—namely the targeted perturbation of pathological rhythms—can be introduced. We therefore performed a modeling study to determine the essential dynamics of neuronal oscillations in the human subthalamic nucleus by phased coupled neuronal oscillators in an aim to identify those phase intervals in which stimulation would maximize the degree of perturbation of pathological rhythms.

Materials and Methods

The present modeling study is based on intraoperative electrophysiological recordings of two PD patients who underwent standard DBS surgery with bilateral electrode implantation in the STN and was performed in accordance with the guidelines of the local ethics committee of the Medical Faculty of the University of Tuebingen. This data represents spontaneous brain activity recorded for about 3 min through the final quadripolar DBS electrode (model 3389, Medtronic, Inc., Minneapolis, MN) which was implanted in one brain hemisphere while electrode insertion was prepared for the second side. The STN was localized via direct targeting on preoperative magnetic resonance imaging (Foltynie et al., 2011) and then during surgery with online electrophysiology (Chen et al., 2006; Holdefer et al., 2010). In accordance with standard operating procedures, dopamine medication was not administered for the last time 12 or more hours prior to surgery to avoid interference with intraoperative recordings and clinical testing (Hammond et al., 2007). In addition, intraoperative propofol medication was discontinued about 30 min before electrophysiological recordings were initiated (Raz et al., 2010). Local field potential (LFP) signals were continuously sampled at 1.4 KHz and amplified by factor 50. The reference electrodes were attached to the ears and the ground was placed on the nasion. The impedance for the intracranial electrode was ~1KΩ. Stereotactic planning (Machado et al., 2006) was based on multimodal preoperative images (1 mm slice thickness) from contrast-enhanced magnetic resonance imaging (MRI) and computer tomography (CT) imaging (Siemens, Erlangen, Germany). Standard electrophysiological recordings (AlphaOmega, Nazareth, Israel) and clinical test stimulation were performed intraoperatively to adjust electrode localization. This was later confirmed by postoperative MR and CT imaging.

Phase Synchronization

The implanted quadripolar electrode (Figure 1A) enabled us to independently record four oscillatory sources, i.e., local field potentials (LFP), in the target region and to quantify their neural interaction. For this purpose, we used the phase approach proposed by Rosenblum and Pikovsky (2001) to reconstruct the phase response curve (PRC) of the LFP signals recorded in the STN. The concept of phase synchronization between coupled oscillators has not been used to derive a phase-response curve of brain signals yet, but has already been applied in several electrophysiological studies in different functional systems (Kralemann et al., 2013; Zhu et al., 2013). However, brain signals are particularly suited for such an approach as filtered LFP signals at a given frequency band have a sinusoidal waveform similar to an oscillator. Therefore, phase modeling might capture the neural network dynamics of coupled oscillators (Wang, 2010).

FIGURE 1

Figure 1. (A) Deep brain stimulation lead positioned in the subthalamic nucleus region with four electrode contacts (0–3) and respective local field potential (LFP) recordings. (B) The LFPs from the electrode contact (black) with the most pronounced oscillatory activity in the beta-band and the electrode contact (red) with the most prominent phase synchronization to the first contact are captured as coupled neural oscillators. In a second step, the phase of these oscillators is extracted by Hilbert transform. (C) Illustration of the typical characteristics of a phase response curve (PRC) with positive (orange) or negative (blue) values indicating an earlier or later start of the perturbated cycle, respectively, i.e., the perturbation could cause a positive (orange) phase shift (phase advance), or a negative (blue) phase shift (phase delay), depending on the timing of the perturbating stimulus. The black curve represents the ongoing oscillation in the absence of any perturbation.

Data analysis of LFPs were performed in Matlab (The Mathworks, Natick, Massachusetts, USA) using custom-made scripts, DAMOCO toolbox (Kralemann et al., 2007, 2008) and Fieldtrip, an open source analysis toolbox (Oostenveld et al., 2011). To reconstruct the PRC, the power of LFP was first computed by the multitaper method. Maximum spectral estimation was provided by a single Hanning taper (Mitra and Pesaran, 1999). This approach enabled us to identify the specific electrode contact that displayed pronounced oscillatory activity in the beta-band—a known pathophysiological marker of motor impairment in PD (Little and Brown, 2012)—and which we used as a reference for the further computing. We next applied the Weighted Phase Lag Index (WPLI) (Vinck et al., 2011) to detect the electrode contact displaying the most prominent phase synchronization with this reference electrode. WPLI is defined as follows:

\begin{matrix} \begin{array}{l} W P L I_{L F P - L F P} \\ = \frac{n^{- 1} \sum_{n = 1}^{N} | i m a g (S_{L F P - L F P n}) | s g n (i m a g (S_{L F P - L F P n}))}{n^{- 1} \sum_{n = 1}^{N} i m a g (S_{L F P - L F P n})} \end{array} & (1) \end{matrix}

where S is the cross-spectrum density matrix, imag the imaginary part and sgn the sign function. Moreover, WPLI is insensitive to volume conduction effects since it eliminates zero phase lag signals. Thereby, it is more sensitive to detect true phase interaction as compared to common phase measurements such as the imaginary part of coherence (Nolte et al., 2004). The value of WPLI was standardized by an estimate of its standard deviation and values beyond threshold of 3 (corresponding p < 0.003) were considered statistically significant (Nolte et al., 2004; Hohlefeld et al., 2013, 2014). This statistically significant frequency range (see Figure 2B) was used for further analysis, i.e., phase extraction.

FIGURE 2

Figure 2. (A) Power spectrum of LFPs recordings at the four different contacts (0–3) for patient 1 and 2, respectively, revealed that the most upper electrode contact (3, black) of the quadripolar lead has the highest spectral power in the beta-frequency band. (B) Computing phase synchronization between electrode contact 3 (black) and the others applying the Weighted Phase Lag Index (WPLI) showed that the most prominent coupling is with contact 2 (red). The statistically significant frequency range for WPLI is indicated in gray. The dashed line represents the significance threshold (3, corresponding p < 0.003).

Physiologically Informed Phase Response Curve

Having determined the significant frequency interval of neural interaction for each patient, we went on to filter the LFP in that particular, i.e., patient-specific, frequency range using a Kaiser FIR filter with the MATLAB filtfilt function to avoid phase distortions. The unwrapped phase of each sample of the filtered LFP was then computed using the Hilbert transform (Figure 1B). The method proposed by Rosenblum and collaborators (Kralemann et al., 2007, 2008) was then used to empirically reconstruct the phase coupling function between the two electrode contacts in the STN that had already been determined. This method adapts the empirical phases to a generic model of two coupled oscillators as follows:

\begin{matrix} \begin{array}{l} {\dot{φ}}_{1} = ω_{1} + F_{2 \to 1} (φ_{1}, φ_{2}) \\ {\dot{φ}}_{2} = ω_{2} + F_{1 \to 2} (φ_{2}, φ_{1}) \end{array} & (2) \end{matrix}

where the dot represents the derivative, while φ and ω are the phase and autonomous frequency, respectively. Since the extracted phase is non-universe θ_1,2, a transformation is needed to provide an invariant description of the coupled dynamics. The transformation from the unwrapped phase to the genuine phase of N observations is defined as following,

\begin{matrix} φ = θ + 2 π \sum_{n \neq 1} \frac{S_{n}}{i n} (e^{i n θ} - 1) & (3) \end{matrix}

where

S_{n} = \frac{1}{N} \sum_{j = 1}^{N} e^{i n θ (j)} .

After this phase correction is performed, the coupling function F of two coupled systems can be approximated by the Fourier series:

\begin{matrix} F = \sum_{n, m} A_{n, m} e^{- i (n φ_{2} + m φ_{1})} & (4) \end{matrix}

with coefficients

A_{n, m} = \int_{0}^{2 π} \int φ_{1} e^{- i n φ_{1} - i m φ_{2}}

where n and m are indices synonymous to the n:m phase locking the indices of two oscillators and A refers to the coefficients of the respective Fourier series (Kralemann et al., 2007, 2008). The resulting PRC reflects the interaction of the oscillators. When the PRC results in positive or negative values (Smeal et al., 2010), it indicates the cycle to start sooner or later, respectively, i.e., a perturbation could cause a positive phase shift (phase advance), or a negative phase shift (phase delay) (Figure 1C).

Results

The recording of local field potentials (LFP) revealed that the most upper electrode contact of the quadripolar lead has the highest spectral power in the beta-frequency band (15–30 Hz, Figure 2A), indicating that it is located in the sensorimotor part of the STN (Figure 2A) (Holdefer et al., 2010; Yoshida et al., 2010; Zaidel et al., 2010; Novak et al., 2011; Deffains et al., 2014). Using this electrode contact as a reference, the neighboring contact revealed the most prominent phase synchronization when applying the Weighted Phase Lag Index (WPLI). This phase coherence showed a significant synchronization in the frequency range between 23–27 and 17–22 Hz for P1 and P2, respectively (Figure 2B), indicating that there is a pathological increase of functional connectivity in the sensorimotor part of the STN (Pogosyan et al., 2010). By applying a modeling approach of two neural oscillators characterized by phase dependency, we succeeded in quantifying this interaction by a phase response curve. The inherent variability of the PCRs resulting from the empirical data was resolved using the similarity method (Kralemann et al., 2013). PRC extracted from both P1 and P2 showed a similar pattern with two different domains (Figure 3): a phase delay in the 0 < φ < ~0.3 and ~0.7 < φ < 1.0 intervals and a strong phase advance in the ~0.3 < φ < 0.7 interval. The strongest phase-coupling effect was observed around φ = 0.6.

FIGURE 3

Figure 3. Empirically determined PRC phase response curves (PRC) of patient 1 and 2, respectively, with the averaged curve indicated in red. The strongest phase-coupling effect in both P1 and P2 was observed around φ = 0.6 with the highest positive value, i.e., the maximum phase advance.

Discussion

The beta-band functional connectivity observed within the STN tallies well-with earlier work, indicating that this synchronization reflects a pathological marker for different motor features of PD (Pogosyan et al., 2010; Hohlefeld et al., 2013). We also located the most prominent phase synchronization in the dorsal, sensorimotor part of the STN (Pogosyan et al., 2010), i.e., the area known to be clinically most effective in suppressing PD symptoms during DBS (Herzog et al., 2004; Yokoyama et al., 2006; Schlaier et al., 2014). It remains to be experimentally disentangled, whether additional physiological markers may capture a broader spectrum of PD symptoms and may therefore be better suited for closed-loop applications. However, more sophisticated approaches of simultaneous sensing and stimulation with online signal processing may decrease the battery life span and should therefore be balanced with potential clinical benefits (Little and Brown, 2014).

On the basis of these findings, phase-specific interventions appear to be the most straightforward approach for specifically perturbating pathological synchronization. The PRC disentangles the timing of the regional interaction in the dorsal STN to reveal the highest synchronization between the two oscillators at phase 0.6 (Figure 4). In fact, the PRC also reveals phase sensitivity of an oscillator to external perturbation (Smeal et al., 2010). We therefore propose that the application of DBS stimuli at this phase, i.e., replacing one of the two oscillators by an external stimulus, maximizes perturbation of the pathological state. Along these lines, the effect of a perturbation on an oscillating system is known to depend on the phase at which the perturbation is applied (Pikovsky et al., 2001; Kralemann et al., 2013). This assumption is further supported by physiologically calibrated modeling results which suggest that precisely timed stimulation pulses could indeed be used to shift the phase of oscillations (Witt et al., 2013).

FIGURE 4

Figure 4. Illustration of a closed-loop modulation protocol, recording and stimulating with two adjacent DBS electrode contacts in the STN, triggered by the ongoing oscillatory phase estimated on the basis of the PRC curve. To determine this precise timing, a closed-loop system will require a physiological calibration for each patient including the estimation of WPLI and PRC based on the individual vivo recordings. After this calibration, i.e., based on the subject-specific frequency range and the precise phase of maximum perturbation, the respective phase information will be the feedback signal to control the stimulation.

Moreover, when stimuli are properly phased with respect to ongoing oscillations, they might even induce long-term potentiation/depotentiation (LTP/LTD) effects (Martin et al., 2000; Kauer and Malenka, 2007): both in vitro (Huerta and Lisman, 1993) and in vivo (Pavlides et al., 1988) experimental studies reported LTD induction when stimuli were applied during the positive phase of the theta rhythm. Similarly, depotentiation of existing LTP was achieved when stimuli were phase-locked to the negative phase of theta (Huerta and Lisman, 1995). Phase-specific pulses may therefore induce bidirectional modifications of synaptic strength (Martin et al., 2000) with the potential to turn a whole network into a synchronized or a desynchronized state (Pfister and Tass, 2010).

Future closed-loop modulation protocols will require simultaneous recording and stimulation. The presented findings suggest that two adjacent DBS electrode contacts in the STN may be used for this purpose, i.e., to record and to stimulate, respectively. In such a scenario the stimulating contact would be triggered by the ongoing oscillatory phase recorded at the adjacent electrode and estimated on the basis of the PRC curve. To determine the precise timing, a physiological calibration of the closed-loop system will be required for each patient beforehand, including the estimation of WPLI and PRC based on the individual vivo recordings. After this calibration, i.e., based on the subject-specific frequency range and the precise phase of maximum perturbation, the respective phase information will be the feedback signal to control the stimulation. This approach can be implemented with online algorithms, i.e., real-time calculation of the instantaneous phase based on the Hilbert transform (Figure 4).

Phase-dependent stimulation protocols may therefore be capable of modifying the effective interactions between oscillators in an entirely new fashion, potentially inducing lasting effects mediated by LTP/LTD. Furthermore, with regard to different neurological and neuropsychiatric disorders, such an approach would transform brain stimulation from a symptomatic intervention and temporary modulation—displaying its effects only for the duration of stimulation—to a treatment option that induces long-term plastic changes and durable effects lasting beyond stimulation.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

RA was supported by the Graduate Training Centre of Neuroscience, International Max Planck Research School, Tuebingen, Germany. AG was supported by grants from the German Research Council (DFG EC 307), and from the Federal Ministry for Education and Research (BFNT 01GQ0761, BMBF 16SV3783, BMBF 03160064B, BMBF V4UKF014). We acknowledge support by Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of University of Tübingen.

References

Chen, C. C., Pogosyan, A., Zrinzo, L. U., Tisch, S., Limousin, P., Ashkan, K., et al. (2006). Intra-operative recordings of local field potentials can help localize the subthalamic nucleus in Parkinson's disease surgery. Exp. Neurol. 198, 214–221. doi: 10.1016/j.expneurol.2005.11.019

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Deffains, M., Holland, P., Moshel, S., de Noriega, F. R., Bergman, H., and Israel, Z. (2014). Higher neuronal discharge rate in the motor area of the subthalamic nucleus of Parkinsonian patients. J. Neurophysiol. 112, 1409–1420. doi: 10.1152/jn.00170.2014

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Deuschl, G., Raethjen, J., Hellriegel, H., and Elble, R. (2011). Treatment of patients with essential tremor. Lancet Neurol. 10, 148–161. doi: 10.1016/S1474-4422(10)70322-7

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Foltynie, T., Zrinzo, L., Martinez-Torres, I., Tripoliti, E., Petersen, E., Holl, E., et al. (2011). MRI-guided STN DBS in Parkinson's disease without microelectrode recording: efficacy and safety. J. Neurol. Neurosurg. Psychiatry 82, 358–363. doi: 10.1136/jnnp.2010.205542

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Hammond, C., Bergman, H., and Brown, P. (2007). Pathological synchronization in Parkinson's disease: networks, models and treatments. Trends Neurosci. 30, 357–364. doi: 10.1016/j.tins.2007.05.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Herzog, J., Fietzek, U., Hamel, W., Morsnowski, A., Steigerwald, F., Schrader, B., et al. (2004). Most effective stimulation site in subthalamic deep brain stimulation for Parkinson's disease. Mov. Disord. 19, 1050–1054. doi: 10.1002/mds.20056

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Hohlefeld, F. U., Huchzermeyer, C., Huebl, J., Schneider, G. H., Brücke, C., Schönecker, T., et al. (2014). Interhemispheric functional interactions between the subthalamic nuclei of patients with Parkinson's disease. Eur. J. Neurosci. 40, 3273–3283. doi: 10.1111/ejn.12686

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Hohlefeld, F. U., Huchzermeyer, C., Huebl, J., Schneider, G. H., Nolte, G., Brücke, C., et al. (2013). Functional and effective connectivity in subthalamic local field potential recordings of patients with Parkinson's disease. Neuroscience 250, 320–332. doi: 10.1016/j.neuroscience.2013.07.028

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Holdefer, R. N., Cohen, B. A., and Greene, K. A. (2010). Intraoperative local field recording for deep brain stimulation in Parkinson's disease and essential tremor. Mov. Disord. 25, 2067–2075. doi: 10.1002/mds.23232

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Huerta, P. T., and Lisman, J. E. (1993). Heightened synaptic plasticity of hippocampal CA1 neurons during a cholinergically induced rhythmic state. Nature 364, 723–725. doi: 10.1038/364723a0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Huerta, P. T., and Lisman, J. E. (1995). Bidirectional synaptic plasticity induced by a single burst during cholinergic theta oscillation in CA1 in vitro. Neuron 15, 1053–1063. doi: 10.1016/0896-6273(95)90094-2

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Kauer, J. A., and Malenka, R. C. (2007). Synaptic plasticity and addiction. Nat. Rev. Neurosci. 8, 844–858. doi: 10.1038/nrn2234

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Kralemann, B., Cimponeriu, L., Rosenblum, M., Pikovsky, A., and Mrowka, R. (2007). Uncovering interaction of coupled oscillators from data. Phys. Rev. E 76, 055201. doi: 10.1103/PhysRevE.76.055201

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Kralemann, B., Cimponeriu, L., Rosenblum, M., Pikovsky, A., and Mrowka, R. (2008). Phase dynamics of coupled oscillators reconstructed from data. Phys. Rev. E 77:066205. doi: 10.1103/PhysRevE.77.066205

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Kralemann, B., Frühwirth, M., Pikovsky, A., Rosenblum, M., Kenner, T., Schaefer, J., et al. (2013). In vivo cardiac phase response curve elucidates human respiratory heart rate variability. Nat. Commun. 4:2418. doi: 10.1038/ncomms3418

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Little, S., and Brown, P. (2012). What brain signals are suitable for feedback control of deep brain stimulation in Parkinson's disease? Ann. N. Y. Acad. Sci. 1265, 9–24. doi: 10.1111/j.1749-6632.2012.06650.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Little, S., and Brown, P. (2014). Focusing brain therapeutic interventions in space and time for Parkinson's disease. Curr. Biol. 24, R898–R909. doi: 10.1016/j.cub.2014.08.002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Little S., Pogosyan A., Neal S., Zavala B., Zrinzo L., Hariz M., et al. (2013). Adaptive deep brain stimulation in advanced Parkinson disease. Ann. Neurol. 74, 449–457. doi: 10.1002/ana.23951

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Machado, A., Rezai, A. R., Kopell, B. H., Gross, R. E., Sharan, A. D., and Benabid, A.-L. (2006). Deep brain stimulation for Parkinson's disease: surgical technique and perioperative management. Mov. Disord. 21, S247–S258. doi: 10.1002/mds.20959

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Martin, S. J., Grimwood, P. D., and Morris, R. G. M. (2000). Synaptic plasticity and memory: an evaluation of the hypothesis. Annu. Rev. Neurosci. 23, 649–711. doi: 10.1146/annurev.neuro.23.1.649

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Mitra, P. P., and Pesaran, B. (1999). Analysis of dynamic brain imaging data. Biophys. J. 76, 691–708. doi: 10.1016/S0006-3495(99)77236-X

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Moro, E., Poon, Y. W., Lozano, A. M., Saint-Cyr, J. A., and Lang, A. E. (2006). Subthalamic nucleus stimulation: improvements in outcome with reprogramming. Arch. Neurol. 63, 1266–1272. doi: 10.1001/archneur.63.9.1266

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Mure, H., Hirano, S., Tang, C. C., Isaias, I. U., Antonini, A., Ma, Y., et al. (2011). Parkinson's disease tremor-related metabolic network: characterization, progression, and treatment effects. NeuroImage 54, 1244–1253. doi: 10.1016/j.neuroimage.2010.09.028

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Nolte, G., Bai, O., Wheaton, L., Mari, Z., Vorbach, S., and Hallett, M. (2004). Identifying true brain interaction from EEG data using the imaginary part of coherency. Clin. Neurophysiol. 115, 2292–2307. doi: 10.1016/j.clinph.2004.04.029

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Novak, P., Przybyszewski, A. W., Barborica, A., Ravin, P., Margolin, L., and Pilitsis, J. G. (2011). Localization of the subthalamic nucleus in Parkinson disease using multiunit activity. J. Neurol. Sci. 310, 44–49. doi: 10.1016/j.jns.2011.07.027

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Oostenveld, R., Fries, P., Maris, E., and Schoffelen, J.-M. (2011). FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Intell. Neurosci. 2011, 1–9. doi: 10.1155/2011/156869

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Pavlides, C., Greenstein, Y. J., Grudman, M., and Winson, J. (1988). Long-term potentiation in the dentate gyrus is induced preferentially on the positive phase of θ-rhythm. Brain Res. 439, 383–387. doi: 10.1016/0006-8993(88)91499-0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Pfister, J.-P., and Tass, P. (2010). STDP in oscillatory recurrent networks: theoretical conditions for desynchronization and applications to deep-brain stimulation. Front. Comput. Neurosci. 4:22. doi: 10.3389/fncom.2010.00022

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Pikovsky, A., Rosenblum, M., and Kurths, J. (2001). Synchronization: A Universal Concept in Nonlinear Sciences. Cambridge: Cambridge University Press.

Pogosyan, A., Yoshida, F., Chen, C. C., Martinez-Torres, I., Foltynie, T., Limousin, P., et al. (2010). Parkinsonian impairment correlates with spatially extensive subthalamic oscillatory synchronization. Neuroscience 171, 245–257. doi: 10.1016/j.neuroscience.2010.08.068

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Raz, A., Eimerl, D., Zaidel, A., Bergman, H., and Israel, Z. (2010). Propofol decreases neuronal population spiking activity in the subthalamic nucleus of parkinsonian patients. Anesth. Anal. 111, 1285–1289. doi: 10.1213/ANE.0b013e3181f565f2

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Rosenblum, M. G., and Pikovsky, A. S. (2001). Detecting direction of coupling in interacting oscillators. Phys. Rev. E 64:045202. doi: 10.1103/PhysRevE.64.045202

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Rosin, B., Slovik, M., Mitelman, R., Rivlin-Etzion, M., Haber, S. N., Israel, Z., et al. (2011). Closed-loop deep brain stimulation is superior in ameliorating parkinsonism. Neuron 72, 370–384. doi: 10.1016/j.neuron.2011.08.023

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Schlaier, J., Hanson, C., Janzen, A., Fellner, C., Hochreiter, A., Proescholdt, M., et al. (2014). Deep brain stimulation in Parkinson's disease: motor effects relative to the MRI-defined STN. Neurosurg. Rev. 37, 461–471. doi: 10.1007/s10143-014-0521-2

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Schuepbach, W. M. M., Rau, J., Knudsen, K., Volkmann, J., Krack, P., Timmermann, L., et al. (2013). Neurostimulation for Parkinson's disease with early motor complications. N. Eng. J. Med. 368, 610–622. doi: 10.1056/NEJMoa1205158

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Smeal, R. M., Ermentrout, G. B., and White, J. A. (2010). Phase-response curves and synchronized neural networks. Philos. Trans. R. Soc. Lond. B Biol. Sci. 365, 2407–2422. doi: 10.1098/rstb.2009.0292

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Vidailhet, M., Vercueil, L., Houeto, J.-L., Krystkowiak, P., Benabid, A.-L., Cornu, P., et al. (2005). Bilateral deep-brain stimulation of the globus pallidus in primary generalized dystonia. N. Eng. J. Med. 352, 459–467. doi: 10.1056/NEJMoa042187

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Vinck, M., Oostenveld, R., van Wingerden, M., Battaglia, F., and Pennartz, C. M. A. (2011). An improved index of phase-synchronization for electrophysiological data in the presence of volume-conduction, noise and sample-size bias. NeuroImage 55, 1548–1565. doi: 10.1016/j.neuroimage.2011.01.055

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Wang, X.-J. (2010). Neurophysiological and computational principles of cortical rhythms in cognition. Physiol. Rev. 90, 1195–1268. doi: 10.1152/physrev.00035.2008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Witt, A., Palmigiano, A., Neef, A., El Hady, A., Wolf, F., and Battaglia, D. (2013). Controlling the oscillation phase through precisely timed closed-loop optogenetic stimulation: a computational study. Front. Neural Circuits 7:49. doi: 10.3389/fncir.2013.00049

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Yokoyama, T., Ando, N., Sugiyama, K., Akamine, S., and Namba, H. (2006). Relationship of stimulation site location within the subthalamic nucleus region to clinical effects on parkinsonian symptoms. Stereotact. Funct. Neurosurg. 84, 170–175. doi: 10.1159/000094956

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Yoshida, F., Martinez-Torres, I., Pogosyan, A., Holl, E., Petersen, E., Chen, C. C., et al. (2010). Value of subthalamic nucleus local field potentials recordings in predicting stimulation parameters for deep brain stimulation in Parkinson's disease. J. Neurol. Neurosurg. Psychiatry 81, 885–889. doi: 10.1136/jnnp.2009.190918

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Zaidel, A., Spivak, A., Grieb, B., Bergman, H., and Israel, Z. (2010). Subthalamic span of β oscillations predicts deep brain stimulation efficacy for patients with Parkinson's disease. Brain 133(Pt 7), 2007–2021. doi: 10.1093/brain/awq144

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Zhu, Y., Hsieh, Y.-H., Dhingra, R. R., Dick, T. E., Jacono, F. J., and Galán, R. F. (2013). Quantifying interactions between real oscillators with information theory and phase models: application to cardiorespiratory coupling. Phys. Rev. E 87:022709. doi: 10.1103/PhysRevE.87.022709

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Keywords: brain state-dependent stimulation, closed-loop stimulation, biophysical model, coupled neuronal oscillators, phase response curve, Parkinson's disease, deep brain stimulation

Citation: Azodi-Avval R and Gharabaghi A (2015) Phase-dependent modulation as a novel approach for therapeutic brain stimulation. Front. Comput. Neurosci. 9:26. doi: 10.3389/fncom.2015.00026

Received: 29 October 2014; Accepted: 09 February 2015;
Published online: 26 February 2015.

Edited by:

Julien Modolo, Western University, Canada

Reviewed by:

Jesus M. Cortes, Ikerbasque. Biocruces Health Research Institute, Spain
Chiun-Fan Chen, Spaulding Neuromodulation Center and Harvard Medical School, USA

Copyright © 2015 Azodi-Avval and Gharabaghi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Alireza Gharabaghi, Division of Functional and Restorative Neurosurgery and Division of Translational Neurosurgery, Department of Neurosurgery, Eberhard Karls University, Otfried-Mueller-Str.45, 72076 Tuebingen, Germany e-mail: alireza.gharabaghi@uni-tuebingen.de

ORIGINAL RESEARCH article

Phase-dependent modulation as a novel approach for therapeutic brain stimulation

Introduction

Materials and Methods

Phase Synchronization

Physiologically Informed Phase Response Curve

Results

Discussion

Conflict of Interest Statement

Acknowledgments

References

This article is part of the Research Topic

People also looked at