Katherine R. Saul

Benchmarking of dynamic simulation predictions in two software platforms using an upper limb musculoskeletal model

Several opensource or commercially available software platforms are widely used to develop dynamic simulations of movement. While computational approaches are conceptually similar across platforms, technical differences in implementation may influence output. We present a new upper limb dynamic model as a tool to evaluate potential differences in predictive behavior between platforms. We evaluated to what extent differences in technical implementations in popular simulation software environments result in differences in kinematic predictions for single and multijoint movements using EMG- and optimization-based approaches for deriving control signals. We illustrate the benchmarking comparison using SIMM–Dynamics Pipeline–SD/Fast and OpenSim platforms. The most substantial divergence results from differences in muscle model and actuator paths. This model is a valuable resource and is available for download by other researchers. The model, data, and simulation results presented here can be used by future researchers to benchmark other software platforms and software upgrades for these two platforms.

Characterizing upper limb muscle volume and strength in older adults: A comparison with young adults

Aging is associated with the loss of muscle volume (MV) and force leading to difficulties with activities of daily living. However, the relationship between upper limb MV and joint strength has not been characterized for older adults. Quantifying this relationship may help our understanding of the functional declines of the upper limb that older adults experience. Our objective was to assess the relationship between upper limb MV and maximal isometric joint moment-generating capacity (IJM) in a single cohort of healthy older adults (age ≥ 65 years) for 6 major functional groups (32 muscles). MV was determined from MRI for 18 participants (75.1±4.3 years). IJM at the shoulder (abduction/adduction), elbow (flexion/extension), and wrist (flexion/extension) was measured. MV and IJM measurements were compared to previous reports for young adults (28.6±4.5 years). On average older adults had 16.5% less total upper limb MV compared to young adults. Additionally, older adult wrist extensors composed a significantly increased percentage of upper limb MV. Older adult IJM was reduced across all joints, with significant differences for shoulder abductors (p<0.0001), adductors (p=0.01), and wrist flexors (p<0.0001). Young adults were strongest at the shoulder, which was not the case for older adults. In older adults, 40.6% of the variation in IJM was accounted for by MV changes (p≤0.027), compared to 81.0% in young adults. We conclude that for older adults, MV and IJM are, on average, reduced but the significant linear relationship between MV and IJM is maintained. These results suggest that older adult MV and IJM cannot be simply scaled from young adults.

Postural dependence of passive tension in the supraspinatus following rotator cuff repair: A simulation analysis☆

BACKGROUND Despite surgical advances, repair of rotator cuff tears is associated with 20-70% incidence of recurrent tearing. The tension required to repair the torn tendon influences surgical outcomes and may be dependent on the gap length from torn tendon that must be spanned by the repair. Detailed understanding of forces throughout the range of motion (ROM) may allow surgeons to make evidence-based recommendations for post-operative care. METHODS We used a computational shoulder model to assess passive tension and total moment-generating capacity in supraspinatus for repairs of gaps up to 3 cm throughout the shoulder (ROM). FINDINGS In 60° abduction, increased gap length from 0.5 cm to 3 cm caused increases in passive force from 3N to 58 N, consistent with those seen during clinical repair. For reduced abduction, passive forces increased substantially. For a 0.5 cm gap, tension throughout the ROM (elevation, plane of elevation, and rotation) is within reasonable limits, but larger gaps are associated with tensions that markedly exceed reported pull-out strength of sutures and anchors. Peak moment for a large 3 cm gap length was 5.09 Nm, a 53% reduction in moment-generating capacity compared to uninjured supraspinatus. INTERPRETATION We conclude that shoulder posture is an important determinant of passive forces during rotator cuff repair surgery. Choosing postures that reduce forces intraoperatively to permit repair of larger gaps may lead to failure postoperatively when the shoulder is mobilized. For larger defects, loss of strength in supraspinatus may be substantial following repair even if retear is prevented.

Biomechanics of the Steindler Flexorplasty Surgery: A Computer Simulation Study

PURPOSE Our goal was to investigate the capacity of a Steindler flexorplasty to restore elbow flexion to persons with C5-C6 brachial plexus palsy. In this procedure the origin of the flexor-pronator mass is moved proximally onto the humeral shaft. We examined how the choice of the proximal attachment site for the flexor-pronator mass affects elbow flexion restoration, especially considering possible side effects including limited wrist and forearm motion owing to passive restraint from stretched muscles. METHODS A computer model of the upper extremity was used to simulate the biomechanical consequences of various surgical alterations. Unimpaired, preoperative, and postoperative conditions were simulated. Seven possible transfer locations were used to investigate the effects of choice of transfer location. RESULTS Each transfer site produced a large increase in elbow flexion strength. Transfer to more proximal attachment sites also produced large increases in passive resistance to wrist extension and forearm supination. CONCLUSIONS To reduce detrimental side effects while achieving clinical goals our theoretical analysis suggests a transfer to the distal limit of the traditional transfer region.

Biomechanical Contributions of Posterior Deltoid and Teres Minor in the Context of Axillary Nerve Injury: A Computational Study

PURPOSE To determine whether transfer to only the anterior branch of the axillary nerve will restore useful function after axillary nerve injury with persistent posterior deltoid and teres minor paralysis. METHODS We used a computational musculoskeletal model of the upper limb to determine the relative contributions of posterior deltoid and teres minor to maximum joint moment generated during a simulated static strength assessment and to joint moments during 3 submaximal shoulder movements. Movement simulations were performed with and without simulated posterior deltoid and teres minor paralysis to identify muscles that may compensate for their paralysis. RESULTS In the unimpaired limb model, teres minor and posterior deltoid accounted for 16% and 14% of the total isometric shoulder extension and external rotation joint moments, respectively. During the 3 movement simulations, posterior deltoid produced as much as 20% of the mean shoulder extension moment, whereas teres minor accounted for less than 5% of the mean joint moment in all directions of movement. When we paralyzed posterior deltoid and teres minor, the mean extension moments generated by the supraspinatus, long head of triceps, latissimus dorsi, and middle deltoid increased to compensate. Compensatory muscles were not fully activated during movement simulations when posterior deltoid and teres minor were paralyzed. CONCLUSIONS Reconstruction of the anterior branch of the axillary nerve only is an appropriate technique for restoring shoulder abduction strength after isolated axillary nerve injury. When shoulder extension strength is compromised by extensive neuromuscular shoulder injury, reconstruction of both the anterior and posterior branches of the axillary nerve should be considered. CLINICAL RELEVANCE By quantifying the biomechanical role of muscles during submaximal movement, in addition to quantifying muscle contributions to maximal shoulder strength, we can inform preoperative planning and permit more accurate predictions of functional outcomes.

A simulation analysis of the combined effects of muscle strength and surgical tensioning on lateral pinch force following brachioradialis to flexor pollicis longus transfer

Biomechanical simulations of tendon transfers performed following tetraplegia suggest that surgical tensioning influences clinical outcomes. However, previous studies have focused on the biomechanical properties of only the transferred muscle. We developed simulations of the tetraplegic upper limb following transfer of the brachioradialis (BR) to the flexor pollicis longus (FPL) to examine the influence of residual upper limb strength on predictions of post-operative transferred muscle function. Our simulations included the transfer, ECRB, ECRL, the three heads of the triceps, brachialis, and both heads of the biceps. Simulations were integrated with experimental data, including EMG and joint posture data collected from five individuals with tetraplegia and BR-FPL tendon transfers during maximal lateral pinch force exertions. Given a measured co-activation pattern for the non-paralyzed muscles in the tetraplegic upper limb, we computed the highest activation for the transferred BR for which neither the elbow nor the wrist flexor moment was larger than the respective joint extensor moment. In this context, the effects of surgical tensioning were evaluated by comparing the resulting pinch force produced at different muscle strength levels, including patient-specific scaling. Our simulations suggest that extensor muscle weakness in the tetraplegic limb limits the potential to augment total pinch force through surgical tensioning. Incorporating patient-specific muscle volume, EMG activity, joint posture, and strength measurements generated simulation results that were comparable to experimental results. Our study suggests that scaling models to the population of interest facilitates accurate simulation of post-operative outcomes, and carries utility for guiding and developing rehabilitation training protocols.

Normal aging alters in vivo passive biomechanical response of the rat gastrocnemius-Achilles muscle-tendon unit

Predisposition to Achilles tendon (AT) ruptures in middle-aged individuals may be associated with age-related changes to inherent passive biomechanical properties of the gastrocnemius-Achilles (GC-AT) muscle-tendon unit, due to known muscle-tendon structural changes in normal aging. The goal of this study was to determine whether the passive biomechanical response of the GC-AT muscle-tendon unit was altered with age in 6 young (8 months) and 6 middle-aged (24 months) F344xBN hybrid rats from the National Institute on Aging colony. Fungs quasilinear viscoelastic (QLV) model was used to determine in vivo history and time-dependent load-relaxation response of the GC-AT. Effective stiffness and modulus were also estimated using linear regression analysis. Fungs QLV revealed a significantly decreased magnitude of the relaxation response (parameter C, p=0.026) in middle-aged animals compared to young animals (0.108±0.007 vs. 0.144±0.015), with similar time-dependent viscous GC-AT properties (τ(1), τ(2)). The product of elastic parameters (A*B), which represents the initial slope of the elastic response, was significantly increased by 50% in middle-aged rats (p=0.014). Estimated GC-AT stiffness increased 28% at peak tensions in middle-aged rats (2.7±0.2 N/mm) compared to young rats (1.9±0.2 N/mm; p=0.036). While the limitations of this animal model must be considered, the changes we describe could be associated with the observation that GC-AT pathology and injury is more common in middle-aged individuals. Further studies are necessary to characterize the load-to-failure behavior of AT in middle-aged compared to young animals.

Computational sensitivity analysis to identify muscles that can mechanically contribute to shoulder deformity following brachial plexus birth palsy.

PURPOSE Two mechanisms, strength imbalance or impaired longitudinal muscle growth, potentially cause osseous and postural shoulder deformity in children with brachial plexus birth palsy. Our objective was to determine which muscles, via either deformity mechanism, were mechanically capable of producing forces that could promote shoulder deformity. METHODS In an upper limb computational musculoskeletal model, we simulated strength imbalance by allowing each muscle crossing the shoulder to produce 30% of its maximum force. To simulate impaired longitudinal muscle growth, the functional length of each muscle crossing the shoulder was reduced by 30%. We performed a sensitivity analysis to identify muscles that, through either simulated deformity mechanism, increased the posteriorly directed, compressive glenohumeral joint force consistent with osseous deformity or reduced the shoulder external rotation or abduction range of motion consistent with postural deformity. RESULTS Most of the increase in the posterior glenohumeral joint force by the strength imbalance mechanism was caused by the subscapularis, latissimus dorsi, and infraspinatus. Posterior glenohumeral joint force increased the most owing to impaired growth of the infraspinatus, subscapularis, and long head of biceps. Through the strength imbalance mechanism, the subscapularis, anterior deltoid, and pectoralis major muscles reduced external shoulder rotation by 28°, 17°, and 10°, respectively. Shoulder motion was reduced by 40° to 56° owing to impaired growth of the anterior deltoid, subscapularis, and long head of triceps. CONCLUSIONS The infraspinatus, subscapularis, latissimus dorsi, long head of biceps, anterior deltoid, pectoralis major, and long head of triceps were identified in this computational study as being the most capable of producing shoulder forces that may contribute to shoulder deformity following brachial plexus birth palsy. CLINICAL RELEVANCE The muscles mechanically capable of producing deforming shoulder forces should be the focus of experimental studies investigating the musculoskeletal consequences of brachial plexus birth palsy and are potentially critical targets for treating shoulder deformity.

The science of rotator cuff tears: translating animal models to clinical recommendations using simulation analysis

PurposeThe purpose of this article is to review basic science studies using various animal models for rotator cuff research and to describe structural, biomechanical, and functional changes to muscle following rotator cuff tears. The use of computational simulations to translate the findings from animal models to human scale is further detailed.MethodsA comprehensive review was performed of the basic science literature describing the use of animal models and simulation analysis to examine muscle function following rotator cuff injury and repair in the ageing population.ResultsThe findings from various studies of rotator cuff pathology emphasize the importance of preventing permanent muscular changes with detrimental results. In vivo muscle function, electromyography, and passive muscle–tendon unit properties were studied before and after supraspinatus tenotomy in a rodent rotator cuff injury model (acute vs chronic). Then, a series of simulation experiments were conducted using a validated computational human musculoskeletal shoulder model to assess both passive and active tension of rotator cuff repairs based on surgical positioning.ConclusionOutcomes of rotator cuff repair may be improved by earlier surgical intervention, with lower surgical repair tensions and fewer electromyographic neuromuscular changes. An integrated approach of animal experiments, computer simulation analyses, and clinical studies may allow us to gain a fundamental understanding of the underlying pathology and interpret the results for clinical translation.

Computer Simulation of Nerve Transfer Strategies for Restoring Shoulder Function After Adult C5 and C6 Root Avulsion Injuries

PURPOSE Functional ability after nerve transfer for upper brachial plexus injuries relies on both the function and magnitude of force recovery of targeted muscles. Following nerve transfers targeting either the axillary nerve, suprascapular nerve, or both, it is unclear whether functional ability is restored in the face of limited muscle force recovery. METHODS We used a computer model to simulate flexing the elbow while maintaining a functional shoulder posture for 3 nerve transfer scenarios. We assessed the minimum restored force capacity necessary to perform the task, the associated compensations by neighboring muscles, and the effect of altered muscle coordination on movement effort. RESULTS The minimum force restored by the axillary, suprascapular, and combined nerve transfers that was required for the model to simulate the desired movement was 25%, 40%, and 15% of the unimpaired muscle force capacity, respectively. When the deltoid was paralyzed, the infraspinatus and subscapularis muscles generated higher shoulder abduction moments to compensate for deltoid weakness. For all scenarios, movement effort increased as restored force capacity decreased. CONCLUSIONS Combined axillary and suprascapular nerve transfer required the least restored force capacity to perform the desired elbow flexion task, whereas single suprascapular nerve transfer required the most restored force capacity to perform the same task. Although compensation mechanisms allowed all scenarios to perform the desired movement despite weakened shoulder muscles, compensation increased movement effort. Dynamic simulations allowed independent evaluation of the effect of restored force capacity on functional outcome in a way that is not possible experimentally. CLINICAL RELEVANCE Simultaneous nerve transfer to suprascapular and axillary nerves yields the best simulated biomechanical outcome for lower magnitudes of muscle force recovery in this computer model. Axillary nerve transfer performs nearly as well as the combined transfer, whereas suprascapular nerve transfer is more sensitive to the magnitude of reinnervation and is therefore avoided.