A neuroprosthesis for tremor management through the control of muscle co-contraction
© Gallego et al.; licensee BioMed Central Ltd. 2013
Received: 24 April 2012
Accepted: 25 March 2013
Published: 15 April 2013
Skip to main content
© Gallego et al.; licensee BioMed Central Ltd. 2013
Received: 24 April 2012
Accepted: 25 March 2013
Published: 15 April 2013
Pathological tremor is the most prevalent movement disorder. Current treatments do not attain a significant tremor reduction in a large proportion of patients, which makes tremor a major cause of loss of quality of life. For instance, according to some estimates, 65% of those suffering from upper limb tremor report serious difficulties during daily living. Therefore, novel forms for tremor management are required. Since muscles intrinsically behave as a low pass filter, and tremor frequency is above that of volitional movements, the authors envisioned the exploitation of these properties as a means of developing a novel treatment alternative. This treatment would rely on muscle co-contraction for tremor management, similarly to the strategy employed by the intact central nervous system to stabilize a limb during certain tasks.
We implemented a neuroprosthesis that regulated the level of muscle co-contraction by injecting current at a pair of antagonists through transcutaneous neurostimulation. Co-contraction was adapted to the instantaneous parameters of tremor, which were estimated from the raw recordings of a pair of solid state gyroscopes with a purposely designed adaptive algorithm. For the experimental validation, we enrolled six patients suffering from parkinsonian or essential tremor of different severity, and evaluated the effect of the neuroprosthesis during standard tasks employed for neurological examination.
The neuroprosthesis attained significant attenuation of tremor (p<0.001), and reduced its amplitude up to a 52.33±25.48%. Furthermore, it alleviated both essential and parkinsonian tremor in spite of their different etiology and symptomatology. Tremor severity was not a limiting factor on the performance of the neuroprosthesis, although there was a subtle trend towards larger attenuation of more severe tremors. Tremor frequency was not altered during neurostimulation, as expected from the central origin of Parkinson’s disease and essential tremor. All patients showed a good tolerance to neurostimulation in terms of comfort and absence of pain, and some spontaneously reported that they felt that tremor was reduced when the neuroprosthesis was activated.
The results presented herein demonstrate that the neuroprosthesis provides systematic attenuation of the two major types of tremor, irrespectively from their severity. This study sets the basis for the validation of the neuroprosthesis as an alternative, non-invasive means for tremor management.
Pathological tremors, understood as the ensemble of tremors that are cause of functional disability, constitute the most prevalent movement disorder. Although obtaining an estimate of its prevalence is challenging , a recent population study showed that tremor affects 15% of those with age ranging between 50 and 89 years . According to some estimates, more than 65% of those who suffer from pathological tremor—referred to as tremor in the remainder of the paper—at the upper limb report serious difficulties when performing their activities of daily living (ADLs) . Therefore, tremor has an enormous impact on the independence and quality of life of many.
The two most relevant conditions, in terms of prevalence, that cause tremor are essential tremor (ET) and Parkinson’s disease (PD). ET is believed to originate by abnormal neuronal oscillations at the cerebellar and thalamocortical pathways [4, 5] and, according to some, constitutes the most common movement disorder in adults  (it affects ∼ 5% of people over age 65 ). Much is yet unknown about ET, for example there is still no consensus on whether it is a neurodegenerative or non-degenerative disease . PD is a neurodegenerative disease that arises from the death of dopaminergic nigrostriatal neurons, causing abnormal oscillations in the loop linking the cortex, basal ganglia and thalamus . PD is estimated to afflict more than ten million people all over the world . In addition to them, up to eight more syndromes are acknowledged to cause tremor . Importantly, no tremor is yet fully understood , which hampers the development of novel therapies, and the refinement of the existing ones.
Tremors are managed either through pharmacotherapy or neurosurgery, the latter consisting in the implantation of a deep brain stimulator (DBS), or the realization of therapeutical lesions in deep brain structures (gamma knife thalamotomy). Nonetheless, both pharmacotherapy and neurosurgery carry drawbacks associated. Drugs often induce side effects [11, 12], and show decreased effectiveness over years of use [13, 14], while DBS is related to increased risk of intracranial haemorrhage  and psychiatric manifestations , and the percentage of eligible patients is extremely low . As for the tremors themselves, the mechanisms accounting for the alleviation of their symptoms by these therapies are not fully elucidated.
Among the alternative approaches for tremor management, the application of mechanical loads, either forces or masses, is regarded as one of the most promising solutions. A number of works show how adequate inertial (mass) loading of the tremulous limb reduces the amplitude of most types of tremor [18, 19], and so does the application of forces with certain characteristics. Force loads typically consist of viscous fields [20, 21], although the effect of added stiffness  and inertia  on tremor has also been investigated. Furthermore, voluntary [24, 25] or artificially elicited  muscle contractions affect the severity of the tremor, by generating forces on the musculoskeletal system.
These findings motivated the development of orthotic devices for tremor attenuation, being the majority of them wheel-chair mounted [27, 28] or fixed to an external frame [29, 30]. To the best of our knowledge, the only ambulatory system for tremor attenuation through mechanical loading extensively validated with patients is the wearable robot WOTAS . This system attained consistent attenuation of tremors over a certain severity but, however, did not fulfill users’ expectations in terms of cosmetics and aesthetics . Besides, it was found that the intrinsic difficulty of transmitting external forces to the skeletal system through the soft tissues hampered the improvement of solutions of this kind. Therefore, here we present a novel approach for tremor management aimed at circumventing these limitations. Our novel solution consists of a neuroprosthesis (NP) that applies forces to the tremulous limb through transcutaneous neurostimulation. We show that by co-contracting the affected muscles it is possible to alter the inherent low pass filter properties of human muscles  and attenuate the tremor, without affecting the concomitant voluntary movement. Our results demonstrate that the NP provides systematic alleviation of tremors of different etiology and characteristics, avoiding discomfort and pain. Notice that transcutaneous neurostimulation for tremor management had been preliminarily evaluated in . There, an alternanting stimulation pattern that opposed the tremor activation bursts was employed to reduce its severity . The system was, however, table mounted.
This paper is organized as follows. First, we describe the NP and the controller it implements, together with the experimental protocol and methods employed for its validation. Afterwards, we review and discuss the results achieved. The manuscript ends up with conclusions summarizing the outcomes of the study.
The rationale behind the use of co-contraction was to increase simultaneously the stiffness and viscosity of the tremulous limb, which in turn would decrease the cut-off frequency (∼ 2–3 Hz [31, 32, 34]) of the inherent low pass filter response of muscles, filtering out the tremor. Indeed, co-contraction is exploited by the intact central nervous system to stabilize the limbs during specific tasks [35, 36]. Furthermore, this could be implemented directly by increasing muscle activation, because stiffness and viscosity are monotonic functions of the contraction level .
where θ(s) is the tremulous component of movement (estimated as described in epigraph ‘Tremor Parameterization’), T(s) represents the torque that generates the tremor, I, D, and K stand for the inertia, viscosity and stiffness of the joint respectively, and G is the magnitude of the response of the resultant system (that is function of its mechanical parameters). The viscosity and stiffness added by the NP as a result of increased muscle contraction are denoted by DNP and KNP, and satisfy DNP,KNP≥0, being DNP,KNP=0 if the NP was not activated.
In our case, the level of artificially elicited co-contraction was adapted to the characteristics—instantaneous amplitude and frequency—of the tremor. An independent controller was implemented for each joint, because the characteristics of tremor differs among them, and show a non-stationary behavior that is normally uncoupled. Neurostimulation was controlled as follows. Tremor frequency defined when the control output was to be updated, while tremor amplitude modulated the amount of current to be injected to each muscle. Tremor parameters estimated during a (tremor) period were employed to generate the control action in the subsequent one, analogously to repetitive control theory . This provided the controller with a certain predictive nature. The reason to modulate the level of neurostimulation independently for each muscle in the pair of antagonists was twofold: first, because of their different electrophysiological response to neurostimulation, and second, because this response also varies with time in a different manner for each muscle. Neurostimulation was modulated in amplitude; current frequency and pulse width were kept constant. Details on the adaptive filter employed to estimate instantaneous tremor amplitude and frequency, and on the control algorithm, are given next.
where y t r e m o r (k) represents the tremor, and y v o l (k) the component of voluntary movement. A two–stage filter , summarized next, separated the voluntary movement and the tremor in order to estimate the instantaneous parameters of the latter. This information was in turn employed by the controller to drive the NP.
By selecting an adequate value of parameter θ we obtained a precise estimation of voluntary movement with negligible delay, which is of foremost importance to get an accurate tracking of tremor parameters. The delay introduced by Butterworth, Chebyshev and elliptic low pass filters was notably higher, which motivated our choice.
Instantaneous tremor frequency and amplitude were subsequently estimated (in the second stage) with a Weighted Frequency Fourier Linear Combiner (WFLC)  and a Kalman Filter  in the following manner: the WFLC tracked tremor frequency, and this value was fed into the Kalman filter, which estimated the instantaneous amplitude of tremor. This cascade architecture permited optimizing the tracking of both parameters .
where A and B are the amplitude terms of the Fourier series, and ω t is the instantaneous frequency derived with the WFLC; t rk,k−1 represents the estimation of tremor based on these three parameters. Measurement noise covariance was defined as , whereas process noise Q(k) was modeled as a diagonal matrix, assuming independence of the state variables. The measurement noise covariance was related to the accuracy of the estimation of yt r e m o r,k with the g−h filter, while the process noise covariance determined the adaptation rate of the amplitude terms of the truncated Fourier series, A and B.
where K p and K i (t rk,k−1), are the controller gains (the value of the latter was switched depending on tremor amplitude, as shown in (11)), and u(k) the control action. Importantly, a saturation was applied to u(k) in order to limit the electrical charge injected to the muscle. This value was defined for each muscle during a calibration phase.
where and represent the integral of the amplitude spectrum in the 0–3 Hz and 3–12 Hz bands. Similar ratios have been employed in other works about tremor, e.g. for EMG analysis .
Details of patients
Tremor frequency (Hz)
Patients wore the NP at the most affected side and were asked to perform the clinical task that made their tremor more evident. Typical exercises for neurological examination were considered: resting the arm on the lap, keeping both arms outstretched, the finger to finger test, and the finger to nose test . Patients were sitting in a comfortable armchair during the whole session. The NP recorded both wrist flexion/extension and elbow flexion/extension with two pairs of solid state gyroscopes (Technaid S.L., Madrid, Spain), although we present results only for the wrist, because it was where the tremor was visible for all patients. Neurostimulation was delivered at the flexor carpi ulnaris and extensor carpi radialis with a multichannel monopolar neurostimulator with charge compensated pulses (UNA Systems, Belgrade, Serbia); the common electrode was located at either the dorsal or volar side of the wrist. Maximum neurostimulation amplitude for each muscle was personalized during an initial calibration phase; pulse width and frequency were set to either 300 μ s and 40 pps or 250 μ s and 30 pps respectively, depending on the patient: we used the first combination of values in patients that needed high current density to elicit a visible muscle contraction. The controller was implemented in a stand alone computer (QNX Software Systems, Ontario, Canada). Figure 1(B) shows a patient wearing the NP.
Each patient performed a number of repetitions of two types of trials; the duration of each was 30 s. The first type of trial consisted of two 15 s sub-periods, during the second of which the NP was activated. In the second type of trial, the NP was never activated. This type of trial was included in the protocol to avoid a possible distortion of the study due to, for example, a natural reduction of tremor amplitude over time. Both types were randomized, and the experimental design was balanced using an optimal algorithm, in which the repetitions were ordered by concatenating latin squares . In total, each patient performed between 6 and 12 repetitions. Patients with mild or moderate tremor were asked to count mentally backwards during the experiments to exacerbate their tremor [51, 52], and the two with PD counted out loud at the beginning of the session as recommended in . Patients with severe or very persistent tremor were not asked to count during the recordings to avoid excessive mental stress, which could have possible ethical implications.
To quantify the effect of the NP on tremor amplitude, we computed the ratio R a t t of the integral of the power spectral density of the tremor during the part of the trial with co-contraction, to the same variable without it . Before, data were split into 1 s non-overlapping windows to minimize the effects of eventual non-stationarities, and zero padded; R a t t was calculated with the median values for both conditions. We computed R a t t also for the trials during which the NP was not activated, in order to investigate whether it had a real effect on the tremor. To this end, we considered equivalent periods for its calculation.
We assessed whether the NP had a real effect on tremor amplitude by comparing R a t t in the ensemble of trials of the two types (pooled trials in which these NP was activated or not) with a Mann-Whitney test, given that the pooled datasets did not conform normality after transformation (Kolmogorov-Smirnov test, p<0.05). We also evaluated whether there was an effect on tremor frequency with a Mann-Whitney test. Throughout the paper, values are reported as mean ± SD.
The parameters for the tremor parameterization algorithms and the controller were set to: i) for the g−h filter, θ= 0.900, ii) for the Kalman filter, 0.01, Q(k)=d i a g(1,1,0), and iii) for the WFLC, μ0= 0.001, μ1= 0.01, μ b = 0.01, M= 1, f0 was initialized to the value computed from the amplitude spectrum at tremor onset (data epoch of 2 s with zero padding). The parameters of these three algorithms were obtained by analyzing offline a previously recorded dataset. Controller gains were selected manually based on the amplitude of the tremor that was observed during calibration and the required neurostimulation amplitude. During this phase, we also defined the saturation level (maximum current amplitude) the controller could deliver at each muscle. Thresholds for tremor onset, frequency estimation, integral reset, and gain switching were t r t h = 0.1 rad/s, T V R t h = 3, t h i n t r e s e t = 0.1 rad/s, and t h i n t g a i n = 0.1 rad/s.
These results demonstrate that the NP here presented constitutes a feasible approach to tremor management. Through the control of muscle co-contraction, the NP systematically alleviated the tremor independently from its characteristics, proving that such approach is a viable alternative for treating upper limb tremor. Interestingly, all patients reported that the sensation generated by the NP was tolerable and not unpleasant, and the overall impression was that they could habituate to it. Furthermore, a few patients spontaneously declared that when the NP was activated they could control better their limbs. Users’ perception of the prototype was generally good, and remarkably better than for previous robotic devices.
The NP attained a reduction of tremor amplitude in patients suffering both from PD and ET, in spite of their inherent differences in underlying pathophysiologic mechanisms  and symptomatology. Our initial concern was that joint rigidity arising from PD  could hinder tremor reduction with the NP in these patients, given that their joints already show increased stiffness and viscosity . The physical reason for this would be a displacement of the natural cut off frequency of their muscles towards lower values (for healthy individuals it is in the range ∼ 2–3 Hz [31, 32, 34]), which could in turn impede attaining the level of co-contraction necessary to attenuate the tremor. However, our results deemed this hypothesis untrue, even in a patient with severe parkinsonian tremor of low frequency (patient 02, who had a tremor frequency of ∼ 3.5 Hz and experienced a tremor reduction with the NP of R a t t =26.79±13.13%). We believe that the effects of levodopa intake might have facilitated this by decreasing limb rigidity. The fact that the NP successfully alleviated low frequency tremors is of great interest, as the frequency of both essential and parkinsonian tremor experience a decline with time , and thus broadens the group of potential users of the system.
In this line, tremor frequency did not to influence the performance of the NP, although our group was biased to patients suffering from low or medium frequency tremor (below ∼ 5 Hz, see Table 1). This characteristic was a consequence of our patient selection because: i) we included two patients suffering from PD, which tremor frequency is typically below 6 Hz , and ii) the severity of ET is inversely correlated with its frequency (the slope between the log displacement and the log amplitude is ∼−4 ), and our group comprised two subjects with moderate and one with severe tremor (patients 01, 03 and 05). Nevertheless, the performance of the NP for the patient with high frequency tremor (patient 06) was good, achieving an attenuation of R a t t =34.45±8.87%. Although this patient exhibited tremor which amplitude decayed over time (R a t t =77.35±29.37% when the NP was not active), the improvement with the NP was clearly noticeable when both conditions are compared. Furthermore, it was not clear a priori whether the NP could attenuate low frequency tremors due to the fact that they manifest close to the natural cut off frequency of muscles, and thus require from a larger alteration of their physiologic response. Given that average attenuation for patient 02 (tremor frequency ∼ 3.5 Hz) was R a t t =26.79±13.13%, and for patient 01 (tremor frequency ∼ 4.2 Hz) was R a t t =58.99±26.81%, we conclude that low frequency tremors can be successfully managed with the NP. These results suggest that tremor frequency is not a criterion that restricts the applicability of the system.
As shown in Figure 7, the NP did not alter tremor frequency in spite of the obvious proprioceptive feedback. This is motivated by the central origin of tremor in PD  and ET [4, 5], which predominates over the other mechanisms that contribute to their genesis, namely long and short latency reflexes, and mechanical factors . Therefore, our results are in line with evidences on the limited role of sensory feedback in the generation, maintenance and modulation of tremor in PD . ET patients however, exhibit a more evident interaction between the stretch reflex and the tremor itself , being the most noticeable example the separation of both components, otherwise entrained, under inertial loading . In this line, visual inspection of Figure 4(F) (bottom plot, the power spectral density) suggests that the central and reflex components were separated in patient 06, both when the NP was activated and not. A more profound investigation of this phenomenon is needed, and would require from the integration of electromyographic (EMG) recordings and advanced signal processing techniques—to remove the artifacts generated by the NP on the EMG signal—in a novel protocol.
Furthermore, given that the NP implements a proportional integral controller, our results suggest that an increase in tremor amplitude was compensated by an increase in the level of muscle co-contraction. As a consequence, we can conjecture that both essential and parkinsonian tremor are reduced when muscle contraction is sufficiently augmented, i.e. when the stiffness and viscosity of the tremulous joint are increased due to their monotonic relationship with the activation level . Comparing this with previous works that assessed the influence of voluntary muscle contraction on the amplitude of tremor, our results are in line with what has been reported for PD , but contradicts a recent study that showed how high contraction intensities caused larger amplitude fluctuations in ET . Remarkably, these studies assessed the effect of volitional muscle activation on tremor, while ours deals with artificial activation with a NP, and thus motor control mechanisms are involved in a lesser extent. Therefore, we believe that the pathophysiologic landmarks of PD and ET may account for the difference. Although the cerebellothalamic pathways play a role in the generation of parkinsonian symptoms , no evidence of cerebellar involvement apart from a compensative hyperactivation —increased as the disease progresses —, has been found. As for ET, on the contrary, all three cerebellar areas are impaired to a certain extent . Thus we believe that cerebellar disfunction could explain the exacerbation of ET during voluntary contractions, possibly through reticulospinal projections to muscles [61, 62].
In a few trials, we observed that tremor migrated towards proximal joints when it was suppressed at the more distal ones; i.e. when it was suppressed at the wrist it appeared or increased at the elbow or at the shoulder. This occurred in patients with both PD and ET, and has been previously reported in two studies dealing with inertial or force loading in tremor. The first was a study on coordination during postural holding in PD, where the authors found that mass loading of the index finger enhanced the tremor at the proximal segments while posture was maintained . The second focused on the attenuation of tremors of different etiology (mainly ET) with a wearable robot that applied a constant viscosity to the affected joints . None of them found a detailed physiologic explanation for tremor migration. Therefore, further research on this topic is needed, given that it constitutes a major aspect when developing NPs or neurorobots for functional compensation of tremor during daily living. Future studies will need to address both the coordination mechanisms per se, and the influence that inertial loads and muscle activation have on them.
When compared with current treatments for ET, the reduction of tremor provided by the NP (R a t t =50.37±29.01%) was similar to that of drugs with proved efficacy (∼50%), but worse than the results attained with DBS (tremor attenuation ∼ 50 – 80%, ). Nevertheless, the latter implies a surgical procedure, while the NP is envisioned as a mock up system that the user may wear underneath his clothes. Furthermore, when compared to pharmacotherapy, the NP reduced tremor amplitude in all 4 ET patients in our group, while drugs are effective in ∼ 50% of patients . No conclusive results can be drawn given the size of our group, although this finding is encouraging. The 2 patients with PD continued their intake of levodopa during the experiments but experienced a reduction of their tremor (R a t t =42.56±24.91%). This suggests that for these patients the NP was an interesting approach to complement pharmacotherapy.
Finally, the limitations of our study need to be acknowledged. First, although our results demonstrate that the NP constitutes a feasible alternative to tremor management, we believe they may be improved by developing an automatic method to select the gains of the controller. This would require the identification of the parameters of a model of muscle response to stimulation, which constitutes our future work. As mentioned above, we believe that better gain selection could have eliminated most trials with tremor exacerbation. Second, for extended use of the NP, it would be necessary to adapt the parameters of the controller when the muscle exhibits accommodation to the ongoing level of neurostimulation. Manual tuning of the parameters by a practitioner, together with adaptive control techniques may provide a solution for this.
This study presented the design, implementation and validation of a novel NP for tremor management. Experimental results with a representative group of patients (mean reduction of tremor amplitude R a t t =52.33±25.48%, p<0.001) proved the concept of tremor attenuation through modulation of muscle co-contraction with an ambulatory NP. Importantly, consistent tremor attenuation was achieved irrespectively from its severity, frequency, and etiology (validation was performed with patients suffering from ET or PD), although there was a trend towards greater reduction of more severe tremors. Importantly, the patients found the sensation induced by the NP not unpleasant, and informed that they could use it during daily living. Furthermore, some of them reported spontaneously that they felt that tremor was reduced when the NP was activated. All patients exhibited a positive response to the NP, which is of great interest given that a significant proportion of those suffering from tremor do not respond to medication (e.g., 50% of ET patients). These results encourage the functional and clinical evaluation of the NP as a non invasive alternative to tremor treatment, either alone, or as a complement to pharmacotherapy.
aThis rating only included those items related to tremor severity at the most affected limb.
Activity of daily living
Deep brain stimulator
Ratio of tremor to voluntary movement.
This work has been supported by the Commission of the European Union through grants ICT2007-224051 (TREMOR) and ICT-2011-287739 (NeuroTREMOR), by the Spanish Ministry of Science and Innovation through grants DPI2008-06772-C03-01 (REHABOT) and CSD2009-00067 (HYPER), and by the IMPIVA through grant IMDEEA/2011/91 (SISTR2). The authors gratefully acknowledge all the volunteers that participated in the study, and also people within Asociación Parkinson Valencia for their cooperation throughout the project. They also thank A. D. Koutsou (Bioengineering Group, Consejo Superior de Investigaciones Científicas), S. Mena (Physiotherapy Department, University of Valencia), and A. Castillo (Servicio de Neurología, Hospital Universitario de Valencia) for their help with the experiments.
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.