|
|
||||||||
REPORT
The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, Georgia
Submitted 5 October 2007; accepted in final form 16 December 2007
|
|
ABSTRACT |
|---|
|
|
|
INTRODUCTION |
|---|
|
Typically, feedback models of human postural control have reproduced joint torques and segmental motions of the body, but not muscle activity. Using single- or multilink inverted-pendulum models, they demonstrate that a set of time-invariant feedback gains can explain joint kinematics during either quiet standing or postural responses to perturbations (Alexandrov et al. 2001a
; Bortolami et al. 2003
; Kiemel et al. 2002
; Kuo 1995
; Park et al. 2004
; Peterka 2000
; Runge et al. 1995
; van der Kooij et al. 1999
). Because feedback loops at each joint are used to generate stabilizing joint torques, these models cannot uniquely specify temporal patterns of muscle activation. Muscles must be explicitly included because the low-pass dynamics of the body introduce redundancy in the temporal domain, whereby different temporal patterns of muscle activation can produce similar kinematic outputs (Gottlieb et al. 1995
; Lockhart and Ting 2007
).
Evidence suggests that muscle activity during human postural responses is dependent on acceleration, velocity, and displacement signals, as previously demonstrated in cats. In response to support-surface translations, temporal patterns of muscle activity in humans and cats have a similar rapid initial rise followed by a longer, sustained plateau region (Macpherson et al. 1989
). In cats, this waveform is due to CoM acceleration, velocity, and displacement feedback (Lockhart and Ting 2007
). Consistent with this feedback model, muscle activity in human postural responses have been shown to be modified by perturbation velocity and total excursion (Diener et al. 1988
), smoothness of the initial perturbation trajectory or acceleration (Brown et al. 2001
; Siegmund et al. 2002
; Szturm and Fallang 1998
), and the deceleration impulse at the end of the perturbation (Bothner and Jensen 2001
; Carpenter et al. 2005
; McIlroy and Maki 1994
).
We hypothesized that the activity of multiple muscles during human postural responses to perturbation is generated by a common delayed feedback law based on CoM motion. As a first step, we scaled the single inverted-pendulum feedback model used in Lockhart and Ting (2007)
to human dimensions (similar to Peterka 2000
) and examined whether this model was capable of reconstructing temporal patterns of muscle activation in proximal and distal muscles. We examined forward and backward support-surface perturbations to standing balance that elicited "ankle strategy" responses (Horak and Nashner 1986
). We demonstrate that a delayed feedback law on CoM acceleration, velocity, and displacement can reconstruct temporal patterns of both muscle activity and CoM kinematics during postural responses to support surface translations.
|
|
METHODS |
|---|
|
Platform acceleration and position, and surface EMG from 11 muscles in the legs and trunk were collected at 1,080 Hz, synchronized with body-segment kinematics collected at 120 Hz (Fig. 1A). Platform signals were low-pass filtered at 30 Hz (third-order zero-lag Butterworth filter). EMGs were collected from the following muscles on the right side of the body: TA, tibialis anterior; MG, medial gastrocnemius; SOL, soleus; VLAT, vastus lateralis; RFEM, rectus femoris; SEMB, semimembranosus; SEMT, semitendinosus; BFLH, long head of biceps femoris; BFSH, short head of biceps femoris; ES, erector spinae; and RA, rectus abdominis. Raw EMG signals were high-pass filtered at 35 Hz (third-order zero-lag Butterworth filter), demeaned, half-wave rectified, and low-pass filtered at 40 Hz (first-order zero-lag Butterworth filter). EMG signals were then normalized to the maximum EMG observed in each muscle over all conditions for each subject. Body-segment kinematics were derived from a custom bilateral Helen–Hay 25-marker set that included head–arms–trunk (HAT), thigh, and shank-foot segments. Center-of-mass motion was calculated from kinematic data as a weighted sum of segmental masses (Winter 2005
).
|
We determined whether our feedback model could reproduce the time course of EMG signals in each subject. The model consisted of a single-link inverted pendulum, with a point mass m (equivalent to each subject's mass) and length h (equal to the height of each subject's CoM during quiet standing) (Fig. 1B). Disturbance torques calculated from experimentally recorded platform accelerations were applied at the ankle to model the effect of support-surface perturbations (Lockhart and Ting 2007
; Peterka 2000
). Delayed feedback of horizontal CoM trajectories [displacement, p(t); velocity, v(t); and acceleration, a(t)] were used to stabilize the inverted pendulum (Fig. 1B). EMG reconstructions (EMGp) were taken as the output of the feedback controller, which was a linear combination of the weighted horizontal CoM kinematic trajectories at a common neural transmission delay (
)
![]() | (1) |
For each muscle in each subject, the feedback gains (kp, kv, ka) and delay (
) that best matched the EMG reconstruction to the measured EMG signal were found. We used an optimization (MATLAB, fmincon.m) to find the values of ki and
using the following cost function
![]() | (2) |
Recorded and reconstructed EMG patterns were compared with those predicted by an optimal control model (Lockhart and Ting 2007
). Using a controller design similar to that of the linear quadratic regulator (He et al. 1991
), this delayed quadratic regulator (DQR) model determined gains for CoM kinematic feedback channels, without a priori knowledge of recorded EMG, through the use of a quadratic cost function and time-delayed feedback. Feedback gains on delayed CoM kinematics (ki) were optimized using the following cost function
![]() | (3) |
= 20, requiring the minimum possible level of muscle activation to achieve the postural task. The final term penalized final pendulum configurations that were not upright with weight
. Because the optimization process consistently selected the minimum allowable feedback delay, this delay was set to 100 ms for all subjects to allow the calculation of an intersubject average of the optimal postural control solution and to facilitate qualitative comparisons with recorded and reconstructed EMG patterns. |
|
RESULTS |
|---|
|
|
|
r2 = 0.00 (P = 0.80),
VAF = 0.01 (P = 0.31); MG:
r2 = 0.05 (P = 0.11),
VAF = 0.02 (P = 0.10)].
Acceleration feedback was required to reconstruct EMG activity using physiological delays. When acceleration feedback was removed, delays shorter than the 55-ms latency of the stretch response during postural perturbations (Diener et al. 1984
) were required (intersubject range = 10–60 ms; Fig. 3D). Without acceleration feedback, the early EMG activity in the initial burst and plateau regions, including the initial slope of the response, were underpredicted (data not shown). Further, the goodness of fit between reconstructed and recorded EMGs was reduced in TA [
r2 = –0.14 (P = 7 x 10–4);
VAF = –0.07 (P = 0.006)], but not MG [
r2 = –0.05 (P = 0.42);
VAF = –0.03 (P = 0.13)]. In both cases, however, the reconstructed EMGs without acceleration feedback were often insufficient to maintain the pendulum in an upright configuration (data not shown).
|
|
DISCUSSION |
|---|
|
The nervous system may take advantage of the naturally occurring physical relationships between acceleration, velocity, and displacement to provide feedback control of the CoM without need for feedforward control mechanisms. Previous studies have observed a positive, phase-leading correlation between muscle activity during quiet stance and CoM motion, suggesting the use of predictive, feedforward control (Fitzpatrick et al. 1992
, 1996
; Gatev et al. 1999
). The phase-lead characteristics of acceleration feedback may serve to explain this observation in the context of feedback control. In our model, the contribution of acceleration feedback is fully reflected in the muscular response before significant displacement-related information becomes available. Moreover, the acceleration component of the reconstructed muscular response leads CoM displacement, but occurs after the CoM acceleration induced by the perturbation. The phase lead of acceleration feedback with respect to CoM displacement in our simulations was about 135 ms, consistent with the 100- to 250-ms phase lead observed experimentally for high-frequency postural sway (Fitzpatrick et al. 1992
). The early burst of muscle activity during postural responses to perturbation, here shown to arise from acceleration feedback, was previously attributed to a feedforward component (Diener et al. 1988
). Consistent with our model, however, the middle portion of the response varies with changes in perturbation velocity, whereas the late response is affected by changes in perturbation displacement (Diener et al. 1988
).
Several other studies provide support for acceleration feedback in postural control. Postural responses have been shown to scale with perturbation acceleration in the neck muscles of seated subjects (Siegmund 2004
; Siegmund et al. 2002
) and in perturbations to arm movements (Soechting and Lacquaniti 1988
). In standing posture, muscle onset latency and total ankle moment are also affected by perturbation acceleration (Brown et al. 2001
; Siegmund et al. 2002
; Szturm and Fallang 1998
). Further, the rate of muscle activity onset during perturbations to treadmill walking has also been related to perturbation acceleration (Dietz et al. 1987
). Several studies during standing postural responses suggest that the termination of the postural response results from feedback on the deceleration impulse (Bothner and Jensen 2001
; Carpenter et al. 2005
; McIlroy and Maki 1994
). Consistent with this finding, in our model, termination of the postural response can also be attributed to the delayed effects of the deceleration impulse (Fig. 3A).
Our study supports the idea that a small set of variables related to task-level goals are used to coordinate multiple muscles throughout the body during postural control and other movements. Activity in muscles crossing the hip, knee, and ankle joints all exhibited temporal patterns that were explained by combinations of the CoM motion as modeled by an inverted pendulum. Although the hip and knee joints did not undergo appreciable joint angle changes (Fig. 1A), proximal muscle activity may be necessary to minimize joint motions from interaction torques generated by ankle muscle activity (van Antwerp et al. 2007
; Zajac and Gordon 1989
). Therefore whenever the ankle muscles are activated, the proximal muscles must also be activated to maintain the postural configuration. We propose that a muscle synergy defining consistent spatial patterns of multiple muscle activity for ankle-strategy responses (Torres-Oviedo and Ting 2007
) may be temporally regulated by feedback signals. The spatiotemporal patterns of muscle activation for postural control could thus be specified by defining a constant set of gains on CoM acceleration, velocity, and displacement for each muscle.
Although we have demonstrated the feasibility of task-level feedback in explaining ankle-strategy responses to support-surface translations, more complex biomechanical models may be necessary to represent the full range of responses—ankle, hip, and mixed strategies—in the postural-control suite (Alexandrov et al. 2001b
; Horak and Nashner 1986
; Runge et al. 1999
). This is especially pertinent for modeling muscular responses to backward translations, as well as to support-surface rotations and upper-body perturbations, where hip-strategy responses produce significant joint motions and muscle activation about the proximal joints (Jo and Massaquoi 2004
; Runge et al. 1999
). Because the hip-strategy response has a distinct muscle synergy pattern that can be decomposed from a mixed response (Torres-Oviedo and Ting 2007
), it is possible that the hip-strategy response is also regulated by a task-level feedback controller that is independent of the ankle-strategy controller.
Comparisons of experimentally recorded EMG with an optimal control solution suggest that the postural responses of our human subjects, although similar to the optimal solution, may not have completely achieved the optimal feedback pattern for responding to support-surface translations during the course of our experiment. In contrast, cats subjected to a similar perturbation protocol exhibited EMG patterns that matched the optimal solution as predicted by the DQR model (Lockhart and Ting 2007
). The cats underwent a rigorous training regimen in which they learned to stand on the perturbation platform over the course of several weeks or months (cf. Macpherson et al. 1987
). Our human subjects, however, were completely naïve to postural perturbation studies and each completed the experimental protocol in less than 1 h. We hypothesize that, during their training regimen, the cats may have slowly adapted their muscular responses toward the optimal control solution for the task. We therefore predict that, with training, human muscle activity during postural responses may more closely match the optimal feedback pattern predicted by our DQR model. Alternately, it may be possible that each human subject used a different set of optimality criteria, which could be modeled either by varying the weights in the cost function (Qu et al. 2007
) or by changing the components of the cost function altogether.
Our feedback model may provide a low-dimensional framework for understanding variability in muscle activation patterns during postural control (Ting 2007
). Extensive intersubject variability in temporal patterns of muscle activity may be accounted for by varying only three feedback gains (Fig. 2, A and B). Rather than performing a point-by-point adjustment of neural activity over time, the CNS may adjust gains to each feedback channel. This differential weighting of feedback channels may explain changes in muscle responses due to habituation and changes in central set (Horak et al. 1989
). For example, when the interval between acceleration and deceleration of translation perturbations is short and predictable, subjects anticipate the deceleration timing (Carpenter et al. 2005
; McIlroy and Maki 1994
). The advance in the timing of response termination might occur due to changes in CoM velocity and displacement feedback gains, which alter the time at which the acceleration feedback triggers the offset of EMG activity.
|
|
GRANTS |
|---|
|
|
|
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: L. H. Ting, The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, 313 Ferst Drive, Atlanta, GA 30332-0535 (E-mail: lting{at}emory.edu)
|
|
REFERENCES |
|---|
|
Alexandrov AV, Frolov AA, Massion J. Biomechanical analysis of movement strategies in human forward trunk bending. II. Experimental study. Biol Cybern 84: 435–443, 2001b.[CrossRef][Web of Science][Medline]
Bortolami SB, DiZio P, Rabin E, Lackner JR. Analysis of human postural responses to recoverable falls. Exp Brain Res 151: 387–404, 2003.[CrossRef][Web of Science][Medline]
Bothner KE, Jensen JL. How do non-muscular torques contribute to the kinetics of postural recovery following a support surface translation? J Biomech 34: 245–250, 2001.[CrossRef][Web of Science][Medline]
Brown LA, Jensen JL, Korff T, Woollacott MH. The translating platform paradigm: perturbation displacement waveform alters the postural response. Gait Posture 14: 256–263, 2001.[CrossRef][Web of Science][Medline]
Carpenter MG, Thorstensson A, Cresswell AG. Deceleration affects anticipatory and reactive components of triggered postural responses. Exp Brain Res 167: 433–445, 2005.[CrossRef][Web of Science][Medline]
Diener HC, Dichgans J, Bootz F, Bacher M. Early stabilization of human posture after a sudden disturbance: influence of rate and amplitude of displacement. Exp Brain Res 56: 126–134, 1984.[Web of Science][Medline]
Diener HC, Horak FB, Nashner LM. Influence of stimulus parameters on human postural responses. J Neurophysiol 59: 1888–1905, 1988.
Dietz V, Quintern J, Sillem M. Stumbling reactions in man: significance of proprioceptive and pre-programmed mechanisms. J Physiol 386: 149–163, 1987.
Fitzpatrick R, Burke D, Gandevia SC. Loop gain of reflexes controlling human standing measured with the use of postural and vestibular disturbances. J Neurophysiol 76: 3994–4008, 1996.
Fitzpatrick RC, Gorman RB, Burke D, Gandevia SC. Postural proprioceptive reflexes in standing human subjects: bandwidth of response and transmission characteristics. J Physiol 458: 69–83, 1992.
Gatev P, Thomas S, Kepple T, Hallett M. Feedforward ankle strategy of balance during quiet stance in adults. J Physiol 514: 915–928, 1999.
Gottlieb GL, Chen CH, Corcos DM. Relations between joint torque, motion, and electromyographic patterns at the human elbow. Exp Brain Res 103: 164–167, 1995.[Web of Science][Medline]
He JP, Levine WS, Loeb GE. Feedback gains for correcting small perturbations to standing posture. IEEE Trans Autom Control 36: 322–332, 1991.[CrossRef]
Horak FB, Diener HC, Nashner LM. Influence of central set on human postural responses. J Neurophysiol 62: 841–853, 1989.
Horak FB, Nashner LM. Central programming of postural movements: adaptation to altered support-surface configurations. J Neurophysiol 55: 1369–1381, 1986.
Jo S, Massaquoi SG. A model of cerebellum stabilized and scheduled hybrid long-loop control of upright balance. Biol Cybern 91: 188–202, 2004.[Web of Science][Medline]
Kiemel T, Oie KS, Jeka JJ. Multisensory fusion and the stochastic structure of postural sway. Biol Cybern 87: 262–277, 2002.[CrossRef][Web of Science][Medline]
Kuo AD. An optimal control model for analyzing human postural balance. IEEE Trans Biomed Eng 42: 87–101, 1995.[CrossRef][Web of Science][Medline]
Lockhart DB, Ting LH. Optimal sensorimotor transformations for balance. Nat Neurosci 10: 1329–1336, 2007.[CrossRef][Web of Science][Medline]
Macpherson JM, Horak FB, Dunbar DC, Dow RS. Stance dependence of automatic postural adjustments in humans. Exp Brain Res 78: 557–566, 1989.[Web of Science][Medline]
Macpherson JM, Lywood DW, Van Eyken A. A system for the analysis of posture and stance in quadrupeds. J Neurosci Methods 20: 73–82, 1987.[CrossRef][Web of Science][Medline]
McIlroy WE, Maki BE. The "deceleration response" to transient perturbation of upright stance. Neurosci Lett 175: 13–16, 1994.[CrossRef][Web of Science][Medline]
Park S, Horak FB, Kuo AD. Postural feedback responses scale with biomechanical constraints in human standing. Exp Brain Res 154: 417–427, 2004.[CrossRef][Web of Science][Medline]
Peterka RJ. Postural control model interpretation of stabilogram diffusion analysis. Biol Cybern 82: 335–343, 2000.[CrossRef][Web of Science][Medline]
Qu X, Nussbaum MA, Madigan ML. A balance control model of quiet upright stance based on an optimal control strategy. J Biomech 40: 3590–3597, 2007.[CrossRef][Web of Science][Medline]
Runge CF, Shupert CL, Horak FB, Zajac FE. Ankle and hip postural strategies defined by joint torques. Gait Posture 10: 161–170, 1999.[CrossRef][Web of Science][Medline]
Runge CF, Zajac FE 3rd, Allum JH, Risher DW, Bryson AE Jr, Honegger F. Estimating net joint torques from kinesiological data using optimal linear system theory. IEEE Trans Biomed Eng 42: 1158–1164, 1995.[CrossRef][Web of Science][Medline]
Siegmund GP. Gradation of neck muscle responses and head/neck kinematics to acceleration and speed change in rear-end collisions. Stapp Car Crash J 48: 419–430, 2004.[Medline]
Siegmund GP, Sanderson DJ, Inglis JT. The effect of perturbation acceleration and advance warning on the neck postural responses of seated subjects. Exp Brain Res 144: 314–321, 2002.[CrossRef][Web of Science][Medline]
Soechting JF, Lacquaniti F. Quantitative evaluation of the electromyographic responses to multidirectional load perturbations of the human arm. J Neurophysiol 59: 1296–1313, 1988.
Szturm T, Fallang B. Effects of varying acceleration of platform translation and toes-up rotations on the pattern and magnitude of balance reactions in humans. J Vestib Res 8: 381–397, 1998.[CrossRef][Web of Science][Medline]
Ting LH. Dimensional reduction in sensorimotor systems: a framework for understanding muscle coordination of posture. In: Computational Neuroscience: Theoretical Insights into Brain Function (Progress in Brain Research), edited by Cisek P, Drew T, Kalaska JF. Amsterdam: Elsevier, 2007, vol. 165, p. 301–325.
Torres-Oviedo G, Ting LH. Muscle synergies characterizing human postural responses. J Neurophysiol 98: 2144–2156, 2007.
van Antwerp KW, Burkholder TJ, Ting LH. Inter-joint coupling effects on muscle contributions to endpoint force and acceleration in a musculoskeletal model of the cat hindlimb. J Biomech 40: 3570–3679, 2007.[CrossRef][Web of Science][Medline]
van der Kooij H, Jacobs R, Koopman B, Grootenboer H. A multisensory integration model of human stance control. Biol Cybern 80: 299–308, 1999.[CrossRef][Web of Science][Medline]
Winter DA. Biomechanics and Motor Control of Human Movement. Hoboken, NJ: Wiley, 2005.
Zajac FE, Gordon ME. Determining muscle's force and action in multi-articular movement. Exerc Sport Sci Rev 17: 187–230, 1989.[Web of Science][Medline]
This article has been cited by other articles:
![]() |
S. Kim, F. B. Horak, P. Carlson-Kuhta, and S. Park Postural Feedback Scaling Deficits in Parkinson's Disease J Neurophysiol, November 1, 2009; 102(5): 2910 - 2920. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. J. Chiel, L. H. Ting, O. Ekeberg, and M. J. Z. Hartmann The Brain in Its Body: Motor Control and Sensing in a Biomechanical Context J. Neurosci., October 14, 2009; 29(41): 12807 - 12814. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. D. Goodworth and R. J. Peterka Contribution of Sensorimotor Integration to Spinal Stabilization in Humans J Neurophysiol, July 1, 2009; 102(1): 496 - 512. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. D. J. Welch and L. H. Ting A Feedback Model Explains the Differential Scaling of Human Postural Responses to Perturbation Acceleration and Velocity J Neurophysiol, June 1, 2009; 101(6): 3294 - 3309. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. D. Loram, M. Lakie, and P. J. Gawthrop Visual control of stable and unstable loads: what is the feedback delay and extent of linear time-invariant control? J. Physiol., March 15, 2009; 587(6): 1343 - 1365. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |