Background: Osteoarthritis is a highly prevalent disease affecting the hip and knee joint and is characterized by load-mediated pain and decreased quality of life. Dependent on involved joint, patients present antalgic movement compensations, aiming to decrease loading on the involved joint. However, the associated alterations in mechanical loading of the ipsi- and contra-lateral lower limb joints, are less documented. Here, we documented the biomechanical fingerprint of end-stage hip and knee osteoarthritis patients in terms of ipsilateral and contralateral hip and knee loading during walking and stair ambulation. Methods: Three-dimensional motion-analysis was performed in 20 hip, 18 knee osteoarthritis patients and 12 controls during level walking and stair ambulation. Joint contact forces were calculated using a standard musculoskeletal modelling workflow in Opensim. Involved and contralateral hip and knee joint loading was compared against healthy controls using independent t-tests (p < 0.05). Findings: Both hip and knee cohorts significantly decreased loading of the involved joint during gait and stair ambulation. Hip osteoarthritis patients presented no signs of ipsilateral knee nor contralateral leg overloading, during walking and stair ascending. However, knee osteoarthritis patients significantly increased loading at the ipsilateral hip, and contralateral hip and knee joints during stair ambulation compared to controls. Interpretation: The biomechanical fingerprint in knee and hip osteoarthritis patients confirmed antalgic movement strategies to unload the involved leg during gait. Only during stair ambulation in knee osteoarthritis patients, movement adaptations were confirmed that induced unbalanced intra- and inter-limb loading conditions, which are known risk factors for secondary osteoarthritis.

Current fixation plates used to operatively stabilize clavicular fractures are suboptimal, leading to reoperation rates up to 53%. Plate irritation, which can be caused by a poor geometric fit and plate thickness, has been found to be an important factor for reoperation. Moreover, muscle attachment sites (MAS) have to be detached occasionally. To improve current surgical clavicle fracture treatment with plate osteosynthesis, a change in plate design is required. The goal of this study was to design a patient-specific clavicle fracture fixation plate that includes geometrical optimization and stiffness determination. Its biomechanical performance has been compared with a commercial plate by examining the geometric fit, anatomical outline, stresses and interfragmentary motion using a finite element analysis with physiological loading and boundary conditions. Evaluation showed a better geometrical fit of the patient-specific plate as well as an improved fracture reduction. Displacements between fracture fragments were lower in case of the patient-specific plate, both when a fracture gap and no fracture gap were present. Results indicate a superior mechanical performance in terms of all investigated outcomes of the patient-specific plate compared to the commercial plate, while better aligning with the patient-specific geometry and without the need for MAS release. Due to the patient-specific geometry and reduced thickness, these fixation plates are expected to decrease the operation time, as intraoperative contouring will become irrelevant, and to decrease reoperation rates as implant irritation will be minimized.

Mechanical performance of clavicle fracture fixation plates is often evaluated using finite element (FE) analysis. Typically, these studies use simplified loading conditions and assume a transversal fracture orientation. However, the loading conditions and fracture orientation influence how the fracture site and thus fixation plate is loaded. In this study, a musculoskeletal model that included the clavicle muscles and scapulohumeral rhythm was defined based on previously published models. The standard OpenSim workflow (inverse kinematics, inverse dynamics, static optimization, and joint reaction analysis) was used to calculate muscle and joint contact forces based on 3D marker data collected in three subjects during seven activities of daily living (ADL). These loading conditions were then applied to a 3D clavicle model with three different fracture orientations and the mean resulting moments on both fragments were calculated to assess fracture stability. Magnitude of glenohumeral contact forces showed good agreement with instrumented shoulder prosthesis data, whereas simulated muscle activations were comparable to experimental EMG data. An oblique fracture orienting from superomedial to inferolateral was the least self-stabilizing. The loading to which the clavicle is exposed during ADL tasks is more complex than the simplified loading conditions typically used as boundary conditions in FE analyses of clavicle fracture fixation plates. Additionally, transversal fractures did not represent the least self-stabilizing fracture orientation, and thus calculated stresses in the plate could be underestimated. Therefore, more complex loading conditions and evaluation of a midshaft fracture running from superomedial to inferolateral is more relevant in FE analyses of midshaft clavicle fracture fixation plates.

Objective: This study proposes a computationally efficient method to quantify the effect of surgical inaccuracies on ligament strain in total knee arthroplasty (TKA). More specifically, this study describes a framework to determine the implant position and required surgical accuracy that results in a ligament balanced post-operative outcome with a probability of 90%. Methods: The response surface method is used to translate uncertainty in the implant position parameters to uncertainty in the ligament strain. The designed uncertainty quantification technique allows for an optimization with feasible computational cost towards the planned implant position and the tolerated surgical error for each of the twelve degrees of freedom of the implant position. Results: It is shown that the error does not allow for a ligament balanced TKA with a probability of 90% using preoperative planning. Six critical implant position parameters can be identified, namely AP translation, PD translation, VV rotation, IE rotation for the femoral component and PD translation, VV rotation for the tibial component. Conclusion: We introduced an optimization process that allows for the computation of the required surgical accuracy for a ligament balanced postoperative outcome using preoperative planning with feasible computational cost. Significance: Towards the research society, the proposed method allows for a computationally efficient uncertainty quantification on a complex model. Towards surgical technique developers, six critical implant position parameters were identified, which should be the focus when refining surgical accuracy of TKA, leveraging better patient satisfaction.

Background: Musculoskeletal modelling is used to assess musculoskeletal loading during gait. Linear scaling methods are used to personalize generic models to each participant's anthropometry. This approach introduces simplifications, especially when used in paediatric and/or pathological populations. This study aimed to compare results from musculoskeletal simulations using various models ranging from linear scaled to highly subject-specific models, i.e., including the participant's musculoskeletal geometry and electromyography data. Methods: Magnetic resonance images (MRI) and gait data of one typically developing child and three children with cerebral palsy were analysed. Musculoskeletal simulations were performed to calculate joint kinematics, joint kinetics, muscle forces and joint contact forces using four modelling frameworks: 1) Generic-scaled model with static optimization, 2) Generic-scaled model with an electromyography-informed approach, 3) MRI-based model with static optimization, and 4) MRI-based model with an electromyography-informed approach. Findings: Root-mean-square-differences in joint kinematics and kinetics between generic-scaled and MRI-based models were below 5° and 0.15 Nm/kg, respectively. Root-mean-square-differences over all muscles was below 0.2 body weight for every participant. Root-mean-square-differences in joint contact forces between the different modelling frameworks were up to 2.2 body weight. Comparing the simulation results from the typically developing child with the results from the children with cerebral palsy showed similar root-mean-square-differences for all modelling frameworks. Interpretation: In our participants, the impact of MRI-based models on joint contact forces was higher than the impact of including electromyography. Clinical reasoning based on overall root-mean-square-differences in musculoskeletal simulation results between healthy and pathological participants are unlikely to be affected by the modelling choice.

Background:Children with cerebral palsy (CP) present with a pathological gait pattern due to musculoskeletal impairments, such as muscle weakness and altered bony geometry. However, the effect of these impairments on gait performance is still unknown. Research aim:This study aimed to explore the effect of hip muscle weakness and femoral deformities on the gait performance of CP and typical developing (TD) subjects. Methods:6400 musculoskeletal models were created by weakening the hip extensors, abductors, adductors and flexors from 0% to 75 % and increasing the femoral anteversion angle (FAA) and neck shaft angle (NSA) from 20° to 60° and 120° to 160°, respectively. One TD and five CP gait patterns were imposed to each model and muscle forces were calculated. The effect of weakness and bony deformities on the capability gap (CG) at the hip, i.e. the lack in hip moment generating capacity to perform the gait pattern, was investigated using regression analysis. Results:The CG of apparent equinus, stiff knee gait, TD gait, jump gait and true equinus increased with 0.080, 0.038, 0.015, 0.023 and 0.005 Nm/kg per 10 percent hip abductor weakness increase, with 0.211, 0.130, 0.056, 0.045 and 0.011 Nm/kg per 10 degrees FAA increase and with 0.163, 0.080, 0.036, 0.043 and 0.011 Nm/kg per 10 degrees NSA increase, respectively. Combined weakness and bony deformities explained 96 %, 85 %, 82 %, 65 %, 40 % and 35 % of the variance in the CG of apparent equinus, TD gait, stiff knee gait, jump gait, true equinus and crouch gait, respectively. Significance:The results suggest that surgical correction of femoral deformities is more likely to be effective than strength training of hip muscles in enhancing CP gait performance. Jump gait, true equinus and especially crouch were more robust, while apparent equinus and stiff knee gait were limited by hip weakness and femoral deformities.

Multi-scale simulations, combining muscle and joint contact force (JCF) from musculoskeletal simulations with adaptive mechanobiological finite element analysis, allow to estimate musculoskeletal loading and predict femoral growth in children. Generic linearly scaled musculoskeletal models are commonly used. This approach, however, neglects subject- and age-specific musculoskeletal geometry, e.g. femoral neck-shaft angle (NSA) and anteversion angle (AVA). This study aimed to evaluate the impact of proximal femoral geometry, i.e. altered NSA and AVA, on hip JCF and femoral growth simulations. Musculoskeletal models with NSA ranging from 120° to 150° and AVA ranging from 20° to 50° were created and used to calculate muscle and hip JCF based on the gait analysis data of a typically developing child. A finite element model of a paediatric femur was created from magnetic resonance images. The finite element model was morphed to the geometries of the different musculoskeletal models and used for mechanobiological finite element analysis to predict femoral growth trends. Our findings showed that hip JCF increase with increasing NSA and AVA. Furthermore, the orientation of the hip JCF followed the orientation of the femoral neck axis. Consequently, the osteogenic index, which is a function of cartilage stresses and defines the growth rate, barely changed with altered NSA and AVA. Nevertheless, growth predictions were sensitive to the femoral geometry due to changes in the predicted growth directions. Altered NSA had a bigger impact on the growth results than altered AVA. Growth simulations based on mechanobiological principles were in agreement with reported changes in paediatric populations.

Background: Apart from biomechanical alterations in movement patterns, it is known that movement limitations in persons with knee osteoarthritis (PwKOA) are related to an individual's perception and belief regarding pain and disability. To gain more insights into the functional movement behaviour of PwKOA in a clinical setting, inertial sensor technology can be applied. This study first aims to evaluate the ability of inertial sensors to discriminate between healthy controls (HC) and PwKOA. Secondly, this study aims to determine the relationship between movement behaviour, pain-related factors and disability scores. Methods: Twelve HC and 19 PwKOA were included. Five repetitions of six functional movement tasks (walking, forward lunge, sideward lunge, ascent and descent stairs, single leg squat and sit-to-stand) were simultaneously recorded by the inertial sensor system and a camera-based motion analysis system. Statistically significant differences in angular waveforms of the trunk, pelvis and lower limb joints between HC and PwKOA were determined using one-dimensional statistical parametric mapping (SPM1D). The Knee injury and Osteoarthritis Outcome Score and TAMPA scale for Kinesiophobia were used to evaluate the relationship between discriminating joint motion, pain-related factors and disability using spearman's correlation coefficients. Results: PwKOA had significantly less trunk rotation, internal pelvis rotation and knee flexion ROM during walking. Additionally, the reduced knee flexion (i.e. at the end of the stance phase and swing phase) was related to increased level of perceived pain. During the sideward lunge, PwKOA had significantly less knee flexion, ankle plantarflexion and hip abduction. This decreased hip abduction (i.e. during stance) was related to higher fear of movement. Finally, PwKOA had significantly less knee flexion during the forward lunge, single leg squat and during ascent and descent stairs. No significant correlations were observed with disability. Conclusions: Inertial sensors were able to discriminate between movement characteristics of PwKOA and HC. Additionally, significant relationships were found between joint motion, perceived pain and fear of movement. Since inertial sensors can be used outside the laboratory setting, these results are promising as they indicate the ability to evaluate movement deviations. Further research is required to enable measurements of small movement deviations in clinically relevant tasks.

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.

Background:The purpose of this study was to compare muscle forces and hip contact forces (HCFs) during squatting in patients with cam-type femoroacetabular impingement (cam-FAI) before and after hip corrective surgery and with healthy control participants.Methods:Ten symptomatic male patients with cam-FAI performed deep squatting preoperatively and at 2 years postoperatively. Patients were matched by age and body mass index to 10 male control participants. Full-body kinematics and kinetics were computed, and muscle forces and HCFs were estimated using a musculoskeletal model and static optimization. Normalized squat cycle (% _SC) trials were compared using statistical nonparametric mapping (SnPM).Results:Postoperatively, patients with cam-FAI squatted down with higher anterior pelvic tilt, higher hip flexion, and greater hip extension moments than preoperatively. Preoperative patients demonstrated lower anterior pelvic tilt and lower hip flexion compared with the participants in the control group. Postoperative patients showed increased semimembranosus force compared with their preoperative values. Preoperative forces were lower than the control group for the adductor magnus, the psoas major, and the semimembranosus; however, the preoperative patients showed greater inferior gluteus maximus forces than the patients in the control group, whereas the postoperative patients did not differ from the control patients. Higher posterior, superior, and resultant HCF magnitudes were identified postoperatively in comparison with the preoperative values. Preoperative posterior HCF was lower than in the control group, whereas the postoperative posterior HCF did not differ from those in the control group.Conclusions:Higher postoperative anterior pelvic tilt was associated with an indication of return to closer to normal pelvic motion, which resembled data from the control group. Lower preoperative anterior pelvic tilt was associated with muscle force imbalance, indicated by decreased semimembranosus and increased gluteus maximus forces. The overall increased postoperative muscle forces were associated with improved pelvic mobility and increased HCFs that were comparable with the control-group standards.Clinical Relevance:Muscle forces and HCFs may be indicative of postoperative joint health restoration and alleviated symptoms.

Background The unipodal stance task is a clinical task that quantifies postural stability and alignment of the lower limb joints, while weight bearing on one leg. As persons with knee osteoarthritis (KOA) have poor postural and knee joint stability, objective assessment of this task might be useful. Objective To investigate the discriminant validity of three-dimensional joint kinematics and centre of mass displacement (COM) between healthy controls and persons with knee KOA, during unipodal stance using inertial sensors. Additionally, the reliability, agreement and construct validity are assessed to determine the reproducibility and accuracy of the discriminating parameters. Methods Twenty healthy controls and 19 persons with unilateral severe KOA were included. Five repetitions of the unipodal stance task were simultaneously recorded by an inertial sensor system and a camera-based system (gold standard). Statistical significant differences in kinematic waveforms between healthy controls and persons with severe knee KOA were determined using one-dimensional statistical parametric mapping (SPM1D). Results Persons with severe knee KOA had more lateral trunk lean towards the contralateral leg, more hip flexion throughout the performance of the unipodal stance task, more pelvic obliquity and COM displacement towards the contralateral side. However, for the latter two parameters the minimum detectable change was greater than the difference between healthy controls and persons with severe knee KOA. The construct validity was good (coefficient of multiple correlation 0.75, 0.83 respectively) and the root mean squared error (RMSE) was low (RMSE <1.5°) for the discriminant parameters. Conclusion Inertial sensor based movement analysis can discriminate between healthy controls and persons with severe knee KOA for lateral trunk lean and hip flexion, but unfortunately not for the knee angles. Further research is required to improve the reproducibility and accuracy of the inertial sensor measurements before they can be used to assess differences in tasks with a small range of motion.

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.

With aging, hand mobility and manual dexterity decline, even under healthy circumstances. To assess how aging affects finger movement control, we compared elderly and young subjects with respect to (1) finger movement independence, (2) neural control of extrinsic finger muscles and (3) finger tendon displacements during single finger flexion. In twelve healthy older (age 68–84) and nine young (age 22–29) subjects, finger kinematics were measured to assess finger movement enslaving and the range of independent finger movement. Muscle activation was assessed using a multi-channel electrode grid placed over the flexor digitorum superficialis (FDS) and the extensor digitorum (ED). FDS tendon displacements of the index, middle and ring fingers were measured using ultrasound. In older subjects compared to the younger subjects, we found: (1) increased enslaving of the middle finger during index finger flexion (young: 25.6 ± 12.4%, elderly: 47.0 ± 25.1%; p = 0.018), (2) a lower range of independent movement of the index finger (young _middle = 74.0%, elderly _middle : 45.9%; p < 0.001), (3) a more evenly distributed muscle activation pattern over the finger-specific FDS and ED muscle regions and (4) a lower slope at the beginning of the finger movement to tendon displacement relationship, presenting a distinct period with little to no tendon displacement. Our study indicates that primarily the movement independence of the index finger is affected by aging. This can partly be attributed to a muscle activation pattern that is more evenly distributed over the finger-specific FDS and ED muscle regions in the elderly.

To improve our understanding on the neuromechanics of finger movements, a comprehensive musculoskeletal model is needed. The aim of this study was to build a musculoskeletal model of the hand and wrist, based on one consistent data set of the relevant anatomical parameters. We built and tested a model including the hand and wrist segments, as well as the muscles of the forearm and hand in OpenSim. In total, the model comprises 19 segments (with the carpal bones modeled as one segment) with 23 degrees of freedom and 43 muscles. All required anatomical input data, including bone masses and inertias, joint axis positions and orientations as well as muscle morphological parameters (i.e. PCSA, mass, optimal fiber length and tendon length) were obtained from one cadaver of which the data set was recently published. Model validity was investigated by first comparing computed muscle moment arms at the index finger metacarpophalangeal (MCP) joint and wrist joint to published reference values. Secondly, the muscle forces during pinching were computed using static optimization and compared to previously measured intraoperative reference values. Computed and measured moment arms of muscles at both index MCP and wrist showed high correlation coefficients (r = 0.88 averaged across all muscles) and modest root mean square deviation (RMSD = 23% averaged across all muscles). Computed extrinsic flexor forces of the index finger during index pinch task were within one standard deviation of previously measured in-vivo tendon forces. These results provide an indication of model validity for use in estimating muscle forces during static tasks.