The relation between neuromechanical parameters and Ashworth score in stroke patients
- Erwin de Vlugt†1Email author,
- Jurriaan H de Groot†2, 3,
- Kim E Schenkeveld2,
- JHans Arendzen2,
- Frans CT van der Helm1 and
- Carel GM Meskers†2, 3
© de Vlugt et al; licensee BioMed Central Ltd. 2010
Received: 16 November 2009
Accepted: 27 July 2010
Published: 27 July 2010
Quantifying increased joint resistance into its contributing factors i.e. stiffness and viscosity ("hypertonia") and stretch reflexes ("hyperreflexia") is important in stroke rehabilitation. Existing clinical tests, such as the Ashworth Score, do not permit discrimination between underlying tissue and reflexive (neural) properties. We propose an instrumented identification paradigm for early and tailor made interventions.
Ramp-and-Hold ankle dorsiflexion rotations of various durations were imposed using a manipulator. A one second rotation over the Range of Motion similar to the Ashworth condition was included. Tissue stiffness and viscosity and reflexive torque were estimated using a nonlinear model and compared to the Ashworth Score of nineteen stroke patients and seven controls.
Ankle viscosity moderately increased, stiffness was indifferent and reflexive torque decreased with movement duration. Compared to controls, patients with an Ashworth Score of 1 and 2+ were significantly stiffer and had higher viscosity and patients with an Ashworth Score of 2+ showed higher reflexive torque. For the one second movement, stiffness correlated to Ashworth Score (r2 = 0.51, F = 32.7, p < 0.001) with minor uncorrelated reflexive torque. Reflexive torque correlated to Ashworth Score at shorter movement durations (r2 = 0.25, F = 11, p = 0.002).
Stroke patients were distinguished from controls by tissue stiffness and viscosity and to a lesser extent by reflexive torque from the soleus muscle. These parameters were also sensitive to discriminate patients, clinically graded by the Ashworth Score. Movement duration affected viscosity and reflexive torque which are clinically relevant parameters. Full evaluation of pathological joint resistance therefore requires instrumented tests at various movement conditions.
Increased mechanical resistance to an imposed movement is common after central nervous system damage, such as stroke and may interfere with function. Its assessment and treatment are therefore major goals in rehabilitation. Main contributors to increased joint resistance are increased viscosity and stiffness of muscle and connective tissue (clinically labeled "hypertonia") and hyperactivity of the stretch reflex (clinically labeled "spasticity") . The Ashworth Score (AS) is a widely used clinical measure of joint resistance . The AS subjectively grades the manual sensation of mechanical resistance experienced by the examiner during a one second (1 s) joint rotation over the full range of motion . The impossibility to discriminate between the underlying mechanisms and the limited reproducibility and resolution have been the motivating challenge to develop an alternative method describing joint resistance in quantitative neuromechanical measures from the torque response . Discerning muscular and connective tissue properties from the neural reflexes would facilitate the diagnosis of the physiological substrate of increased joint resistance and the subsequent indication for treatment.
Quantitative studies focused on the characteristics of the torque response signals, either versus time or joint angle [2, 5–7]. Peak torque, rate of change and offset of the torque were found to correlate with AS but did not allow for discrimination between individual components of joint resistance. Alternatively, computational models allowed for simultaneous estimation of viscosity, stiffness and reflex torque [8–11]. Critical in such model-based system identification is the structure of the model comprising the relevant neuromechanical components. As in almost any biological system, joint mechanical behavior is highly nonlinear for substantial changes of states, i.e. joint position and velocity, as is the case during e.g. an Ashworth test [12–14]. This implies that a specific linear model structure that is valid for one combination of states will be invalid for almost any other combination. As a consequence, results obtained from small amplitude models [8, 14] may not be generalized to large amplitude conditions. For large amplitude joint rotations, important nonlinear properties such as e.g. the joint angle-dependent stiffness may not be neglected . It is therefore not surprising that different and sometimes conflicting results were reported from different models and types of joint movements [2, 8, 9]. For a valid description of joint neuromechanical behavior during large angular excursions, nonlinear modeling is thus required.
To what extent does duration of an imposed movement affect neuromechanical parameters, i.e. stiffness, viscosity and reflexive torque, in chronic stroke patients and healthy subjects?
Do neuromechanical parameters discriminate between stroke patients and healthy subjects?
Do neuromechanical parameters correlate to disorder severity as graded by the AS?
The clinical relevance of the instrumented identification is to directly attain patients to the appropriate treatment and to be able to quantify the effects of treatment.
Subjects & patients
Seven healthy subjects (mean age 55.4, SD 10.3 years) were recruited as a control group. The medical ethics committee of Leiden University Medical Center approved the study. All participants gave their written informed consent prior to the experimental procedure.
Muscle activation of the tibialis anterior (TA), gastrocnemius lateralis (GL), soleus (SL) and gastrocnemius medialis (GM) was measured by electromyography (EMG) using a Delsys Bagnoli 4 system. Inter electrode distance was 10 mm. EMG signals were sampled at 2500 Hz, on-line band pass filtered (20-450 Hz) and off-line rectified and integrated by low pass filtering (3th-order Butterworth) at 20 Hz (IEMG). Reaction torque and ankle angle were sampled at 250 Hz. Angular velocity and acceleration were derived by single and double differentiation of the recorded angle signal respectively. To avoid amplification of noise due to differentiation, angle and force signals were low pass filtered at 20 Hz (3th-order Butterworth).
1. Clinical test
Measurements were performed on the affected ankle of each patient and at the right ankle in case of controls. The Ashworth Score (AS) of the affected ankle  was assessed by an experienced physician [HA]. In order to avoid obtaining a biased and a study-specific Ashworth test, the physician was instructed to perform the Ashworth test as he would perform as usual in the clinic. Total time to perform the Ashworth test including positioning and instructing of the patient was about 5 minutes. The instrumented rotation measurements were performed by an experimenter [KS] who was blind to the clinical outcome. Judgment on the validity of the model was solely based on the recorded signals (internal validity). For the control group, only the instrumented measurements were performed. All measurements were completed within a single session of approximately one hour.
2. Instrumented joint rotation
The ankle angle was defined as the position of the foot with respect to the lower leg; the perpendicular position was defined as zero degrees or central position. Maximum dorsiflexion angle was assessed by a monotonically in- and decreasing dorsiflexion torque (100 s up, 100 s down) imposed by the manipulator from zero to a maximum value of 15 Nm resulting in slow rotations of approximately 0.5 deg/s. The angle before onset of the dorsiflexion torque was taken as the maximal plantar flexion angle. The angular excursion in plantar flexion direction was limited to -30 degrees, which was the maximal angle of the manipulator. RoM was defined as the difference between the maximum dorsiflexion and plantar flexion angle and used as boundary for the subsequent RaH rotations. At 15 Nm the foot was approximately at a perpendicular angle with respect to the horizontal for all subjects. Consequently, the variability in torque introduced by gravity around the maximal dorsiflexion angles could be considered negligible and thus there was no need to compensate for gravity during these tests.
RaH rotations were performed by the manipulator through the full RoM at four different durations of 0.25, 0.5, 1 and 2 s. As RoM differed between subjects while durations were fixed, rotation velocities were different between subjects. Prior to each RaH rotation, the ankle was moved from central position to the maximal plantar flexion angle in 2 s time. Subsequently, at a random time instant but within 3 to 4 s, the RaH rotation was started. In all cases, the RaH rotation ended at the maximal dorsiflexion angle. The hold phase lasted for 4 s after which the ankle was moved back again to the central position. Time to cover a complete movement profile did not exceed 15 s. Rest periods of 30 s were maintained between each movement profile which is sufficient for full recovery of passive stiffness . All movement profiles were performed twice to test for repeatability of the estimation procedure. Subjects were asked to remain maximally relaxed during the entire experiment and not actively resist any motions. Level of relaxation was checked off-line from EMG activity of all muscles prior to the RaH rotation. When IEMG was larger than three times standard deviation for longer than 1 s the observation was discarded from the analysis.
Neuromechanical model, parameter estimation and internal validity
(mean ± 1 s.d.)
mass (ankle + footplate)
1.86 ± 0.42
1.28 ± 1.08
26.4 ± 15.4
muscle length shift
-0.0081 ± 0.0023
muscle force shift
-21.2 ± 9.6
e 1 , e 2 ,
e 3 , e 4
EMG weighting factors
3.5 ± 1.05, 2.0 ± 0.96,
3.1 ± 0.77, 2.6 ± 1.1 (× 105)
activation cutoff frequency
1.28 ± 0.34
where N is the number of time samples used for estimation of the parameters, J the Jacobian matrix, and e the 1 × N error vector. The Jacobian is a N × n p matrix, with n p = 10 the number of estimated parameters, containing first derivatives of the (final) error to each parameter.
Two different type of indicators were derived from the covariance matrix. The first is the interdependence of the parameters for which the auto-covariance (diagonal terms of P) of each parameter was compared to the cross-covariance (off-diagonal terms of P) between the one parameter and all the others. If the auto-covariance was higher than all cross-covariances, the corresponding parameter was estimated/assumed independently and its estimated value was assumed to be reliable. The second measure is the sensitivity of the parameters for which the auto-covariance value on itself is representative. High sensitivity means that the parameter has an observable contribution in the system's response (i.e. the ankle torque in this study) and therefore can be estimated with certain accuracy. The square root of the auto-covariance, such as obtained from P in the above expression, is the standard error of the mean (SEM) of the parameter estimation . For high sensitivity, the SEM needs to be low compared to the corresponding parameter value.
For visual inspection, we have normalized the covariance matrix by dividing each i,j-th element by (i, j from 1 to n p ) such that all diagonal terms equal to one. SEM values were normalized to their corresponding parameter values and subsequently averaged over all trials and subjects.
with T meas (t) the measured ankle reaction torque and T mod (t) the estimated ankle torque from the model (Eq. A1, Appendix) over the time frame used for parameterization.
and similarly for the reflex torque of the tibialis anterior, with n indicating the time sample of the identification time frame [1 ... N]. The r.m.s. value is a common way to denote the energy of a signal.
The model parameters were defined on the (metric linear) muscle level while for interpretation and analysis of the results, viscosity and stiffness were expressed in the (angular) joint domain according to Eqs A10 and A11 (Appendix). Viscosity and stiffness increase exponentially with joint angle (muscle length). Because of the exponential relationship, both viscosity and stiffness could only be compared at the same joint angle, θcomp, for all subjects (controls and patients). θcomp was determined by the smallest maximal dorsiflexion angle amongst all subjects. Any differences in viscosity and/or stiffness between subjects and patients was largest at θcomp. Statistical testing of viscosity and stiffness at smaller joint angles was therefore considered less meaningful, hence not performed.
For statistical analysis, a disease gradation was defined, ranging from healthy subjects to patients graded by AS. Thus, within the tested population, four groups were discerned, i.e. controls (C), a clinically unaffected patient group: AS0; a mildly affected patient group: AS1; and a severely affected patient group, i.e. the patients exhibiting an AS of 2 and higher: AS2+.
To test the differences in RoM between patients graded by AS and controls, a one way ANOVA was used with a Bonferroni post hoc test. Movement duration and velocity were separately related with the RoM. As RoM differed between subjects, duration and velocity were not interchangeable. Movement duration was standardized and thus the factor duration (not velocity) was applied in the analysis. To test the effects of movement duration and disease gradation, a Linear Mixed Model was used with disease gradation as fixed and movement duration as repeated factor. In case of significant effects of either factor, a Bonferroni post hoc test was used to specify the differences between the groups. Correlation between relevant neuromechanical parameters and AS was assessed using linear regression. All statistical testing was performed using SPSS 16.0, SPSS Inc. at an alpha of 0.05.
Both Controls and Patients could perform the tests. No problems were observed with cognitive or language deficits interfering with the comprehension of instructions required to participate in the study. A total of 10 trials from three healthy subjects were removed from the analysis because of sudden and large IEMG bursts of all muscles before the onset of the RaH movements, indicating insufficient relaxation.
Range of Motion (RoM)
All patients and controls reached to the maximal plantarflexion angle of -30 degrees, which was the limit of the manipulator. Consequently, all the observed loss in RoM was accounted for by the reduced dorsiflexion.
To check for stretch induced muscle activity that might have affected the RoM measurement, the mean IEMG at zero torque (before dorsiflexion torque was imposed) was compared to the mean IEMG at the maximal dorsiflexion torque. Mean IEMG was taken over a 1 s interval and was larger at 15 Nm than at zero torque for almost all subjects. However, the increments were small (0.5-1%) relative to the magnitude of the IEMG responses observed during the RaH movements (see further). Therefore, the small IEMG increment during the RoM measurements were considered to have a negligible effect on the reported RoM values.
Torque response to ramp-and-hold movement
Amongst all muscles, the soleus showed the highest activity in response to the imposed movements. Muscle activity emerged in brief bursts that increased in magnitude with shorter movement duration.
The composition of the net muscle activity from the individual IEMG signals is presented in Figure 5 (same subjects and conditions as in Figure 4; recordings in grey and model estimates in black). TA activity was absent. For the stroke patient, soleus activity showed distinct bursts and dominated the net estimated activity of the triceps surae. The estimated contribution of the three calf muscles to the total estimated reflexive torque (Figure 5 M), as obtain from the optimized weighting factors (e 2, e 3 and e 4 ) was 3%, 91% and 6% for the GL, SL and GM respectively. Comparable distribution of muscle torque amongst the triceps surae was found for all other subjects and patients.
Model validity and parameter accuracy
On the average, the SEM was less than 10% except for the IEMG weighting factors (Figure 6, bottom). The weighting factors of both gastrocnemii (e 2 and e 4 ) were least sensitive.
Estimated mass (1.86 ± 0.42 kg), muscle length shift (-0.0081 ± 0.0023 m), muscle force shift (-21.2 ± 9.6 N) and activation cut-off frequency (1.28 ± 0.34 Hz) did not change significantly with movement duration and also were not different between the patients and the control group. Viscosity and stiffness coefficients and reflex torque markedly differed as described in the following sections. Table 2 summarizes the initial and averaged (optimal) estimated values of all model parameters.
Influence of movement duration
Difference between patients and controls
Ankle viscosity (F = 20.2, p < 0.0001), stiffness (F = 19.5, p < 0.0001) and reflexive torque of the triceps surae (F = 5.8, p = 0.003) differed with disease grade. Post hoc testing revealed that for ankle viscosity and stiffness, control subjects could be discerned from stroke patients with an AS of 1 and higher; for reflexive torque, controls differed significantly from patients with an AS2+.
Interaction of disease grade and test condition
Reflexive torque of the triceps surae decreased with duration and this effect was stronger for patients with higher AS (Figure 9, top, interaction term F = 2.91, p = 0.013). At the 1 s movement duration, stiffness significantly related to AS (r2 = 0.51, F = 32.7, p < 0.001) while reflex torque did not (r2 = 0.09, F = 3.22, p = 0.08). At shorter durations, reflex torque significantly related to disease grade (r2 = 0.25, F = 11, p = 0.002). Reflex torque from tibialis anterior did not relate to movement duration nor to AS.
The overall aim of this study was to estimate neuromechanical parameters at the ankle joint in stroke patients during ramp-and-hold (RaH) rotations with different duration using a nonlinear dynamic ankle model. The experiments included the Ashworth test condition: a typical 1 s rotation over the full range of motion, which is clinically used to judge joint resistance in spasticity.
Influence of movement duration on neuromuscular properties
Stretch reflex torque from the triceps surae showed a marked threshold in the movement duration in between 0.5 - 1.0 s, above which there was no substantial reflex response observed (Figure 9, top). The increase of reflexive torque from the triceps surae with movement duration beyond the threshold was expected for it is consistent with the well known velocity dependence of the stretch reflex .
The only other parameter that was influenced by movement duration, albeit slightly, was joint viscosity (Figure 8, top). The slower the joint was rotated the larger its viscosity (velocity to force relation). The increased viscosity was significant only for the longest (2 s) duration indicating to a nonlinear relationship.
Difference between controls and patients
Stiffness, viscosity and reflexive torque from the triceps surae significantly differed between controls and the stroke patients with an AS of one and higher. Increased stiffness was not significantly higher for patients with AS0 compared to controls, indicating small differences with a statistical problem of power.
Although subjects were instructed to relax and not react to the RaH movements, stroke patients may have exhibited an increased ankle torque due to a possible higher background activity of the muscles at rest, as was reported by . Also, an increase in stiffness from within the interior of the muscle cell was found in spastic muscle tissue and which is believed to originate from altered strain properties of intracellular proteins like titin [19, 20]. We assumed that the increased stiffness in the stroke patients as found in this study was mainly from intracellular tissues since the observed stiffness behavior was well described by an exponential force-length relationship (Eq. A9) that is typical for passive tissues [13, 21–23]. Increased stiffness at joint positions beyond the 'relaxed' position is believed to underlie contractures (muscle shortening) as observed in spastic patients [19, 20].
Disease severity is expressed by tissue stiffness in stroke
Intrinsic ankle stiffness was responsible for the increased AS in stroke patients. This means that joint resistance, as was indicated by the AS, is accounted for by the physical property 'stiffness', which is most likely originating from passive tissues. For the extent that AS provides a measure of disease severity, at least for the changes within the mechanical condition of the joint secondary to the neural disorder, we now may state that stiffness of the passive tissues increases with disease severity in stroke.
Ashworth Scale does not comprises the stretch reflex response
Mechanical joint resistance is never determined by passive stiffness only, since reflexive torque was present during all applied RaH movements. However, for the two longest movement durations lasting 1 s, i.e. the Ashworth test duration, and 2 s the reflexive contributions were small. At shorter movement durations of 0.5 s and 0.25 s, the reflex torque from the triceps surae increased with AS.
Ashworth test versus instrumented ramp-and-hold movements
It is important to realize that the manual performance of the Ashworth test may differ from the instrumented ramp-and-hold movements as applied in the present study. The instrumented conditions were of a constant velocity (ramp phase) whereas imposed manual manipulations may result in a bell-shaped velocity profile . Therefore, the instrumented tests in this study are to be considered as separate tests next to the Ashworth test. Direct comparison to the Ashworth test must be taken with some care, but only for those properties that appeared to be dependent on movement velocity being joint viscosity and the stretch reflex torque, as was discussed above.
For the sake of direct comparison to the AS, movement duration was chosen to be the independent controlled variable, but resulted in different velocities between patients and controls. Thus, a structural bias with higher Controls velocities (because of increase RoM) was included in the inter-subject analysis of viscosity and triceps surae reflex torque. If velocity was controlled for, viscosity would likely exhibit less differences between controls and patients and less interaction with disease grade (AS). For the triceps surae reflex torque, the opposite would occur: differences between controls and patients, and in between AS groups, would be larger if velocity was controlled for. Although viscous torques have a marginal contribution to the overall joint torque in comparison to the stiffness and reflex torques, the bias problem requires the inter-subjective significance of (only) the tissue viscosity to be taken with care.
However, discrepancy between the description of the Ashworth test and the actual manual performance underlines the necessity of applying controlled test conditions to obtain reliable and valid outcome parameters.
Validity of the method
The full model consisted of 10 parameters that were estimated reliably as indicated by the low interdependence values (Figure 6, top). For all but the IEMG weighting factors, the sensitivity was high (low SEM values). The combination of low interdependence and high sensitivity indicates that these parameters were estimated reliably and accurately.
Both viscosity and stiffness coefficients decreased 17% and 12% on average respectively with repetition. This decrease in passive joint visco-elasticity with ongoing loading was previously reported in both the normal and spastic case, e.g. [23, 25, 26]. Also, the length shift parameter was 9% larger with repetition which means that the ankle joint angle beyond which the visco-elastic torque started to increase shifted to dorsiflexion, which is probably related to the decrease in visco-elasticity.
The force shift coefficient (F 0 ) had influence on all parameters (last column of the covariance matrix). Based on the small interdependence amongst most other parameters, it follows that the estimation of F 0 was influenced by the other parameters to some extent. The prime role of F 0 was to shift the exponential force-length characteristic to have more flexibility in describing the ankle stiffness but perhaps it was also used to account for small model remnants.
The IEMG weighting factors, in particular these for the gastrocnemii muscles, were least sensitive while their interdependence was exceptionally small. This means that the contributions of the gastrocnemii could be estimated independently but their estimated contributions to reflexive torque were far less compared to the soleus muscle. The intertrial difference for the soleus was smallest (12 ± 20%) which confirmed its dominant contribution to triceps surae reflex torque compared to the gastrocnemii muscles.
Because the gain (participation) of each EMG channel was also estimated, the method was free to select which muscles contributed and to what extent. Any cross-talk between agonists (soleus and both gastrocnemii) was therefore of no problem. Cross-talk between antagonistic muscles may have disturbed the selection between muscles. However, it has been shown that there is 5% cross-talk from the tibialis to the soleus at most and under supra maximal stimulation . It was not likely that supra maximal activation occurred during our experiment so any effect of cross-talk was most likely very small.
In our model, the Achilles and tibialis tendons were taken as infinitely stiff. Over all subjects and patients plantarflexion torque never exceeded 30 Nm. In normal subjects maximal voluntary contraction (MVC) produces about 150-225 Nm (female-male) of plantarflexion torque at 10 degrees dorsiflexion . Thus, plantarflexion torque was in the range of 13-20% MVC of normal, resulting in a maximal tendon elongation of 0.4-0.6 cm respectively . The total muscle-tendon length change followed directly from the ankle angle and moment arm and varied in the range of 3.5-4.7 cm, which means that 17% of the muscle-tendon length change would be from the Achilles tendon at most. As the consequence of omission of the Achilles tendon in our model, joint stiffness and viscosity values may be slightly underestimated since we assumed one element instead of two elements in series. An infinitely stiff Achilles tendon has also been assumed in previous studies that estimated neuromechanical properties of the elbow and ankle joint [8, 9, 14, 30].
Overall model validity is illustrated by high "goodness of fit" (VAF) values that were above 90% for all movement conditions tested. Together with the low interaction between the parameters and high sensitivity of the model parameters we therefore may conclude that the underlying neuromechanical behavior of the ankle joint was well quantified by the model for all conditions tested.
Comparison to the literature
Increased stiffness was also observed in a comparable study  in the paretic limb of stroke patients, but which did not increase with AS. In that particular study, continuous (> 30 s) small amplitude joint rotations (1.5 degrees) were applied at high speeds. Since the joint system (as any biological system) is highly nonlinear [13, 31] and varies as a function of time, long lasting small amplitude behavior cannot be generalized or extrapolated to brief (< 2 s) large amplitude behavior (> 15 degrees) as used in an Ashworth test and applied in the current study.
Ankle joint viscosity had a mean value of 0.69 Nms/rad and 1.14 Nms/rad for the control and patient group respectively, which are in the same ranges as found previously by . Mean ankle joint stiffness was 14 Nm/rad for the control group and 31 Nm/rad for the stroke patients, which are both a factor 3 to 4 lower than found by  and for the controller group a factor 3 lower than found by . The discrepancy can be explained from the usage of much smaller displacements (several degrees) in [13, 14] as passive joint stiffness strongly increases with decreasing amplitude of displacement [13, 31].
The mean estimated mass was 1.86 kg and modeled as a point mass at a fixed distance of 0.15 m from the rotation center of the ankle joint, i.e. the inertia was 0.042 Nms2/rad. The inertia of the footplate was 0.032 Nms2/rad such that the mean foot inertia was 0.010 Nms2/rad, which is only slightly higher than the range of 0.007-0.009 Nms2/rad as previously reported [32, 33].
In previous studies, reflex contributions to ankle torque were estimated by quantification of parameters from a feedback model representing the functioning of the muscle spindles and Golgi tendon organs [8, 14, 34]. However, the inputs to proprioceptive sensors such as length (velocity) and force of the muscle fibers cannot be measured during movement experiments accurately to date. This study directly estimates reflex torque from measured muscle activity (IEMG) and therewith no assumptions about the functioning of the proprioceptive sensors were required, allowing for direct estimations of the net reflex torque.
Ankle stretch reflexes in stroke patients were previously found above angular velocities of 80 deg/s , which is comparable to the RaH rotations of 0.5 s duration (≈ 40 deg. ROM in 0.5 s) in this study. Rotation velocity during the Ashworth test, performed as a 1 s full RoM movement , is therefore assumed to be sub-threshold not evoking stretch reflexes.
The muscle length shift parameter x 0 was used for shifting the exponential stiffness function with muscle length and can be interpreted as the muscle length at which the passive elastic force starts to increase substantially. The shift was -8.1·10-3 m on average. In , the ankle angle at which passive plantarflexion torque started to increase rapidly was approximately 0.4 rad plantarflexion (-23 deg). For our model, the angle at which the passive stiffness torque started to increase was for that muscle length x where the exponential power term x - x 0 (Eq. A8) was zero, that is for x = -x 0 = 8.1·10-3 m. From Eqs. A4 and A5 it follows that this value for x corresponded to an angle of -0.43 rad, which is close to the referred value above. The shift parameter can be interpreted as a physiological meaningful parameter describing the passive elasticity property of the triceps surae and was not different for the stroke patients compared to the controls. Apparently, the increase in passive tissue stiffness in the stroke patients was fully described by the (increased) curvature parameter k of the stiffness force-length relationship (Eq. A8).
Cut-off frequencies of second order models describing muscle activation dynamics have been reported in several previous studies. Most of these studies found values ranging from 1 to 3.3 Hz. In  maximal values around 7 Hz were found for the ankle triceps, which seems too high to our opinion. The mean value of 1.28 Hz as reported in our study is within the range of 1.0 - 1.4 Hz as found by  and somewhat lower than the cut-off frequencies found for the trunk (2.0 - 3.3 Hz)  and for elbow muscles (1.9 - 2.8 Hz) , likely because the soleus muscle is composed largely of slow twitch muscle fibers.
Muscle activity in response to the imposed (fastest) movements was observed as distinctive bursts (Figure 5: I) and were more pronounced for the stroke patients (not shown). Similar bursts of activity during comparable joint movements were reported by others [23, 38, 39]. Likely, the motoneurons in stroke patients tend to fire in a more synchronized way in response to afferent input from the stretch receptors that may be the result of decreased motoneuron thresholds  or increased sensitivity of afferent inputs .
In , a similar nonlinear relationship was found for the ankle in SCI patients with largest increase in viscosity below 20 deg./s, which is in the same range as the velocities during the 2 s movements (~ 40 deg.) in the current study. Viscous behavior of connective tissues (intra and extra muscular)  and a possible small amount of actin-myosin cross-bridges in the resting muscle  may have contributed to the velocity dependent behavior. The relationship between movement velocity and joint viscosity remains to be solved and may be important for understanding energy dissipation in functional tasks, e.g. during walking.
The current findings that joint viscosity and reflex torque depended on the duration, and thus the velocity of movement, implicate that for unambiguous assessment of joint resistance the Ashworth test should be performed in a strictly standardized way, actually according to a prescribed velocity instead of a 1 s movement. However, stretch velocity is difficult to standardize in manual testing. Instrumented evaluation comprising extended experimental conditions in combination with nonlinear computational modeling may prove to be a powerful tool to evaluate joint function.
Instrumented tests, like the one applied here, facilitate assessment of quantitative and objective ranges of neuromechanical properties correlating to disorder severity and may guide the clinician in optimal treatment planning e.g. choosing a stiffness reducing strategy instead of reducing reflex activity.
Functional evaluation, e.g. during walking, is compulsory for treatment guidance which can not be extrapolated from passive movements as studied here. We prepare for a larger study to compare neuromuscular properties as measured during static (sitting) and dynamic (walking) conditions to determine to what extent static measures can be used to predict functional improvement during dynamic conditions.
Contribution to joint stiffness and viscosity from any muscle background activity could not be explicitly separated by the current model. That is, all angular velocity and angle related intrinsic torques where lumped together into a viscous and a stiffness torque component respectively. A further division between passive visco-elastic torque and torque emerging from (constant) muscle activation is planned for future studies. To determine its clinical value, we plan to apply the current method to a larger cohort of patients to study the effect of different interventions on neuromuscular properties of the ankle joint.
This study demonstrated a new measurement technique for quantification of neuromechanical parameters of the individual ankle joint from a single dorsiflexion movement. Tissue and reflex torque were most sensitive parameters to discriminate stroke patients from healthy control subjects and also "grade" patients. Stroke patients exhibited increased ankle stiffness and viscosity with AS. For movement durations shorter than 1 s stroke patients also showed increased reflex torque with AS. Joint resistance observed during the 1 s movement over its RoM originated mainly from increased tissue stiffness. Correlations of relevant parameters to AS were assessed on group level and the relatively high standard deviations illustrate the difficulty experienced in discrimination between AS grades in the clinical practice.
The developed model fully covered the observed neuromechanical behavior of the ankle joint. It provides a basis for further dividing the visco-elasticity into contributions from connective and (active) muscle tissue, and the reflex torque into contributions from muscle spindles and Golgi tendon organs. The present study was primarily aimed at development of the method. Inclusion of larger and more divergent patient groups will demonstrate whether clinical phenotypes can be identified in (combinations of) abnormal system properties, such as enhanced stiffness and reflex torque. This may then be the foundation for therapy guidance, e.g. splinting, casting or surgery versus botulinum toxin. Establishing the sensitivity to interventions is a first step towards therapy evaluation.
We conclude that the combination of instrumented evaluation including multiple experimental conditions and nonlinear computational modeling is a powerful tool to quantitatively assess joint resistance. Objective and high resolution identification of neuromuscular parameters will be of use in daily clinical practice.
Appendix 1: Neuromuscular model
where t is the independent time variable [s], T mod the modeled ankle reaction torque [Nm], the ankle angular acceleration [rad/s2], I the inertia of ankle plus footplate [kg.m2], T tri the torque generated by the plantar flexion muscles (GL, SL, GM), or triceps surae [Nm], T tib the torque generated by the dorsiflexion muscle (TA) [Nm], and T grav the torque due to gravity [Nm].
Although the TA was not substantially stretched during the ramp phases in the current experiment there was considerable reflex activity during some RaH movements. For these reasons, viscosity and stiffness were modeled for the plantar flexor muscles only and reflexive force was included for both plantar and dorsiflexion muscles.
where x is the (change) of muscle length (linear displacement) [m], the rate of change of muscle length [m/s], F visc the velocity related muscle force from tissue viscosity [Ns/m], F stiff the length related muscle force from tissue stiffness [N/m], F reflex,tri and F reflex, tib the reflexive muscle forces from the triceps surae and TA respectively [N], and r achil (θ) and r tib (θ) the angle dependent moment arm [m] of the Achilles tendon and tibialis anterior tendon respectively.
where θ fgnd represents the angle of the foot with respect to the horizontal (ground) at zero degrees ankle angle [rad]. Here, θ fgnd equals to the angle the point mass had (around the ankle rotation axis) with the horizontal in central position and which was fixed at 10 degrees, and g is the gravitational acceleration (g = 9.8 m/s2).
Force due to stiffness (Eq. A9) exponentially increases with ankle angle corresponding to the length tension properties of ligamentous, tendinous and muscular elastic tissues [12, 23, 47–49]. Increased tissue stiffness, as often seen in spasticity , can be described by Eq. A9 as a steeper (or shifted) force-length relationship. We assume viscous forces of tissues along the ankle joint to relate to compression (shear forces), which increase with tension. Therefore, both viscous and stiffness force scale with position (Eq. A8, A9). Exponential increase in viscous force with joint angle was also derived from . The exponential curvature is shaped by k [1/m], called the stiffness coefficient, while the amount of viscosity is obtained by multiplication the same curvature with the viscosity coefficient b [Ns/m]. Two shift parameters are included in Eqs A8 and A9 such that the viscous and stiffness forces can be shifted in two dimensions, that is, in length by x 0 and in force by F 0 . The muscle length beyond which the force starts to increase exponentially is determined by the shift parameter x 0 , a known property of passive muscle stiffness. The offset force term, F 0 , served purely as a shaping parameter for the stiffness model.
where xcomp the triceps surae muscle length corresponding to θcomp.
where u tib and u tri the neural activity for the tibialis and triceps surae respectivly, e 1 - e 4 are weighting factors [N/Volts], numerical subscripts (1 - 4) correspond to the IEMG signals of the four muscles as referred to by subscripts TA, GL, SL and GM respectively.
and similarly for the tibialis anterior, where α tri the active state of the triceps surae, ω0 = 2π f 0 the cutoff frequency of the activation filter, and s the Laplace operator denoting the first time derivative. The relative damping β was set to one (critically damped) .
and similarly for the tibialis anterior. For a full description of the structure and the parameters of the force-velocity and force-length relationships of the model (Eq. A15) we refer to . The most important parameter values were: optimum muscle length triceps surae (tibialis) 3.5 (4.6) cm occurring at central ankle angle; maximum shortening velocity 8 (8) times optimum muscle length; maximum eccentric force was 1.5 (1.5) times the isometric force, and the isometric force was normalized to 1 since scaling of force was fully determined by the weighting factors e 1 -e 4 .
This study was performed as part of the Dutch TREND project (Trauma RElated Neuronal Dysfunction), supported by the Dutch Government (grant nr. BSIK03016).
- Katz RT, Rymer WZ: Spastic hypertonia: mechanisms and measurement. Arch Phys Med Rehabil 1989, 70: 144-155.PubMedGoogle Scholar
- Damiano DL, Quinlivan JM, Owen BF, Payne P, Nelson KC, Abel MF: What does the Ashworth scale really measure and are instrumented measures more valid and precise? Dev Med Child Neurol 2002, 44: 112-118. 10.1017/S0012162201001761View ArticlePubMedGoogle Scholar
- Ashworth B: Preliminary Trial of Carisoprodol in Multiple Sclerosis. Practitioner 1964, 192: 540-542.PubMedGoogle Scholar
- Pandyan AD, Johnson GR, Price CI, Curless RH, Barnes MP, Rodgers H: A review of the properties and limitations of the Ashworth and modified Ashworth Scales as measures of spasticity. Clin Rehabil 1999, 13: 373-383. 10.1191/026921599677595404View ArticlePubMedGoogle Scholar
- Pandyan AD, Price CI, Barnes MP, Johnson GR: A biomechanical investigation into the validity of the modified Ashworth Scale as a measure of elbow spasticity. Clin Rehabil 2003, 17: 290-293. 10.1191/0269215503cr610oaView ArticlePubMedGoogle Scholar
- Starsky AJ, Sangani SG, McGuire JR, Logan B, Schmit BD: Reliability of biomechanical spasticity measurements at the elbow of people poststroke. Arch Phys Med Rehabil 2005, 86: 1648-1654. 10.1016/j.apmr.2005.03.015View ArticlePubMedGoogle Scholar
- Chung SG, van Rey E, Bai Z, Rymer WZ, Roth EJ, Zhang LQ: Separate quantification of reflex and nonreflex components of spastic hypertonia in chronic hemiparesis. Arch Phys Med Rehabil 2008, 89: 700-710. 10.1016/j.apmr.2007.09.051View ArticlePubMedGoogle Scholar
- Alibiglou L, Rymer WZ, Harvey RL, Mirbagheri MM: The relation between Ashworth scores and neuromechanical measurements of spasticity following stroke. J Neuroeng Rehabil 2008, 5: 18. 10.1186/1743-0003-5-18PubMed CentralView ArticlePubMedGoogle Scholar
- McCrea PH, Eng JJ, Hodgson AJ: Linear spring-damper model of the hypertonic elbow: reliability and validity. J Neurosci Methods 2003, 128: 121-128. 10.1016/S0165-0270(03)00169-9PubMed CentralView ArticlePubMedGoogle Scholar
- de Vlugt E, Schouten AC, van der Helm FC: Adaptation of reflexive feedback during arm posture to different environments. Biol Cybern 2002, 87: 10-26. 10.1007/s00422-002-0311-8View ArticlePubMedGoogle Scholar
- van der Helm FC, Schouten AC, de Vlugt E, Brouwn GG: Identification of intrinsic and reflexive components of human arm dynamics during postural control. J Neurosci Methods 2002, 119: 1-14. 10.1016/S0165-0270(02)00147-4View ArticlePubMedGoogle Scholar
- Harlaar J, Becher JG, Snijders CJ, Lankhorst GJ: Passive stiffness characteristics of ankle plantar flexors in hemiplegia. Clin Biomech (Bristol, Avon) 2000, 15: 261-270. 10.1016/S0268-0033(99)00069-8View ArticleGoogle Scholar
- Weiss PL, Kearney RE, Hunter IW: Position dependence of ankle joint dynamics--I. Passive mechanics. J Biomech 1986, 19: 727-735. 10.1016/0021-9290(86)90196-XView ArticlePubMedGoogle Scholar
- Mirbagheri MM, Barbeau H, Ladouceur M, Kearney RE: Intrinsic and reflex stiffness in normal and spastic, spinal cord injured subjects. Exp Brain Res 2001, 141: 446-459. 10.1007/s00221-001-0901-zView ArticlePubMedGoogle Scholar
- Gajdosik RLaL, DJ, McFarley DC, Meyer KM, Riggin TJ: Dynamic elastic and static viscoelastic stress-relaxation properties of the calf muscle-tendon unit of men and women. Isokinetics and Exercise Science 2006, 14: 33-44.Google Scholar
- Ljung L: System Identification - Theory for the User. second edition. New Jersey: Prentice Hall; 1999.Google Scholar
- Houk JC, Rymer WZ, Crago PE: Dependence of dynamic response of spindle receptors on muscle length and velocity. J Neurophysiol 1981, 46: 143-166.PubMedGoogle Scholar
- Burne JA, Carleton VL, O'Dwyer NJ: The spasticity paradox: movement disorder or disorder of resting limbs? J Neurol Neurosurg Psychiatry 2005, 76: 47-54. 10.1136/jnnp.2003.034785PubMed CentralView ArticlePubMedGoogle Scholar
- Lieber RL, Steinman S, Barash IA, Chambers H: Structural and functional changes in spastic skeletal muscle. Muscle Nerve 2004, 29: 615-627. 10.1002/mus.20059View ArticlePubMedGoogle Scholar
- O'Dwyer NJ, Ada L, Neilson PD: Spasticity and muscle contracture following stroke. Brain 1996,119(Pt 5):1737-1749. 10.1093/brain/119.5.1737View ArticlePubMedGoogle Scholar
- Gajdosik RL: Influence of a low-level contractile response from the soleus, gastrocnemius and tibialis anterior muscles on viscoelastic stress-relaxation of aged human calf muscle-tendon units. Eur J Appl Physiol 2006, 96: 379-388. 10.1007/s00421-005-0091-7View ArticlePubMedGoogle Scholar
- Gajdosik RL, Vander Linden DW, McNair PJ, Riggin TJ, Albertson JS, Mattick DJ, Wegley JC: Viscoelastic properties of short calf muscle-tendon units of older women: effects of slow and fast passive dorsiflexion stretches in vivo. Eur J Appl Physiol 2005, 95: 131-139. 10.1007/s00421-005-1394-4View ArticlePubMedGoogle Scholar
- Singer BJ, Dunne JW, Singer KP, Allison GT: Velocity dependent passive plantarflexor resistive torque in patients with acquired brain injury. Clin Biomech (Bristol, Avon) 2003, 18: 157-165. 10.1016/S0268-0033(02)00173-0View ArticleGoogle Scholar
- Rabita G, Dupont L, Thevenon A, Lensel-Corbeil G, Perot C, Vanvelcenaher J: Differences in kinematic parameters and plantarflexor reflex responses between manual (Ashworth) and isokinetic mobilisations in spasticity assessment. Clin Neurophysiol 2005, 116: 93-100. 10.1016/j.clinph.2004.07.029View ArticlePubMedGoogle Scholar
- McNair PJ, Hewson DJ, Dombroski E, Stanley SN: Stiffness and passive peak force changes at the ankle joint: the effect of different joint angular velocities. Clin Biomech (Bristol, Avon) 2002, 17: 536-540. 10.1016/S0268-0033(02)00062-1View ArticleGoogle Scholar
- Yeh CY, Chen JJ, Tsai KH: Quantifying the effectiveness of the sustained muscle stretching treatments in stroke patients with ankle hypertonia. J Electromyogr Kinesiol 2007, 17: 453-461. 10.1016/j.jelekin.2006.07.001View ArticlePubMedGoogle Scholar
- De Luca CJ, Merletti R: Surface myoelectric signal cross-talk among muscles of the leg. Electroencephalogr Clin Neurophysiol 1988, 69: 568-575. 10.1016/0013-4694(88)90169-1View ArticlePubMedGoogle Scholar
- Gravel D, Richards CL, Filion M: Angle dependency in strength measurements of the ankle plantar flexors. Eur J Appl Physiol Occup Physiol 1990, 61: 182-187. 10.1007/BF00357596View ArticlePubMedGoogle Scholar
- Muramatsu T, Muraoka T, Takeshita D, Kawakami Y, Hirano Y, Fukunaga T: Mechanical properties of tendon and aponeurosis of human gastrocnemius muscle in vivo. J Appl Physiol 2001, 90: 1671-1678.PubMedGoogle Scholar
- Hidler JM, Rymer WZ: Limit cycle behavior in spasticity: analysis and evaluation. IEEE Trans Biomed Eng 2000, 47: 1565-1575. 10.1109/10.887937View ArticlePubMedGoogle Scholar
- Kearney RE, Hunter IW: Dynamics of human ankle stiffness: variation with displacement amplitude. J Biomech 1982, 15: 753-756. 10.1016/0021-9290(82)90090-2View ArticlePubMedGoogle Scholar
- Hidler JM, Rymer WZ: A simulation study of reflex instability in spasticity: origins of clonus. IEEE Trans Rehabil Eng 1999, 7: 327-340. 10.1109/86.788469View ArticlePubMedGoogle Scholar
- Mirbagheri MM, Barbeau H, Kearney RE: Intrinsic and reflex contributions to human ankle stiffness: variation with activation level and position. Exp Brain Res 2000, 135: 423-436. 10.1007/s002210000534View ArticlePubMedGoogle Scholar
- de Vlugt E, Schouten AC, van der Helm FC: Quantification of intrinsic and reflexive properties during multijoint arm posture. J Neurosci Methods 2006, 155: 328-349. 10.1016/j.jneumeth.2006.01.022View ArticlePubMedGoogle Scholar
- Olney SJ, Winter DA: Predictions of knee and ankle moments of force in walking from EMG and kinematic data. J Biomech 1985, 18: 9-20. 10.1016/0021-9290(85)90041-7View ArticlePubMedGoogle Scholar
- Potvin JR, Norman RW, McGill SM: Mechanically corrected EMG for the continuous estimation of erector spinae muscle loading during repetitive lifting. Eur J Appl Physiol Occup Physiol 1996, 74: 119-132. 10.1007/BF00376504View ArticlePubMedGoogle Scholar
- Bobet J, Norman RW: Least-squares identification of the dynamic relation between the electromyogram and joint moment. J Biomech 1990, 23: 1275-1276. 10.1016/0021-9290(90)90386-HView ArticlePubMedGoogle Scholar
- Pisano F, Miscio G, Del Conte C, Pianca D, Candeloro E, Colombo R: Quantitative measures of spasticity in post-stroke patients. Clin Neurophysiol 2000, 111: 1015-1022. 10.1016/S1388-2457(00)00289-3View ArticlePubMedGoogle Scholar
- Sinkjaer T, Magnussen I: Passive, intrinsic and reflex-mediated stiffness in the ankle extensors of hemiparetic patients. Brain 1994,117(Pt 2):355-363. 10.1093/brain/117.2.355View ArticlePubMedGoogle Scholar
- Gorassini MA, Knash ME, Harvey PJ, Bennett DJ, Yang JF: Role of motoneurons in the generation of muscle spasms after spinal cord injury. Brain 2004, 127: 2247-2258. 10.1093/brain/awh243View ArticlePubMedGoogle Scholar
- Nielsen JB, Crone C, Hultborn H: The spinal pathophysiology of spasticity--from a basic science point of view. Acta Physiol (Oxf) 2007, 189: 171-180. 10.1111/j.1748-1716.2006.01652.xView ArticleGoogle Scholar
- Lamontagne A, Malouin F, Richards CL, Dumas F: Impaired viscoelastic behaviour of spastic plantarflexors during passive stretch at different velocities. Clin Biomech (Bristol, Avon) 1997, 12: 508-515. 10.1016/S0268-0033(97)00036-3View ArticleGoogle Scholar
- Bartoo ML, Linke WA, Pollack GH: Basis of passive tension and stiffness in isolated rabbit myofibrils. Am J Physiol 1997, 273: C266-276.PubMedGoogle Scholar
- Proske U, Morgan DL: Do cross-bridges contribute to the tension during stretch of passive muscle? J Muscle Res Cell Motil 1999, 20: 433-442. 10.1023/A:1005573625675View ArticlePubMedGoogle Scholar
- Maganaris CN: In vivo measurement-based estimations of the moment arm in the human tibialis anterior muscle-tendon unit. J Biomech 2000, 33: 375-379. 10.1016/S0021-9290(99)00188-8View ArticlePubMedGoogle Scholar
- Maganaris CN, Baltzopoulos V, Sargeant AJ: Changes in the tibialis anterior tendon moment arm from rest to maximum isometric dorsiflexion: in vivo observations in man. Clin Biomech (Bristol, Avon) 1999, 14: 661-666. 10.1016/S0268-0033(99)00018-2View ArticleGoogle Scholar
- Thelen DG: Adjustment of muscle mechanics model parameters to simulate dynamic contractions in older adults. J Biomech Eng 2003, 125: 70-77. 10.1115/1.1531112View ArticlePubMedGoogle Scholar
- Gajdosik RL: Passive extensibility of skeletal muscle: review of the literature with clinical implications. Clin Biomech (Bristol, Avon) 2001, 16: 87-101. 10.1016/S0268-0033(00)00061-9View ArticleGoogle Scholar
- Esteki A, Mansour JM: An experimentally based nonlinear viscoelastic model of joint passive moment. J Biomech 1996, 29: 443-450. 10.1016/0021-9290(95)00081-XView ArticlePubMedGoogle Scholar
- Genadry WF, Kearney RE, Hunter IW: Dynamic relationship between EMG and torque at the human ankle: variation with contraction level and modulation. Med Biol Eng Comput 1988, 26: 489-496. 10.1007/BF02441916View ArticlePubMedGoogle Scholar
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.