12 Degrees of Freedom Muscle Force Driven Fibril-Reinforced Poroviscoelastic Finite Element Model of the Knee Joint.

Accurate knowledge of the joint kinematics, kinetics, and soft tissue mechanical responses is essential in the evaluation of musculoskeletal (MS) disorders. Since in vivo measurement of these quantities requires invasive methods, musculoskeletal finite element (MSFE) models are widely used for simulations. There are, however, limitations in the current approaches. Sequentially linked MSFE models benefit from complex MS and FE models; however, MS model's outputs are independent of the FE model calculations. On the other hand, due to the computational burden, embedded (concurrent) MSFE models are limited to simple material models and cannot estimate detailed responses of the soft tissue. Thus, first we developed a MSFE model of the knee with a subject-specific MS model utilizing an embedded 12 degrees of freedom (DoFs) knee joint with elastic cartilages in which included both secondary kinematic and soft tissue deformations in the muscle force estimation (inverse dynamics). Then, a muscle-force-driven FE model with fibril-reinforced poroviscoelastic cartilages and fibril-reinforced poroelastic menisci was used in series to calculate detailed tissue mechanical responses (forward dynamics). Second, to demonstrate that our workflow improves the simulation results, outputs were compared to results from the same FE models which were driven by conventional MS models with a 1 DoF knee, with and without electromyography (EMG) assistance. The FE model driven by both the embedded and the EMG-assisted MS models estimated similar results and consistent with experiments from literature, compared to the results estimated by the FE model driven by the MS model with 1 DoF knee without EMG assistance.

Botulinum toxin injections minimally affect modelled muscle forces during gait in children with cerebral palsy.

BACKGROUND: Children with cerebral palsy (CP) present altered gait patterns and electromyography (EMG) activity compared to typically developing children. To temporarily reduce muscular activity and to correct the abnormal muscle force balance, Botulinum Toxin type A (BTX-A) injections are used. RESEARCH QUESTION: What is the effect of BTX-A injections on dynamic muscle forces during gait, when calculated using an EMG-constrained approach?. METHODS: Retrospective data of ten typically developing (TD) and fourteen children with spastic diplegic CP were used for musculoskeletal modeling and dynamic simulations of gait, before and after BTX-A treatment. Individual muscle forces were calculated using an EMG-constrained optimization, in which EMG of eight muscles was used as muscle excitation signal to constrain the muscle activation patterns. Paired t-tests were used to compare average modelled muscle forces in different phases of the gait cycle pre- and post-BTX-A, summarized in the muscle profile score. Two-sample t-tests were used to determine significant differences between TD and pre- and post-BTX-A modelled muscle forces. RESULTS: For most muscles, the force was decreased in CP compared to TD children in all phases of the gait cycle, both before and after BTX-A treatment. Differences in muscle forces before and after BTX-A treatment were limited, with only few significant differences between pre- and post-BTX-A. Compared to a standard static optimization approach, imposing the EMG activity increased modelled muscle forces for most muscles. SIGNIFICANCE: Our findings indicate that BTX-A treatment has a limited effect on the muscle balance in CP children. Besides that, the use of EMG-constrained optimization is recommended when studying muscle balance in children with CP.

Single-event multilevel surgery, but not botulinum toxin injections normalize joint loading in cerebral palsy patients.

BACKGROUND: Many patients with cerebral palsy present a pathologic gait pattern, which presumably induces aberrant musculoskeletal loading that interferes with natural bone growth, causing bone deformations on the long term. Botulinum toxin interventions and single-event multilevel surgeries are used to restore the gait pattern, assuming that a normal gait pattern restores musculoskeletal loading and thus prevents further bone deformation. However, it is unknown if these interventions are able to restore musculoskeletal loading. Hence, we investigated the impact of botulinum toxin injections and single-event multilevel surgery on musculoskeletal loading. METHODS: Gait data collected in 93 children with bilateral cerebral palsy, which included pre- and post multi-level botulinum toxin (49 children) and single-event multilevel surgery (44 children) assessments, and 15 typically developing children were retrospectively processed using a musculoskeletal modelling workflow to calculate joint angles, moments, muscle and joint contact force magnitudes and orientations. Differences from the typically developing waveform were expressed by a root-mean square difference were compared using paired t-tests for each intervention separately (alpha <0.05). FINDINGS: Botulinum toxin induced significant changes in the joint angles, but did not improve the muscle and joint contact forces. Single-event multilevel surgery induced significant kinematic and kinetic changes, which were associated with improved muscle and joint contact forces. INTERPRETATION: The present results indicate that botulinum toxin injections were not able to restore normal gait kinematics nor musculoskeletal loading, whereas single-event multilevel surgery did successfully restore both. Therefore, single-event multilevel surgery might be protective against the re-occurrence of bone deformation on the longer term.

Differences in knee adduction moment between healthy subjects and patients with osteoarthritis depend on the knee axis definition.

OBJECTIVE: This study, firstly, investigates the effect of using an anatomical versus a functional axis of rotation (FAR) on knee adduction moment (KAM) in healthy subjects and patients with knee osteoarthritis (KOA). Secondly, this study reports KAM for models with FAR calculated using weight-bearing and non-weight-bearing motion. DESIGN: Three musculoskeletal models were created using OpenSim with different knee axis of rotation (AR): transepicondylar axis (TEA); FAR calculated based on SARA algorithm using a weight-bearing motion (wFAR) and a non-weight-bearing motion (nwFAR). KAM were calculated during gait in fifty-nine subjects (n=20 healthy, n=16 early OA, n=23 established OA) for all models and groups. RESULTS: Significant differences between the three groups in the first peak KAM were found when TEA was used (p=0.038). However, these differences were no longer present when using FAR. In subjects with established OA, KAMs were significantly reduced when using nwFAR compared to TEA models but also compared to wFAR models. CONCLUSION: The presence of excessive KAM in subjects with established KOA showed to be dependent on the definition of the AR: anatomical versus functional. Therefore, caution should be accounted when comparing KAM in different studies on KOA patients. In patients with end-stage knee OA where increased passive knee laxity is likely to exist, the use of weight-bearing motions should be considered to avoid increased variability in the location and orientation of a FAR obtained from activities with only limited joint loading.