skip to main content


Title: Dynamic optimization of stimulation frequency to reduce isometric muscle fatigue using a modified Hill‐Huxley model
Abstract Introduction

Optimal frequency modulation during functional electrical stimulation (FES) may minimize or delay the onset of FES‐induced muscle fatigue.

Methods

An offline dynamic optimization method, constrained to a modified Hill‐Huxley model, was used to determine the minimum number of pulses that would maintain a constant desired isometric contraction force.

Results

Six able‐bodied participants were recruited for the experiments, and their quadriceps muscles were stimulated while they sat on a leg extension machine. The force–time (F–T) integrals and peak forces after the pulse train was delivered were found to be statistically significantly greater than the force–time integrals and peak forces obtained after a constant frequency train was delivered.

Discussion

Experimental results indicated that the optimized pulse trains induced lower levels of muscle fatigue compared with constant frequency pulse trains. This could have a potential advantage over current FES methods that often choose a constant frequency stimulation train.Muscle Nerve57: 634–641, 2018

 
more » « less
NSF-PAR ID:
10041886
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Muscle & Nerve
Volume:
57
Issue:
4
ISSN:
0148-639X
Page Range / eLocation ID:
p. 634-641
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Functional electrical stimulation (FES) is commonly used for people with neurological conditions. As the muscle geometry changes (i.e., muscle lengthening/ shortening), the force induced by static electrode placement may also change. Experimental results indicate that muscle forces can be increased by spatially switching stimulation as the muscle geometry changes with joint angle. In this letter, an electric field is switched between multiple electrodes placed across the biceps brachii to track a desired trajectory. A switched systems approach is used to develop a position-based switching law, including a switched robust sliding mode controller that successfully tracks the desired angular trajectory about the elbow, despite changes in muscle geometry. Lyapunovbased methods for switched systems are used to prove global exponential tracking. Experimental results from nine able-bodied subjects are presented and the developed control system achieves an average position and velocity error of −0.21 ± 1.17 deg and −0.43 ± 5.38 deg/s, respectively, and, on average, reduces fatigue by 13.6%, as compared to traditional single-electrode methods, demonstrating the performance of the uncertain nonlinear switched control system. 
    more » « less
  2. Abstract

    Objective.Transcutaneous electrical stimulation of peripheral nerves is a common technique to assist or rehabilitate impaired muscle activation. However, conventional stimulation paradigms activate nerve fibers synchronously with action potentials time-locked with stimulation pulses. Such synchronous activation limits fine control of muscle force due to synchronized force twitches. Accordingly, we developed a subthreshold high-frequency stimulation waveform with the goal of activating axons asynchronously.Approach.We evaluated our waveform experimentally and through model simulations. During the experiment, we delivered continuous subthreshold pulses at frequencies of 16.67, 12.5, or 10 kHz transcutaneously to the median and ulnar nerves. We obtained high-density electromyographic (EMG) signals and fingertip forces to quantify the axonal activation patterns. We used a conventional 30 Hz stimulation waveform and the associated voluntary muscle activation for comparison. We modeled stimulation of biophysically realistic myelinated mammalian axons using a simplified volume conductor model to solve for extracellular electric potentials. We compared the firing properties under kHz and conventional 30 Hz stimulation.Main results.EMG activity evoked by kHz stimulation showed high entropy values similar to voluntary EMG activity, indicating asynchronous axon firing activity. In contrast, we observed low entropy values in EMG evoked by conventional 30 Hz stimulation. The muscle forces evoked by kHz stimulation also showed more stable force profiles across repeated trials compared with 30 Hz stimulation. Our simulation results provide direct evidence of asynchronous firing patterns across a population of axons in response to kHz frequency stimulation, while 30 Hz stimulation elicited synchronized time-locked responses across the population.Significance.We demonstrate that the continuous subthreshold high-frequency stimulation waveform can elicit asynchronous axon firing patterns, which can lead to finer control of muscle forces.

     
    more » « less
  3. Abstract

    Despite its evident importance to learning theory and models, the manner in which the lateral perforant path (LPP) transforms signals from entorhinal cortex to hippocampus is not well understood. The present studies measured synaptic responses in the dentate gyrus (DG) of adult mouse hippocampal slices during different patterns of LPP stimulation. Theta (5 Hz) stimulation produced a modest within‐train facilitation that was markedly enhanced at the level of DG output. Gamma (50 Hz) activation resulted in a singular pattern with initial synaptic facilitation being followed by a progressively greater depression. DG output was absent after only two pulses. Reducing release probability with low extracellular calcium instated frequency facilitation to gamma stimulation while long‐term potentiation, which increases release by LPP terminals, enhanced within‐train depression. Relatedly, per terminal concentrations of VGLUT2, a vesicular glutamate transporter associated with high release probability, were much greater in the LPP than in CA3–CA1 connections. Attempts to circumvent the potent gamma filter using a series of short (three‐pulse) 50 Hz trains spaced by 200 ms were only partially successful: composite responses were substantially reduced after the first burst, an effect opposite to that recorded in field CA1. The interaction between bursts was surprisingly persistent (>1.0 s). Low calcium improved throughput during theta/gamma activation but buffering of postsynaptic calcium did not. In all, presynaptic specializations relating to release probability produce an unusual but potent type of frequency filtering in the LPP. Patterned burst input engages a different type of filter with substrates that are also likely to be located presynaptically.image

    Key points

    The lateral perforant path (LPP)–dentate gyrus (DG) synapse operates as a low‐pass filter, where responses to a train of 50 Hz, γ frequency activation are greatly suppressed.

    Activation with brief bursts of γ frequency information engages a secondary filter that persists for prolonged periods (lasting seconds).

    Both forms of LPP frequency filtering are influenced by presynaptic, as opposed to postsynaptic, processes; this contrasts with other hippocampal synapses.

    LPP frequency filtering is modified by the unique presynaptic long‐term potentiation at this synapse.

    Computational simulations indicate that presynaptic factors associated with release probability and vesicle recycling may underlie the potent LPP–DG frequency filtering.

     
    more » « less
  4. Introduction:Current brain-computer interfaces (BCIs) primarily rely on visual feedback. However, visual feedback may not be sufficient for applications such as movement restoration, where somatosensory feedback plays a crucial role. For electrocorticography (ECoG)-based BCIs, somatosensory feedback can be elicited by cortical surface electro-stimulation [1]. However, simultaneous cortical stimulation and recording is challenging due to stimulation artifacts. Depending on the orientation of stimulating electrodes, their distance to the recording site, and the stimulation intensity, these artifacts may overwhelm the neural signals of interest and saturate the recording bioamplifiers, making it impossible to recover the underlying information [2]. To understand how these factors affect artifact propagation, we performed a preliminary characterization of ECoG signals during cortical stimulation.Materials/Methods/ResultsECoG electrodes were implanted in a 39-year old epilepsy patient as shown in Fig. 1. Pairs of adjacent electrodes were stimulated as a part of language cortical mapping. For each stimulating pair, a charge-balanced biphasic square pulse train of current at 50 Hz was delivered for five seconds at 2, 4, 6, 8 and 10 mA. ECoG signals were recorded at 512 Hz. The signals were then high-pass filtered (≥1.5 Hz, zero phase), and the 5-second stimulation epochs were segmented. Within each epoch, artifact-induced peaks were detected for each electrode, except the stimulating pair, where signals were clipped due to amplifier saturation. These peaks were phase-locked across electrodes and were 20 ms apart, thus matching the pulse train frequency. The response was characterized by calculating the median peak within the 5-second epochs. Fig. 1 shows a representative response of the right temporal grid (RTG), with the stimulation channel at RTG electrodes 14 and 15. It also shows a hypothetical amplifier saturation contour of an implantable, bi-directional, ECoG-based BCI prototype [2], assuming the supply voltage of 2.2 V and a gain of 66 dB. Finally, we quantify the worstcase scenario by calculating the largest distance between the saturation contour and the midpoint of each stimulating channel.Discussion:Our results indicate that artifact propagation follows a dipole potential distribution with the extent of the saturation region (the interior of the white contour) proportional to the stimulation amplitude. In general, the artifacts propagated farthest when a 10 mA current was applied with the saturation regions extending from 17 to 32 mm away from the midpoint of the dipole. Consistent with the electric dipole model, this maximum spread happened along the direction of the dipole moment. An exception occurred at stimulation channel RTG11-16, for which an additional saturation contour emerged away from the dipole contour (not shown), extending the saturation region to 41 mm. Also, the worst-case scenario was observed at 6 mA stimulation amplitude. This departure could be a sign of a nonlinear, switch-like behavior, wherein additional conduction pathways could become engaged in response to sufficiently high stimulation.Significance:While ECoG stimulation is routinely performed in the clinical setting, quantitative studies of the resulting signals are lacking. Our preliminary study demonstrates that stimulation artifacts largely obey dipole distributions, suggesting that the dipole model could be used to predict artifact propagation. Further studies are necessary to ascertain whether these results hold across other subjects and combinations of stimulation/recording grids. Once completed, these studies will reveal practical design constraints for future implantable bi-directional ECoG-based BCIs. These include parameters such as the distances between and relative orientations of the stimulating and recording electrodes, the choice of the stimulating electrodes, the optimal placement of the reference electrode, and the maximum stimulation amplitude. These findings would also have important implications for the design of custom, low-power bioamplifiers for implantable bi-directional ECoG-based BCIs.References:[1] Hiremath, S. V., et al. "Human perception of electrical stimulation on the surface of somatosensory cortex." PloS one 12.5 (2017): e0176020.[2] Rouse, A. G., et al. "A chronic generalized bi-directional brain-machine interface." Journal of Neural Engineering 8.3 (2011): 036018 
    more » « less
  5. Functional electrical stimulation (FES) is a potential neurorehabilitative intervention to enable functional movements in persons with neurological conditions that cause mobility impairments. However, the quick onset of muscle fatigue during FES is a significant challenge for sustaining the desired functional movements for more extended periods. Therefore, a considerable interest still exists in the development of sensing techniques that reliably measure FES-induced muscle fatigue. This study proposes to use ultrasound (US) imaging-derived echogenicity signal as an indicator of FES-induced muscle fatigue. We hypothesized that the US-derived echogenicity signal is sensitive to FES-induced muscle fatigue under isometric and dynamic muscle contraction conditions. Eight non-disabled participants participated in the experiments, where FES electrodes were applied on their tibialis anterior (TA) muscles. During a fatigue protocol under either isometric and dynamic ankle dorsiflexion conditions, we synchronously collected the isometric dorsiflexion torque or dynamic dorsiflexion angle on the ankle joint, US echogenicity signals from TA muscle, and the applied stimulation intensity. The experimental results showed an exponential reduction in the US echogenicity relative change (ERC) as the fatigue progressed under the isometric (R2=0.891±0.081) and dynamic (R2=0.858±0.065) conditions. The experimental results also implied a strong linear relationship between US ERC and TA muscle fatigue benchmark (dorsiflexion torque or angle amplitude), with R2 values of 0.840±0.054 and 0.794±0.065 under isometric and dynamic conditions, respectively. The findings in this study indicate that the US echogenicity signal is a computationally efficient signal that strongly represents FES-induced muscle fatigue. Its potential real-time implementation to detect fatigue can facilitate an FES closed-loop controller design that considers the FES-induced muscle fatigue. 
    more » « less