METHOD AND SYSTEM FOR OPTIMISATION OF DBS PROGRAMMING

Abstract

A method and system are described for, based upon a plurality of previously-acquired directional LFP signals measured in a plurality of different directions at a directional sensor lead located in a predetermined region of a patient's brain, determining optimised patient-specific programming parameters for programming a directional stimulation lead with parameters for stimulating the said region. The method comprises a first step of determining, over at least one predetermined frequency range, a power-frequency variation curve of each of the directional LFP signals, a second step of identifying frequency peaks in the power-frequency variation curves, a third step of detecting one of the identified frequency peaks at which a maximum difference in signal power between the directional LFP signals occurs, and a fourth step of calculating a plurality of directional stimulation weighting factors on the basis of the relative signal powers of the directional LFP signals at the detected frequency peak.

Claims

1. A computer-implemented method of, based upon a plurality of directional LFP signals measured in a plurality of different directions at a directional sensor located in a predetermined region of a patient's brain, determining optimized patient-specific programming parameters for programming a directional multi-electrode stimulation lead with parameters for stimulating the region, the method comprising the steps of: determining, over at least one predetermined frequency range, a power-frequency variation curve for each of the directional LFP signals; identifying at least one frequency peak in the power-frequency variation curves; detecting one of the at least one identified frequency peaks at which a maximum difference in signal power between the directional LFP signals occurs; and calculating a plurality of directionally weighted stimulation indices for the directional stimulation lead on the basis of the relative signal powers of the directional LFP signals at the frequency of the detected frequency peak.

2. The method of claim 1, wherein the frequency range comprises a beta frequency range.

3. The method of claim 1, wherein the step of determining a power-frequency variation curve comprises determining power-frequency curves for a plurality of frequency ranges, the frequency ranges being usable alone or in combination to guide stimulation.

4. The method of claim 1, wherein the step of calculating a plurality of directionally weighted stimulation indices comprises: selecting a directionality indicating frequency range (DIFR), based on the detected frequency peak; and determining a directionality index (DI) for each of a plurality of the electrodes of the directional stimulation lead and/or for each of a plurality of montages of the electrodes of the directional stimulation lead.

5. The method of claim 4, wherein the DIFR is a predetermined constant frequency range centered on the frequency of the detected frequency peak, or wherein the width of the DIFR is selected in dependence on the frequency of the detected frequency peak.

6. The method of claim 4, wherein the step of determining the directionality index comprises performing a normalization and/or a baseline correction of an average amplitude of each of the power-frequency variation curves in the DIFR.

7. The method of claim 6, wherein the baseline correction comprises subtracting a predetermined proportion of the average amplitude from an amplitude of each of the curves.

8. The method of claim 7, further comprising summing the average amplitude of each curve over the DIFR to give a proportional distribution for each direction associated with each of the curves.

9. A system for determining optimized patient-specific programming parameters for programming a directional multi-electrode stimulation lead with parameters for stimulating a predetermined region of a patient's brain, the system comprising: means for receiving a plurality of previously-acquired directional LFP signals measured in a plurality of different directions at a directional sensor lead located in the predetermined region of the patient's brain; spectral analysis calculation means configured to determine, over at least one predetermined frequency range, a power-frequency variation curve for each of the directional LFP signals; peak detection means for identifying at least one frequency peak in the power-frequency variation curves determined by the spectral analysis calculation means; difference detection means configured for detecting one of the at least one identified frequency peaks at which a maximum difference in signal power between the directional LFP signals occurs; and directionality determining means for calculating a plurality of directionally weighted stimulation parameters for the directional stimulation lead on the basis of relative signal powers of the directional UP signals at a frequency of the detected frequency peak.

10. The system of claim 9, further comprising an implantable pulse generator device, wherein one or more of the spectral analysis calculation means, the peak detection means, the difference detection means, and the directionality determining means are integrated into the pulse generator device for driving the directional stimulation lead.

11. The system of claim 10, further comprising a sensor lead for directional LFP recording.

12. The system of claim 9, further comprising a stimulation lead for performing directional stimulation based on the directionally weighted stimulation parameters.

13. The system of claim 11, wherein the sensor lead and the stimulation lead are the same device.

14. The system of claim 9, further comprising a programmable implantable pulse generator device for driving the stimulation lead, wherein the programmable implantable pulse generator device comprises means for acquiring or recording LFP signal data from the sensor lead.

15. The method of claim 1, further comprising the step of non-transitory computer-readable medium storing processor executable instructions on a computing device for carrying out the method.

Description

DETAILED DESCRIPTION OF THE INVENTION

[0011] The invention and its advantages will be explained in greater detail with reference to the accompanying drawings, in which:

[0012] FIG. 1 shows a simplified schematic illustration of an example of a directional multi-contact DBS lead arrangement which may be used to provide LFP signal data to the method and system of the invention.

[0013] FIG. 2 shows an example of frequency-power spectra of LFP signal data acquired from the lead of FIG. 1 for a particular DBS site of a particular patient.

[0014] FIGS. 3a and 3b show schematic graphical representations of DBS directional weightings which may be determined by the inventive method and system from the signal data of FIG. 2.

[0015] FIG. 4 shows a simplified top-level abstraction of an example of a clinical procedure incorporating the method of the invention.

[0016] FIG. 5 shows a more detailed view of the steps carried out in a method according to the invention.

[0017] FIG. 6 shows in even greater detail the steps and substeps which may be carried out in a method according to the invention.

[0018] The drawings are intended merely as exemplary illustrations, for the purpose of understanding certain principles underlying the invention, and are not to be construed as limiting the scope of the invention. Where the same reference numerals are used in different drawings, these reference numerals are intended to refer to the same or corresponding features. However, the use of different reference numerals should not necessarily be taken as an indication that the referenced features differ in any particular respect.

[0019] To identify the target cerebral structure for chronic stimulation, an intraoperative recording of the neuronal activity may be performed, for example using a directional multi-contact lead such as the one shown in FIG. 1. Such a recording may reveal evidence of exaggerated oscillatory synchronisation in the STN, GPi or thalamus, which may be detected in the local field potentials (LFP). Two broad types of oscillatory activity in patients with PD have been identified: beta band activity (13-30 Hz) and gamma band activity (65-80 Hz). The power of the beta band activity without medical treatment (ie levodopa) and the degree of suppression of this activity by levodopa and DBS correlate with clinical signs of Parkinsonism (bradykinesia and rigidity) and the degree of clinical improvement with treatment, respectively. Gamma band activity, by contrast, has been shown to increase during improvement of bradykinesia and rigidity. The method and system of the invention make use of the fact that LFP signals in specific frequency bands, especially the beta band activity, are linked to the motor signs and may be used as biomarkers to optimise deep brain stimulation. LFPs recorded from directional DBS leads are screened for individual direction and depth-specific differences in the power frequency spectrum in disease-related frequency bands (eg beta band). A specific algorithm is then used to process these signal data in order to output parameters for the optimum contact configuration for programming the subsequent DBS stimulation regime. The method and system of the invention may be employed for automatic programming of the IPG connected to a directional DBS lead. Alternatively, the method and system of the invention may be employed to provide optimised programming parameters to a clinician, thereby providing a supportive tool for the clinician for programming of the IPG for the directional DBS lead.

[0020] FIG. 2 shows an example set of frequency-power spectra of LPF signal data acquired from the lead 1 of FIG. 1, for a particular DBS site of a particular patient. Curves 12, 13, 14, 15, 16, and 17 show the spectra of the LFP signal data acquired in the beta band from contacts 3.sub.1, 3.sub.2, 3.sub.3, 4.sub.1, 4.sub.2 and 4.sub.3 respectively. A local peak 10 can be seen in the curve 14 for contact 4.sub.3 at approximately 20 Hz. As will be described below, this peak 10 can be detected using an algorithm for detecting peaks in each spectral power curve and determining at which of the peaks the difference between the amplitudes (ie the spread of the amplitude values) among one set of directional contacts (ie 3.sub.1, 3.sub.2 and 3.sub.3 or 4.sub.1, 4.sub.2 and 4.sub.3) are greatest. As will be described below, this maximum difference calculation may be based on normalized and/or baseline-corrected amplitude values of the frequency curves within a frequency window (referred to as the DIFR) at the peak under consideration. Signal differences and signal comparability can be improved by carrying out normalization and/or baseline correction steps in order to increase the signal-to-noise ratio. Normalization steps may include various spectral components, and may involve for example contrasting the the peak amplitude with the mean amplitude of the whole disease-related frequency band.

[0021] Note that the example given in this description assumes that the geometrical arrangement of the sensor lead contacts corresponds directly to that of the stimulation lead, however this need not necessarily be the case. Stimulation weightings may be mapped using appropriate transformation from the geometric configuration of the directional sensor(s) to the geometric configuration of the directional stimulation electrode(s).

[0022] Similarly, the method steps are described in the context of monopolar measurements (ie assessing each contact separately, with a common reference). In the monopolar method, power frequency curves are derived for each contact, and the magnitude of the disease-related spectral component is ranked for the contacts individually. However, it should be understood that the same techniques may be applied to more than single contact at once, in a bipolar or multipolar fashion, such that so-called montages or arrangements of multiple contacts may be assessed and ranked. Method steps applied to individual contacts in this description should be understood to encompass the application of the steps to montages of two or more contacts. Determining the magnitude of the disease-related spectral components in such a group-wise, combinatorial fashion greatly increases the number of possible choices to be ranked, and improves the signal-to-noise ratio of the derived results (power frequency curves and detected peaks).

[0023] FIGS. 3a and 3b show how the spectral power curves of FIG. 2 may be used to derive directional weightings 6 and 7 for the stimulation amplitudes at contacts 3.sub.1, 3.sub.2 and 3.sub.3 (FIG. 3a) and contacts 4.sub.1, 4.sub.2 and 4.sub.3 (FIG. 3b). In FIG. 3a, the three electrodes 3.sub.1, 3.sub.2 and 3.sub.3 show almost zero directional weighting (8%, 9% and 8% respectively), whereas FIG. 3b shows that the three electrodes 4.sub.1, 4.sub.2 and 4.sub.3 are markedly weighted (16%, 23% and 35% respectively). In this example, the inventive method predicts that the DBS should be weighted in favour of the contact 4.sub.3 in order to achieve optimum targeting of the structure to be stimulated. A manual contact testing was carried out blind on the patient whose LFP signals were recorded to generate the spectral curves of FIG. 2, and the manual contact testing produced the following recommended directional stimulation currents for each of the directional contacts of the DBS lead of FIG. 1 as follows:

[0024] 3.sub.1: 4.0 mA

[0025] 3.sub.2: 3.5 mA

[0026] 3.sub.3: 3.0 mA

[0027] 4.sub.1: 2.5 mA

[0028] 4.sub.2: 2.5 mA

[0029] 4.sub.3: 2.0 mA

[0030] As can be seen from the above, the manual contact testing approach confirmed the prediction, made by the method of the invention, that contact 4.sub.3 would provide the best result (lowest stimulation current for effective symptom relief). However, the manual contact testing took six hours of intensive, fatiguing clinical work, while the recommendation from the inventive method was delivered almost instantly.

[0031] Note that the graphical representation of FIGS. 3a and 3b are examples of how the directional contact weighting may be presented graphically to a clinician so that he or she may adjust the DBS programming accordingly. Alternatively, the calculated weightings may be used directly in an automated programming of the IPG and DBS electrode.

[0032] FIG. 4 shows an overview of the context of the method and system of the invention in the overall DBS process. In step 20, directional LFP signals are acquired from directional sensing electrode 1, such as the lead shown in FIG. 1. The acquired signals are processed in step 21, using to a method according to the invention, to generate directionally optimised programming parameters which are then used in step 22 to program an IPG driving a directional DBS lead 8 (which may typically be the lead 1, ie the DBS lead serves as sensor and stimulation lead).

[0033] The LFP signal acquisition step 20 may comprise recording LFPs from the directional lead after its placement in the definitive position within the target structure (eg STN, GPi or thalamus). During the recording, the patient must be withdrawn from dopaminergic medication, the patient must be awake and in a resting position without any voluntary movement. Recording can be performed intraoperatively or postoperatively before the electrode is connected to the implantable pulse generator (IPG). Alternatively, LFPs can be recorded at any time point from the IPG itself, if the IPG is capable of LFP recording.

[0034] The method steps underlying step 21 of FIG. 4 are described below, with particular reference to FIGS. 5 and 6. As described above, the method of the invention comprises the steps of [0035] spectral analysis 30 of the LFP signal data in order to generate spectral power curves for each directional electrode on the sensor lead 1; [0036] peak detection 31 to identify any local peaks in each of the spectral power curves; [0037] difference detection 32 for determining the peak at which there is the greatest difference in amplitude among the directional electrodes of a section of the lead; and [0038] directionality calculation 33, for deriving weightings for stimulation signals applied to each of the directional electrodes of the DBS lead 8.

[0039] FIG. 6 shows an example flow diagram of a method according to the invention. In step 40, one or more suitable frequency bands are selected, according to current knowledge, for analysing the LFP signals associated with a particular disease. The frequency band(s) may advantageously comprise the beta frequency range (13-35 Hz). In principle, several frequency bands present in the LFP may be used as biomarkers (alone or in combination) to guide stimulation, however on the premise of an existing relationship to the disease. Such a relationship could be positive (associated with symptom alleviation) or negative (associated with symptom deterioration). In step 41, a suitable sampling frequency is determined for the subsequent spectral analysis, the sampling frequency being preferably at least twice the Nyquist limit for the LFP signal data. In step 42, the power-frequency-spectrum (spectral analysis with frequency decomposition) for each directional contact is calculated, at eg 1 Hz resolution. In step 43, local peaks in the power-frequency curves are identified. Within the disease-related frequency band, the frequency peak showing the largest differences among directional contacts is identified in step 44. In step 45, this peak frequency is transformed into a frequency range (DIFR) by adding and subtracting a predetermined spectral bandwidth, such as 3 Hz, to give a total width of 7 Hz. The width of DIFR (eg peak frequency 3 Hz) can be set to maximize directional information between the contacts. The width of the DFIR may be predetermined, or it may be dynamically adapted, for example in dependence on a characteristic of the power-frequency curves. In step 46, the average amplitude for each directional contact over the DIFR is baseline corrected by subtracting 90% of the mean amplitude of the directional contact with the lowest amplitude in the DFIR. This step helps to reduce the effect of volume conduction and to improve the visualisation of the spectral differences between the contacts. The baseline corrected average amplitudes for each contact are then summed in step 47, and percentage distribution over the directional contacts is calculated to determine a directionality index (DI). The DI of each directional contact allows contacts to be ranked according to the likelihood that they will afford the best clinical response on stimulation, which is defined as the lowest stimulation intensity required for sufficient clinical response. The percentage value itself does not indicate a percentage clinical response. Where two or more of the directional contacts in a particular lead section have similar DI values, they may be considered in combination for stimulation. The directionality index is superimposed graphically as a vector and as percentage value for each directional contact on a visual model illustrating the stimulation contacts and contact levels. This enables the clinician to immediately understand which section and which directional contact are likely to produce the best DBS results.

[0040] The parameters determined in the method described above may be exported (manually or automatically) into the IPG which drives the DBS stimulation device. The clinician should therefore start to deliver the stimulation on the contact with the highest directionality index, adjusting other stimulation parameters (current, frequency, pulse width etc) accordingly and, if necessary, move on to a different contact, or a different combination of contacts, suggested by the directionality index results if the stimulation effect is clinically not sufficient or if side effects occur. In IPGs with LFP recording capabilities and integrated analysis modules, the inventive method, and the system which embodies it could be fully integrated and automated as an internal feature of the IPG.

[0041] The method described above may preferably be implemented as instructions stored on non-transitory computer-readable media, and/or in a system comprising one or more specially configured or programmed electronic circuits.