Benjamin R. Shuman

Flow mechanotransduction regulates traction forces, intercellular forces, and adherens junctions

Endothelial cells respond to fluid shear stress through mechanotransduction responses that affect their cytoskeleton and cell-cell contacts. Here, endothelial cells were grown as monolayers on arrays of microposts and exposed to laminar or disturbed flow to examine the relationship among traction forces, intercellular forces, and cell-cell junctions. Cells under laminar flow had traction forces that were higher than those under static conditions, whereas cells under disturbed flow had lower traction forces. The response in adhesion junction assembly matched closely with changes in traction forces since adherens junctions were larger in size for laminar flow and smaller for disturbed flow. Treating the cells with calyculin-A to increase myosin phosphorylation and traction forces caused an increase in adherens junction size, whereas Y-27362 cause a decrease in their size. Since tugging forces across cell-cell junctions can promote junctional assembly, we developed a novel approach to measure intercellular forces and found that these forces were higher for laminar flow than for static or disturbed flow. The size of adherens junctions and tight junctions matched closely with intercellular forces for these flow conditions. These results indicate that laminar flow can increase cytoskeletal tension while disturbed flow decreases cytoskeletal tension. Consequently, we found that changes in cytoskeletal tension in response to shear flow conditions can affect intercellular tension, which in turn regulates the assembly of cell-cell junctions.

Repeatability of muscle synergies within and between days for typically developing children and children with cerebral palsy

Muscle synergies are typically calculated from electromyographic (EMG) signals using nonnegative matrix factorization. Synergies identify weighted groups of muscles that are commonly activated together during a task, such as walking. Synergy analysis has become an emerging tool to evaluate neuromuscular control; however, the repeatability of synergies between trials and days has not been evaluated. The goal of this study was to evaluate the repeatability of synergy complexity and structure in unimpaired individuals and individuals with cerebral palsy (CP). EMG data were collected from eight lower-limb muscles during gait for six typically developing (TD) children and five children with CP on two separate days, over three walking speeds. To evaluate synergy complexity, we calculated the total variance accounted for by one synergy (tVAF1). On a given day, the average range in tVAF1 between gait cycles was 18.2% for TD and 19.1% for CP. The average standard deviation in tVAF1 between gait cycles was 4.9% for TD and 5.0% for CP. Average tVAF1 calculated across gait cycles was not significantly different between days for TD or CP participants. Comparing synergy structure, the average (standard deviation) within day correlation coefficients of synergy weights for two or more synergies were 0.89 (0.15) for TD and 0.88 (0.15) for CP. Between days, the average correlation coefficient of synergy weights for two or more synergies was greater than 0.89 for TD and 0.74 for CP. These results demonstrate that synergy complexity and structure averaged over multiple gait cycles are repeatable between days in both TD and CP groups.

Electromyography Data Processing Impacts Muscle Synergies during Gait for Unimpaired Children and Children with Cerebral Palsy

Muscle synergies calculated from electromyography (EMG) data identify weighted groups of muscles activated together during functional tasks. Research has shown that fewer synergies are required to describe EMG data of individuals with neurologic impairments. When considering potential clinical applications of synergies, understanding how EMG data processing impacts results and clinical interpretation is important. The aim of this study was to evaluate how EMG signal processing impacts synergy outputs during gait. We evaluated the impacts of two common processing steps for synergy analyses: low pass (LP) filtering and unit variance scaling. We evaluated EMG data collected during barefoot walking from five muscles of 113 children with cerebral palsy (CP) and 73 typically-developing (TD) children. We applied LP filters to the EMG data with cutoff frequencies ranging from 4 to 40 Hz (reflecting the range reported in prior synergy research). We also evaluated the impact of normalizing EMG amplitude by unit variance. We found that the total variance accounted for (tVAF) by a given number of synergies was sensitive to LP filter choice and decreased in both TD and CP groups with increasing LP cutoff frequency (e.g., 9.3 percentage points change for one synergy between 4 and 40 Hz). This change in tVAF can alter the number of synergies selected for further analyses. Normalizing tVAF to a z-score (e.g., dynamic motor control index during walking, walk-DMC) reduced sensitivity to LP cutoff. Unit variance scaling caused comparatively small changes in tVAF. Synergy weights and activations were impacted less than tVAF by LP filter choice and unit variance normalization. These results demonstrate that EMG signal processing methods impact outputs of synergy analysis and z-score based measures can assist in reporting and comparing results across studies and clinical centers.

Muscle recruitment and coordination with an ankle exoskeleton

Exoskeletons have the potential to assist and augment human performance. Understanding how users adapt their movement and neuromuscular control in response to external assistance is important to inform the design of these devices. The aim of this research was to evaluate changes in muscle recruitment and coordination for ten unimpaired individuals walking with an ankle exoskeleton. We evaluated changes in the activity of individual muscles, cocontraction levels, and synergistic patterns of muscle coordination with increasing exoskeleton work and torque. Participants were able to selectively reduce activity of the ankle plantarflexors with increasing exoskeleton assistance. Increasing exoskeleton net work resulted in greater reductions in muscle activity than increasing exoskeleton torque. Patterns of muscle coordination were not restricted or constrained to synergistic patterns observed during unassisted walking. While three synergies could describe nearly 95% of the variance in electromyography data during unassisted walking, these same synergies could describe only 85-90% of the variance in muscle activity while walking with the exoskeleton. Synergies calculated with the exoskeleton demonstrated greater changes in synergy weights with increasing exoskeleton work versus greater changes in synergy activations with increasing exoskeleton torque. These results support the theory that unimpaired individuals do not exclusively use central pattern generators or other low-level building blocks to coordinate muscle activity, especially when learning a new task or adapting to external assistance, and demonstrate the potential for using exoskeletons to modulate muscle recruitment and coordination patterns for rehabilitation or performance.

Repeatability of electromyography recordings and muscle synergies during gait among children with cerebral palsy

BACKGROUND Clinical gait analysis is commonly used in the evaluation and treatment of children with cerebral palsy (CP). While the repeatability of kinematic and kinetic measures of gait has previously been evaluated, the repeatability of electromyography (EMG) recordings or measures calculated from EMG data, such as muscle synergies, remains unclear for this population. RESEARCH QUESTION Are EMG recordings and muscle synergies from clinical gait analysis repeatable between visits for children with CP? METHODS We recruited 20 children with bilateral CP who had been referred for clinical gait analysis. The children completed two visits less than six weeks apart with EMG data collected bilaterally from five muscles (rectus femoris, medial hamstrings, vastus lateralis, anterior tibialis, and medial gastrocnemius). Variance ratio and cosine similarity were used to evaluate repeatability of EMG waveforms between visits. Nonnegative matrix factorization was used to calculate synergies from EMG data at each visit to compare synergy weights and activations. RESULTS & SIGNIFICANCE The inter-visit variance ratios of EMG data for children with CP were similar to previously reported results for typically-developing children and unimpaired adults (range: 0.39 for vastus lateralis to 0.66 for rectus femoris). The average cosine similarity of the EMG waveforms between visits was greater than 0.9 for all muscles, while synergy weights and activations also had high similarity - greater than 0.8 and 0.9 between visits, respectively. These results demonstrate that EMG repeatability between visits during clinical gait analysis for children with CP is similar to unimpaired individuals. These results provide a baseline for evaluating whether observed changes in EMG recordings between visits reflect real changes in muscle activity or are within the range of inter-visit variability.

Non-neural Muscle Weakness Has Limited Influence on Complexity of Motor Control during Gait

Cerebral palsy (CP) and Duchenne muscular dystrophy (DMD) are neuromuscular disorders characterized by muscle weakness. Weakness in CP has neural and non-neural components, whereas in DMD, weakness can be considered as a predominantly non-neural problem. Despite the different underlying causes, weakness is a constraint for the central nervous system when controlling gait. CP demonstrates decreased complexity of motor control during gait from muscle synergy analysis, which is reflected by a higher total variance accounted for by one synergy (tVAF1). However, it remains unclear if weakness directly contributes to higher tVAF1 in CP, or whether altered tVAF1 reflects mainly neural impairments. If muscle weakness directly contributes to higher tVAF1, then tVAF1 should also be increased in DMD. To examine the etiology of increased tVAF1, muscle activity data of gluteus medius, rectus femoris, medial hamstrings, medial gastrocnemius, and tibialis anterior were measured at self-selected walking speed, and strength data from knee extensors, knee flexors, dorsiflexors and plantar flexors, were analyzed in 15 children with CP [median (IQR) age: 8.9 (2.2)], 15 boys with DMD [8.7 (3.1)], and 15 typical developing (TD) children [8.6 (2.7)]. We computed tVAF1 from 10 concatenated steps with non-negative matrix factorization, and compared tVAF1 between the three groups with a Mann-Whiney U-test. Spearmans rank correlation coefficients were used to determine if weakness in specific muscle groups contributed to altered tVAF1. No significant differences in tVAF1 were found between DMD [tVAF1: 0.60 (0.07)] and TD children [0.65 (0.07)], while tVAF1 was significantly higher in CP [(0.74 (0.09)] than in the other groups (both p < 0.005). In CP, weakness in the plantar flexors was related to higher tVAF1 (r = −0.72). In DMD, knee extensor weakness related to increased tVAF1 (r = −0.50). These results suggest that the non-neural weakness in DMD had limited influence on complexity of motor control during gait and that the higher tVAF1 in children with CP is mainly related to neural impairments caused by the brain lesion.

Associations between muscle synergies and treatment outcomes in cerebral palsy are robust across clinical centers

OBJECTIVE To determine whether patient-specific differences in motor control quantified using muscle synergy analysis were associated with changes in gait after treatment of cerebral palsy (CP) across 2 clinical centers with different treatments and clinical protocols. DESIGN Retrospective cohort study. SETTING Clinical medical center. PARTICIPANTS Center 1: children with CP (n=473) and typically developing (TD) children (n=84). Center 2: children with CP (n=163) and TD children (n=12). INTERVENTIONS Standard clinical care at each center. MAIN OUTCOME MEASURES The Dynamic Motor Control Index During Walking (walk-DMC) was computed from electromyographic data during gait using muscle synergy analysis. Regression analysis was used to evaluate whether pretreatment walking speed or kinematics, muscle synergies, treatment group, prior treatment, or age were associated with posttreatment changes in gait at both clinical centers. RESULTS Walk-DMC was significantly associated with changes in speed and kinematics after treatment with similar regression models at both centers. Children with less impaired motor control were more likely to have improvements in walking speed and gait kinematics after treatment, independent of treatment group. CONCLUSIONS Dynamic motor control evaluated with synergy analysis was associated with changes in gait after treatment at both centers, despite differences in treatments and clinical protocols. This study further supports the finding that walk-DMC provides additional information, not captured in traditional gait analysis, that may be useful for treatment planning.

Crouch severity is a poor predictor of elevated oxygen consumption in cerebral palsy

Children with cerebral palsy (CP) expend more energy to walk compared to typically-developing peers. One of the most prevalent gait patterns among children with CP, crouch gait, is often singled out as especially exhausting. The dynamics of crouch gait increase external flexion moments and the demand on extensor muscles. This elevated demand is thought to dramatically increase energy expenditure. However, the impact of crouch severity on energy expenditure has not been investigated among children with CP. We evaluated oxygen consumption and gait kinematics for 573 children with bilateral CP. The average net nondimensional oxygen consumption during gait of the children with CP (0.18±0.06) was 2.9 times that of speed-matched typically-developing peers. Crouch severity was only modestly related to oxygen consumption, with measures of knee flexion angle during gait explaining only 5-20% of the variability in oxygen consumption. While knee moment and muscle activity were moderately to strongly correlated with crouch severity (r2=0.13-0.73), these variables were only weakly correlated with oxygen consumption (r2=0.02-0.04). Thus, although the dynamics of crouch gait increased muscle demand, these effects did not directly result in elevated energy expenditure. In clinical gait analysis, assumptions about an individuals energy expenditure should not be based upon kinematics or kinetics alone. Identifying patient-specific factors that contribute to increased energy expenditure may provide new pathways to improve gait for children with CP.

Multistep model for predicting upper-limb 3D isometric force application from pre-movement electrocorticographic features

Neural correlates of movement planning onset and direction may be present in human electrocorticography in the signal dynamics of both motor and non-motor cortical regions. We use a three-stage model of jPCA reduced-rank hidden Markov model (jPCA-RR-HMM), regularized shrunken-centroid discriminant analysis (RDA), and LASSO regression to extract direction-sensitive planning information and movement onset in an upper-limb 3D isometric force task in a human subject. This mode achieves a relatively high true positive force-onset prediction rate of 60% within 250ms, and an above-chance 36% accuracy (17% chance) in predicting one of six planned 3D directions of isometric force using pre-movement signals. We also find direction-distinguishing information up to 400ms before force onset in the pre-movement signals, captured by electrodes placed over the limb-ipsilateral dorsal premotor regions. This approach can contribute to more accurate decoding of higher-level movement goals, at earlier timescales, and inform sensor placement. Our results also contribute to further understanding of the spatiotemporal features of human motor planning.Neural correlates of movement planning onset and direction may be present in human electrocorticography in the signal dynamics of both motor and non-motor cortical regions. We use a three-stage model of jPCA reduced-rank hidden Markov model (jPCA-RR-HMM), regularized shrunken-centroid discriminant analysis (RDA), and LASSO regression to extract direction-sensitive planning information and movement onset in an upper-limb 3D isometric force task in a human subject. This mode achieves a relatively high true positive force-onset prediction rate of 60% within 250ms, and an above-chance 36% accuracy (17% chance) in predicting one of six planned 3D directions of isometric force using pre-movement signals. We also find direction-distinguishing information up to 400ms before force onset in the pre-movement signals, captured by electrodes placed over the limb-ipsilateral dorsal premotor regions. This approach can contribute to more accurate decoding of higher-level movement goals, at earlier timescales, and inform sensor placement. Our results also contribute to further understanding of the spatiotemporal features of human motor planning.