Hip contact forces can be predicted with a neural network using only synthesised key points and electromyography in people with hip osteoarthritis

OBJECTIVE

To develop and validate a neural network to estimate hip contact forces (HCF), and lower body kinematics and kinetics during walking in individuals with hip osteoarthritis (OA) using synthesised anatomical key points and electromyography. To assess the capability of the neural network to detect directional changes in HCF resulting from prescribed gait modifications.

DESIGN

A calibrated electromyography-informed neuromusculoskeletal model was used to compute lower body joint angles, moments, and HCF for 17 participants with mild-to-moderate hip OA. Anatomical key points (e.g., joint centres) were synthesised from marker trajectories and augmented with bias and noise expected from computer vision-based pose estimation systems. Temporal convolutional and long short-term memory neural networks (NN) were trained using leave-one-subject-out validation to predict neuromusculoskeletal modelling outputs from the synthesised key points and measured electromyography data from 5 hip-spanning muscles.

RESULTS

HCF was predicted with an average error of 13.4 ± 7.1% of peak force. Joint angles and moments were predicted with an average root-mean-square-error of 5.3 degrees and 0.10 Nm/kg, respectively. The NN could detect changes in peak HCF that occur due to gait modifications with good agreement with neuromusculoskeletal modelling (r2 = 0.72) and a minimum detectable change of 9.5%.

CONCLUSION

The developed neural network predicted HCF and lower body joint angles and moments in individuals with hip OA using noisy synthesised key point locations with acceptable errors. Changes in HCF magnitude due to gait modifications were predicted with high accuracy. These findings have important implications for implementation of load-modification based gait retraining interventions for people with hip OA in a natural environment (i.e., home, clinic).

Effect of different constraining boundary conditions on simulated femoral stresses and strains during gait

Finite element analysis (FEA) is commonly used in orthopaedic research to estimate localised tissue stresses and strains. A variety of boundary conditions have been proposed for isolated femur analysis, but it remains unclear how these assumed constraints influence FEA predictions of bone biomechanics. This study compared the femoral head deflection (FHD), stresses, and strains elicited under four commonly used boundary conditions (fixed knee, mid-shaft constraint, springs, and isostatic methods) and benchmarked these mechanics against the gold standard inertia relief method for normal and pathological femurs (extreme anteversion and retroversion, coxa vara, and coxa valga). Simulations were performed for the stance phase of walking with the applied femoral loading determined from patient-specific neuromusculoskeletal models. Due to unrealistic biomechanics observed for the commonly used boundary conditions, we propose a novel biomechanical constraint method to generate physiological femur biomechanics. The biomechanical method yielded FHD (< 1 mm), strains (approaching 1000 µε), and stresses (< 60 MPa), which were consistent with physiological observations and similar to predictions from the inertia relief method (average coefficient of determination = 0.97, average normalized root mean square error = 0.17). Our results highlight the superior performance of the biomechanical method compared to current methods of constraint for  both healthy and pathological femurs.

Muscle Morphology Does Not Solely Determine Knee Flexion Weakness After Anterior Cruciate Ligament Reconstruction with a Semitendinosus Tendon Graft: A Combined Experimental and Computational Modeling Study

The distal semitendinosus tendon is commonly harvested for anterior cruciate ligament reconstruction, inducing substantial morbidity at the knee. The aim of this study was to probe how morphological changes of the semitendinosus muscle after harvest of its distal tendon for anterior cruciate ligament reconstruction affects knee flexion strength and whether the knee flexor synergists can compensate for the knee flexion weakness. Ten participants 8-18 months after anterior cruciate ligament reconstruction with an ipsilateral distal semitendinosus tendon autograft performed isometric knee flexion strength testing (15°, 45°, 60°, and 90°; 0° = knee extension) positioned prone on an isokinetic dynamometer. Morphological parameters extracted from magnetic resonance images were used to inform a musculoskeletal model. Knee flexion moments estimated by the model were then compared with those measured experimentally at each knee angle position. A statistically significant between-leg difference in experimentally-measured maximal isometric strength was found at 60° and 90°, but not 15° or 45°, of knee flexion. The musculoskeletal model matched the between-leg differences observed in experimental knee flexion moments at 15° and 45° but did not well estimate between-leg differences with a more flexed knee, particularly at 90°. Further, the knee flexor synergists could not physiologically compensate for weakness in deep knee flexion. These results suggest additional factors other than knee flexor muscle morphology play a role in knee flexion weakness following anterior cruciate ligament reconstruction with a distal semitendinosus tendon graft and thus more work at neural and microscopic levels is required for informing treatment and rehabilitation in this demographic.

Joint contact forces during semi-recumbent seated cycling

Semi-recumbent cycling performed from a wheelchair is a popular rehabilitation exercise following spinal cord injury (SCI) and is often paired with functional electrical stimulation. However, biomechanical assessment of this cycling modality is lacking, even in unimpaired populations, hindering the development of personalised and safe rehabilitation programs for those with SCI. This study developed a computational pipeline to determine lower limb kinematics, kinetics, and joint contact forces (JCF) in 11 unimpaired participants during voluntary semi-recumbent cycling using a rehabilitation ergometer. Two cadences (40 and 60 revolutions per minute) and three crank powers (15 W, 30 W, and 45 W) were assessed. A rigid body model of a rehabilitation ergometer was combined with a calibrated electromyogram-informed neuromusculoskeletal model to determine JCF at the hip, knee, and ankle. Joint excursions remained consistent across all cadence and powers, but joint moments and JCF differed between 40 and 60 revolutions per minute, with peak JCF force significantly greater at 40 compared to 60 revolutions per minute for all crank powers. Poor correlations were found between mean crank power and peak JCF across all joints. This study provides foundation data and computational methods to enable further evaluation and optimisation of semi-recumbent cycling for application in rehabilitation after SCI and other neurological disorders.

Hip Contact Forces During Sprinting in Femoroacetabular Impingement Syndrome

PURPOSE

Sprinting often provokes hip pain in individuals with femoroacetabular impingement syndrome (FAIS). Asphericity of the femoral head-neck junction (cam morphology) characteristic of FAIS can increase the risk of anterior-superior acetabular cartilage damage. This study aimed to 1) compare hip contact forces (magnitude and direction) during sprinting between individuals with FAIS, asymptomatic cam morphology (CAM), and controls without cam morphology, and 2) identify the phases of sprinting with high levels of anteriorly directed hip contact forces.

METHODS

Forty-six recreationally active individuals with comparable levels of physical activity were divided into three groups (FAIS, 14; CAM, 15; control, 17) based on their history of hip/groin pain, results of clinical impingement tests, and presence of cam morphology (alpha angle >55°). Three-dimensional marker trajectories, ground reaction forces, and electromyograms from 12 lower-limb muscles were recorded during 10-m overground sprinting trials. A linearly scaled electromyogram-informed neuromusculoskeletal model was used to calculate hip contact force magnitude (resultant, anterior-posterior, inferior-superior, medio-lateral) and angle (sagittal and frontal planes). Between-group comparisons were made using two-sample t -tests via statistical parametric mapping ( P < 0.05).

RESULTS

No significant differences in magnitude or direction of hip contact forces were observed between FAIS and CAM or between FAIS and control groups during any phase of the sprint cycle. The highest anteriorly directed hip contact forces were observed during the initial swing phase of the sprint cycle.

CONCLUSIONS

Hip contact forces during sprinting do not differentiate recreationally active individuals with FAIS from asymptomatic individuals with and without cam morphology. Hip loading during early swing, where peak anterior loading occurs, may be a potential mechanism for cartilage damage during sprinting-related sports in individuals with FAIS and/or asymptomatic cam morphology.

Moment arm and torque generating capacity of semitendinosus following tendon harvesting for anterior cruciate ligament reconstruction: A simulation study

Altered semitendinosus (ST) morphology and distal tendon insertion following anterior cruciate ligament reconstruction (ACLR) may reduce knee flexion torque generating capacity of the hamstrings via impaired ST force generation and/or moment arm. This study used a computational musculoskeletal model to simulate mechanical consequences of tendon harvest for ACLR on ST function by modeling changes in ST muscle tendon insertion point, moment arm, and torque generating capacity across a physiological range of motion. Simulated ST function was then compared between ACLR and uninjured contralateral limbs. Magnetic resonance imaging from 18 individuals with unilateral history of ACLR involving a hamstring autograft was used to analyse bilateral hamstring muscle (ST, semimembranosus, bicep femoris long head and short head) morphology and distal ST tendon insertion. The ACLR cohort was sub-grouped into those with and without ST regeneration. For each participant with ST regeneration (n = 7), a personalized musculoskeletal model was created including postoperative remodeling of ST using OpenSim 4.1. Knee flexion and internal rotation moment arms and torque generating capacities of hamstrings were evaluated. Bilateral differences were calculated with an asymmetry index (%) ([unaffected limb-affected limb]/[unaffected limb + affected limb]*100%). Smaller moment arms or knee torques within injured compared to uninjured contralateral limbs were considered a deficit. Compared to uninjured contralateral limbs, ACLR limbs with tendon regeneration (n = 7) had minor reductions in knee flexion (5.80% [95% confidence interval (CI) = 3.97-7.62]) and internal rotation (4.92% [95% CI = 2.77-7.07]) moment arms. Decoupled from muscle morphology, altered ST moment arms in ACLR limbs with tendon regeneration resulted in negligible deficits in knee flexion (1.20% [95% CI = 0.34-2.06]) and internal rotation (0.24% [95% CI = 0.22-0.26]) torque generating capacity compared to uninjured contralateral limbs. Coupled with muscle morphology, ACLR limbs with tendon regeneration had substantial deficits in knee flexion (19.32% [95% CI = 18.35-20.28]) and internal rotation (15.49% [95% CI = 14.56-16.41]) torques compared to uninjured contralateral limbs. Personalized musculoskeletal models with measures of ST distal insertion and muscle morphology provided unique insights into post-ACLR ST and hamstring function. Deficits in knee flexor and internal rotation moment arms and torque generating capacities were evident in those with ACLR even when tendon regeneration occurred. Future studies may wish to implement this framework in personalized musculoskeletal models following ACLR to better understand individual muscle function for injury prevention and treatment evaluation.

Hamstring harvest results in significantly reduced knee muscular protection during side-step cutting two years after anterior cruciate ligament reconstruction

The purpose of this study was to determine the effect of donor muscle morphology following tendon harvest in anterior cruciate ligament (ACL) reconstruction on muscular support of the tibiofemoral joint during sidestep cutting. Magnetic resonance imaging (MRI) was used to measure peak cross-sectional area (CSA) and volume of the semitendinosus (ST) and gracilis (GR) muscles and tendons (bilaterally) in 18 individuals following ACL reconstruction. Participants performed sidestep cutting tasks in a biomechanics laboratory during which lower-limb electromyography, ground reaction loads, whole-body motions were recorded. An EMG driven neuro-musculoskeletal model was subsequently used to determine force from 34 musculotendinous units of the lower limb and the contribution of the ST and GR to muscular support of the tibiofemoral joint based on a normal muscle-tendon model (Standard model). Then, differences in peak CSA and volume between the ipsilateral/contralateral ST and GR were used to adjust their muscle-tendon parameters in the model followed by a recalibration to determine muscle force for 34 musculotendinous units (Adjusted model). The combined contribution of the donor muscles to muscular support about the medial and lateral compartments were reduced by 52% and 42%, respectively, in the adjusted compared to standard model. While the semimembranosus (SM) increased its contribution to muscular stabilisation about the medial and lateral compartment by 23% and 30%, respectively. This computer simulation study demonstrated the muscles harvested for ACL reconstruction reduced their support of the tibiofemoral joint during sidestep cutting, while the SM may have the potential to partially offset these reductions. This suggests donor muscle impairment could be a factor that contributes to ipsilateral re-injury rates to the ACL following return to sport.

A Digital Twin Framework for Precision Neuromusculoskeletal Health Care: Extension Upon Industrial Standards

There is a powerful global trend toward deeper integration of digital twins into modern life driven by Industry 4.0 and 5.0. Defense, agriculture, engineering, manufacturing, and urban planning sectors have thoroughly incorporated digital twins to great benefit across their respective product lifecycles. Despite clear benefits, a digital twin framework for health and medical sectors is yet to emerge. This paper proposes a digital twin framework for precision neuromusculoskeletal health care. We build upon the International Standards Organization framework for digital twins for manufacturing by presenting best available computational models within a digital twin framework for clinical application. We map a use case for modeling Achilles tendon mechanobiology, highlighting how current modeling practices align with our proposed digital twin framework. Similarly, we map a use case for advanced neurorehabilitation technology, highlighting the role of a digital twin in control of systems where human and machine are interfaced. Future work must now focus on creating an informatic representation to govern how digital data are passed to, from, and within the digital twin, as well as specific standards to declare which measurement systems and modeling methods are acceptable to move toward widespread use of the digital twin framework for precision neuromusculoskeletal health care.

Running Mechanics After Repeated Sprints in Femoroacetabular Impingement Syndrome, Cam Morphology, and Controls

BACKGROUND

People with femoroacetabular with femoroacetabular impingement syndrome (FAIS) often report pain during sports involving repeated sprinting. It remains unclear how sports participation influences running biomechanics in individuals with FAIS.

HYPOTHESIS

Changes in running biomechanics and/or isometric hip strength after repeated sprint exercise would be greatest in individuals with FAIS compared with asymptomatic individuals with (CAM) and without cam morphology (Control).

STUDY DESIGN

Controlled laboratory study.

LEVEL OF EVIDENCE

Level 3.

METHODS

Three-dimensional hip biomechanics during maximal running (10 m) and hip strength were measured in 49 recreationally active individuals (FAIS = 15; CAM = 16; Control = 18) before and after repeated sprint exercise performed on a nonmotorized treadmill (8-16 × 30 m). Effects of group and time were assessed for biomechanics and strength variables with repeated-measures analyses of variance. Relationships between hip pain (Copenhagen Hip and Groin Outcome Score) and changes in hip moments and strength after repeated sprint exercise were determined using Spearman's correlation coefficients (ρ).

RESULTS

Running speed, hip flexion angles, hip flexion and extension moments, and hip strength in all muscle groups were significantly reduced from pre to post. No significant between-group differences were observed before or after repeated sprint exercise. No significant relationships (ρ = 0.04-0.30) were observed between hip pain and changes in hip moments or strength in the FAIS group.

CONCLUSION

Changes in running biomechanics and strength after repeated sprint exercise did not differ between participants with FAIS and asymptomatic participants with and without cam morphology. Self-reported pain did not appear to influence biomechanics during running or strength after repeated sprint exercise in participants with FAIS.

CLINICAL RELEVANCE

A short bout of repeated sprinting may not elicit changes in running biomechanics in FAIS beyond what occurs in those without symptoms. Longer duration activities or activities requiring greater hip flexion angles may better provoke pathology-related changes in running biomechanics in people with FAIS.

Individual muscle contributions to the acceleration of the centre of mass during gait in people with mild-to-moderate hip osteoarthritis

BACKGROUND

People with mild-to-moderate hip osteoarthritis (OA) exhibit hip muscle weakness, alterations in hip kinematics and kinetics and hip contact forces during gait compared to healthy controls. However, it is unclear if those with hip OA use different motor control strategies to coordinate the motion of the centre of mass (COM) during gait. Such information could provide further critical assessment of conservative management strategies implemented for people with hip OA.

RESEARCH QUESTION

Do muscle contributions to the acceleration of the COM during walking differ between individuals with mild-to-moderate hip OA and controls?

METHODS

Eleven individuals with mild-to-moderate hip OA and 10 healthy controls walked at a self-selected speed while whole-body motion and ground reaction forces were measured. Muscle forces during gait were obtained using static optimisation and an induced acceleration analysis was performed to determine individual muscle contributions to the acceleration of the COM during single-leg stance (SLS). Between-group comparisons were made using independent t-tests via Statistical Parametric Modelling.

RESULTS

There were no between-group differences in spatial-temporal gait parameters or three-dimensional whole-body COM acceleration. The rectus femoris, biceps femoris, iliopsoas and gastrocnemius muscles in the hip OA group contributed less to the fore-aft accelerations of the COM (p < 0.05), and more to the vertical COM acceleration with the gluteus maximus (p < 0.05), during SLS, compared to the control group.

SIGNIFICANCE

Subtle differences exist in the way people with mild-to-moderate hip OA use their muscles to accelerate the whole-body centre of mass during the SLS phase of walking relative to healthy controls. These findings improve understanding of the complex functional consequences of hip OA and enhance our understanding of how to monitor the effectiveness of an intervention on biomechanical changes in gait in people with hip OA.

The effect of surgical change to hip geometry on hip biomechanics after primary total hip arthroplasty

The aim of this study was to determine the effect of surgical change to the acetabular offset and femoral offset on the abductor muscle and hip contact forces after primary total hip arthroplasty (THA) using computational methods. Thirty-five patients undergoing primary THA were recruited. Patients underwent a computed tomography scan of their pelvis and hip, and underwent gait analysis pre- and 6-months postoperatively. Surgically induced changes in acetabular and femoral offset were used to inform a musculoskeletal model to estimated abductor muscle and hip joint contact forces. Two experiments were performed: (1) influence of changes in hip geometry on hip biomechanics with preoperative kinematics; and (2) influence of changes in hip geometry on hip biomechanics with postoperative kinematics. Superior and medial placement of the hip centre of rotation during THA was most influential in reducing hip contact forces, predicting 63% of the variance (p < 0.001). When comparing the preoperative geometry and kinematics model, with postoperative geometry and kinematics, hip contact forces increased after surgery (0.68 BW, p = 0.001). Increasing the abductor lever arm reduced abductor muscle force by 28% (p < 0.001) and resultant hip contact force by 17% (0.6 BW, p = 0.003), with both preoperative and postoperative kinematics. Failure to increase abductor lever arm increased resultant hip contact force 11% (0.33 BW, p < 0.001). In conclusion, increasing the abductor lever arm provides a substantial biomechanical benefit to reduce hip abductor and resultant hip joint contact forces. The magnitude of this effect is equivalent to the average increase in hip contact force seen with improved gait from pre-to post-surgery.

Maintaining soldier musculoskeletal health using personalised digital humans, wearables and/or computer vision

OBJECTIVES

The physical demands of military service place soldiers at risk of musculoskeletal injuries and are major concerns for military capability. This paper outlines the development new training technologies to prevent and manage these injuries.

DESIGN

Narrative review.

METHODS

Technologies suitable for integration into next-generation training devices were examined. We considered the capability of technologies to target tissue level mechanics, provide appropriate real-time feedback, and their useability in-the-field.

RESULTS

Musculoskeletal tissues' health depends on their functional mechanical environment experienced in military activities, training and rehabilitation. These environments result from the interactions between tissue motion, loading, biology, and morphology. Maintaining health of and/or repairing joint tissues requires targeting the "ideal" in vivo tissue mechanics (i.e., loading and strain), which may be enabled by real-time biofeedback. Recent research has shown that these biofeedback technologies are possible by integrating a patient's personalised digital twin and wireless wearable devices. Personalised digital twins are personalised neuromusculoskeletal rigid body and finite element models that work in real-time by code optimisation and artificial intelligence. Model personalisation is crucial in obtaining physically and physiologically valid predictions.

CONCLUSIONS

Recent work has shown that laboratory-quality biomechanical measurements and modelling can be performed outside the laboratory with a small number of wearable sensors or computer vision methods. The next stage is to combine these technologies into well-designed easy to use products.

Morphology of proximal and distal human semitendinosus compartments and the effects of distal tendon harvesting for anterior cruciate ligament reconstruction

The human semitendinosus muscle is characterized by a tendinous inscription separating proximal and distal neuromuscular compartments. As each compartment is innervated by separate nerve branches, potential exists for independent operation and control of compartments. However, the morphology and function of each compartment have not been thoroughly examined in an adult human population. Further, the distal semitendinosus tendon is typically harvested for use in anterior cruciate ligament reconstruction surgery, which induces long-term morphological changes to the semitendinosus muscle-tendon unit. It remains unknown if muscle morphological alterations following anterior cruciate ligament reconstruction are uniform between proximal and distal semitendinosus compartments. Here, we performed magnetic resonance imaging on 10 individuals who had undergone anterior cruciate ligament reconstruction involving an ipsilateral distal semitendinosus tendon graft 14 ± 4 months prior, extracting morphological parameters of the whole semitendinosus muscle and each individual compartment from both the (non-injured) contralateral and surgical legs. In the contralateral leg, volume and length of the proximal compartment were smaller than the distal compartment. No between-compartment differences in volume or length were found for anterior cruciate ligament reconstructed legs, likely due to greater shortening of the distal compared to the proximal compartment after anterior cruciate ligament reconstruction. The maximal anatomical cross-sectional area of both compartments was substantially smaller on the anterior cruciate ligament reconstructed leg but did not differ between compartments on either leg. The absolute and relative between-leg differences in proximal compartment morphology on the anterior cruciate ligament reconstructed leg were strongly correlated with the corresponding between-leg differences in distal compartment morphological parameters. Specifically, greater between-leg morphological differences in one compartment were highly correlated with large between-leg differences in the other compartment, and vice versa for smaller differences. These relationships indicate that despite the heterogeneity in compartment length and volume, compartment atrophy is not independent or random. Further, the tendinous inscription endpoints were generally positioned at the same proximodistal level as the compartment maximal anatomical cross-sectional areas, providing a wide area over which the tendinous inscription could mechanically interact with compartments. Overall, results suggest the two human semitendinosus compartments are not mechanically independent.

Sex-Specific Physical Performance Adaptive Responses Are Elicited After 10 Weeks of Load Carriage Conditioning

INTRODUCTION

The purpose of this study was to identify and characterize sex-specific physical and psychophysical performance adaptations in response to a novel 10-week training program.

MATERIALS AND METHODS

Fifteen males and thirteen females completed a standardized load carriage task (5 km at 5.5 km.h-1, wearing a 23 kg torso-borne vest) before and after 10 weeks of resistance and load carriage training. Psychophysical responses (i.e., heart rate and ratings of perceived exertion) were measured throughout the load carriage task. Physical performance (i.e., countermovement and squat jumps, push-ups, sit-ups, and beep test) was measured at before, mid-way, and after the training program (weeks 0, 6, and 11, respectively).

RESULTS

Training elicited significant improvements in squat jump maximal force, push-ups, and beep test performance (P < .05). Males outperformed females in all performance measures, with interactions (time, sex) for push-ups, sit-ups, and beep test performance. After training, aerobic capacity improved by 5.4% (42.9 mL· kg-1· min-1 to 45.2 mL· kg-1· min-1) in males but did not improve in females. Psychophysical responses decreased for both sexes (P < .05) during the load carriage task post-training.

CONCLUSION

While 10 weeks of standardized training elicited positive adaptations in both physical and psychophysical performance, sex-specific differences were still evident. To lessen these differences, sex-specific training should be considered to optimize load carriage performance.

Limiting the Use of Electromyography and Ground Reaction Force Data Changes the Magnitude and Ranking of Modelled Anterior Cruciate Ligament Forces

Neuromusculoskeletal models often require three-dimensional (3D) body motions, ground reaction forces (GRF), and electromyography (EMG) as input data. Acquiring these data in real-world settings is challenging, with barriers such as the cost of instruments, setup time, and operator skills to correctly acquire and interpret data. This study investigated the consequences of limiting EMG and GRF data on modelled anterior cruciate ligament (ACL) forces during a drop–land–jump task in late-/post-pubertal females. We compared ACL forces generated by a reference model (i.e., EMG-informed neural mode combined with 3D GRF) to those generated by an EMG-informed with only vertical GRF, static optimisation with 3D GRF, and static optimisation with only vertical GRF. Results indicated ACL force magnitude during landing (when ACL injury typically occurs) was significantly overestimated if only vertical GRF were used for either EMG-informed or static optimisation neural modes. If 3D GRF were used in combination with static optimisation, ACL force was marginally overestimated compared to the reference model. None of the alternative models maintained rank order of ACL loading magnitudes generated by the reference model. Finally, we observed substantial variability across the study sample in response to limiting EMG and GRF data, indicating need for methods incorporating subject-specific measures of muscle activation patterns and external loading when modelling ACL loading during dynamic motor tasks.

A digital twin framework for robust control of robotic-biological systems

Medical device regulatory standards are increasingly incorporating computational modelling and simulation to accommodate advanced manufacturing and device personalization. We present a method for robust testing of engineered soft tissue products involving a digital twin paradigm in combination with robotic systems. We developed and validated a digital twin framework for calibrating and controlling robotic-biological systems. A forward dynamics model of the robotic manipulator was developed, calibrated, and validated. After calibration, the accuracy of the digital twin in reproducing the experimental data improved in the time domain for all fourteen tested configurations and improved in frequency domain for nine configurations. We then demonstrated displacement control of a spring in lieu of a soft tissue element in a biological specimen. The simulated experiment matched the physical experiment with 0.09 mm (0.001%) root-mean-square error for a 2.9 mm (5.1%) length change. Finally, we demonstrated kinematic control of a digital twin of the knee through 70-degree passive flexion kinematics. The root-mean-square error was 2.00°, 0.57°, and 1.75° degrees for flexion, adduction, and internal rotations, respectively. The system well controlled novel mechanical elements and generated accurate kinematics in silico for a complex knee model. This calibration method could be applied to other situations where the specimen is poorly represented in the model environment (e.g., human or animal tissues), and the control system could be extended to track internal parameters such as tissue strain (e.g., control knee ligament strain). Further development of this framework can facilitate medical device testing and innovative biomechanics research.

Additively manufactured polyethylene terephthalate scaffolds for scapholunate interosseous ligament reconstruction

The regeneration of the ruptured scapholunate interosseous ligament (SLIL) represents a clinical challenge. Here, we propose the use of a Bone-Ligament-Bone (BLB) 3D-printed polyethylene terephthalate (PET) scaffold for achieving mechanical stabilisation of the scaphoid and lunate following SLIL rupture. The BLB scaffold featured two bone compartments bridged by aligned fibres (ligament compartment) mimicking the architecture of the native tissue. The scaffold presented tensile stiffness in the range of 260 ± 38 N/mm and ultimate load of 113 ± 13 N, which would support physiological loading. A finite element analysis (FEA), using inverse finite element analysis (iFEA) for material property identification, showed an adequate fit between simulation and experimental data. The scaffold was then biofunctionalized using two different methods: injected with a Gelatin Methacryloyl solution containing human mesenchymal stem cell spheroids (hMSC) or seeded with tendon-derived stem cells (TDSC) and placed in a bioreactor to undergo cyclic deformation. The first approach demonstrated high cell viability, as cells migrated out of the spheroid and colonised the interstitial space of the scaffold. These cells adopted an elongated morphology suggesting the internal architecture of the scaffold exerted topographical guidance. The second method demonstrated the high resilience of the scaffold to cyclic deformation and the secretion of a fibroblastic related protein was enhanced by the mechanical stimulation. This process promoted the expression of relevant proteins, such as Tenomodulin (TNMD), indicating mechanical stimulation may enhance cell differentiation and be useful prior to surgical implantation. In conclusion, the PET scaffold presented several promising characteristics for the immediate mechanical stabilisation of disassociated scaphoid and lunate and, in the longer-term, the regeneration of the ruptured SLIL.

Neuromusculoskeletal model calibration accounts for differences in electromechanical delay and maximum isometric muscle force

Electromechanical delay (EMD) and maximum isometric muscle force (F_o^M) are important parameters for joint contact force calculation with EMG-informed neuromusculoskeletal (NMS) models. These parameters can vary between tasks (EMD) and individuals (EMD and F_o^M), making it challenging to establish representative values. One promising approach is to personalise candidate parameters to the participant (e.g., F_o^M by regression equation) and then adjust all parameters within a calibration (i.e., numerical optimisation) to minimise error between corresponding pairs of experimental measures and model-predicted values. The purpose of this study was to determine whether calibration of an NMS model resulted in consistent joint contact forces, regardless of EMD value or personalisation of F_o^M. Hip, knee, and ankle contact forces were predicted for 28 participants using EMG-informed NMS models. Differences in joint contact forces with EMD were examined in six models, calibrated with EMD from 15 to 110 ms. Differences in joint contact forces with personalisation of F_o^M were examined in two models, both calibrated with the same initial EMD (50 ms), one with generic and one with personalised values for F_o^M. For all models, joint contact force peaks during the first and second halves of stance were extracted and compared using a repeated-measures analysis of variance. Calibrated models with EMD set between 35 and 70 ms produced similar magnitude and timing of peak joint contact forces. Compared with generic values, personalising and then calibrating F_o^M resulted in comparable peak contact forces at hip, but not knee or ankle, while also producing muscle-specific tensions similar to reported literature. Overall, EMD between 35 and 70 ms and personalised initial values of F_o^M before calibration are advised for EMG-informed NMS modelling.

Shape differences in the semitendinosus following tendon harvesting for anterior cruciate ligament reconstruction

Following hamstring autograft anterior cruciate ligament reconstruction (ACLR), muscle length, cross-sectional area, and volume are reduced. However, these discrete measures of morphology do not account for complex three-dimensional muscle shape. The primary aim of this study was to determine between-limb semitendinosus (ST) shape and regional morphology differences in young adults following tendon harvest for ACLR and to compare these differences with those in healthy controls. In this cross-sectional study, magnetic resonance imaging was performed on 18 individuals with unilateral ACLR and 18 healthy controls. Bilaterally, ST muscles were segmented, and shape differences assessed between limbs and compared between groups using Jaccard index (0-1) and Hausdorff distance (mm). Length (cm), peak cross-sectional area (cm² ), and volume (cm³ ) were measured for the entire muscle and proximal, middle, and distal regions, and compared between limbs and groups. Compared to healthy controls, the ACLR group had significantly (p < 0.001, Cohen's d = -2.33) lower bilateral ST shape similarity and shape deviation was significantly (p < 0.001, d = 2.12) greater. Shape deviation was greatest within the distal region of the ACLR (Hausdorff: 23.1 ± 8.68 mm). Compared to both the uninjured contralateral limb and healthy controls, deficits in peak cross-sectional area and volume in ACLR group were largest in proximal (p < 0.001, d = -2.52 to -1.28) and middle (p < 0.001, d = -1.81 to -1.04) regions of the ST. Overall, shape analysis provides unique insight into regional adaptations in ST morphology post-ACLR. Findings highlight morphological features in distal ST not identified by traditional discrete morphology measures. Clinical significance: Following ACLR, risk of a secondary knee or primary hamstring injury has been reported to be between 2-to-5 times greater compared to those without ACLR. Change in semitendinosus (ST) shape following ACLR may affect force transmission and distribution within the hamstrings and might contribute to persistent deficits in knee flexor and internal rotator strength.

Semitendinosus muscle morphology in relation to surface electrode placement in anterior cruciate ligament reconstructed and contralateral legs

The semitendinosus tendon is commonly harvested as graft tissue for anterior cruciate ligament reconstruction (ACLR). Although the semitendinosus tendon can regenerate following harvesting, ACLR results in substantial reductions in semitendinosus muscle size and length, potentially complicating electrode placement for electromyography. The purpose of this study was to assess whether the most commonly used electrode placement [recommended by the “Surface Electromyography for Non-Invasive Assessment of Muscles” (SENIAM) project] is appropriate for measuring semitendinosus electromyograms after ACLR. In nine participants (unilateral ACLR with a semitendinosus graft), B-mode ultrasonography was used to bilaterally determine (i) the semitendinosus muscle-tendon junction position and the state of tendon regeneration (latter for the ACLR leg only) and (ii) the anatomical cross-sectional area (ACSA) of the semitendinosus muscle at the SENIAM-recommended electrode placement site at rest and during isometric maximal voluntary contraction (MVC) at two knee joint angles. Depending on the contraction state and joint angle, the semitendinosus muscle had retracted past the recommended placement site in 33–78% of ACLR legs, but not in any contralateral legs. The ACSA of semitendinosus was smaller both at rest and MVC in the ACLR compared to contralateral leg. The ACSA for both legs decreased at MVC compared to rest and at deep compared to shallow knee flexion angles, likely due to sliding of the muscle under the skin. These results suggest SENIAM guidelines are likely unsuitable for recording surface electromyograms from the semitendinosus muscle after tendon harvesting for ACLR as the muscle of interest may not be within the electrode detection volume.

Game-play affects hamstring but not adductor muscle fibre mechanics in elite U20 basketball athletes

Muscle tendon unit fibre mechanics of hamstring and adductor strain injuries are not well studied, with factors such as fatigue promoted as risk factors in the absence of mechanistic evidence. In this study, musculoskeletal modelling was used to estimate fibre mechanics of four hamstring (biceps femoris long head, biceps femoris short head, semimembranosus and semitendinosus) and four adductor (adductor brevis, adductor longus, adductor magnus and gracilis) muscles during an anticipated cut task. The cut task was performed by 10 healthy elite male U20 basketball players both before and immediately after they played in one (of four) competitive basketball game. Biceps femoris long head produced significantly lower (p = 0.032) submaximal force post-game in the latter part of swing (30.7% to 35.0% of stride), though its peak force occurred later (37%) and remained unchanged. Semimembranosus produced significantly lower (p = 0.006) force post-game (32.9% to 44.9% of stride), which encompassed the instance of peak force (39%). Neither fibre velocity nor fibre length of the investigated muscles were significantly affected by game-play. These finding suggest that if fatigue is a factor in hamstring and adductor muscle strain injuries and is brought about by game-play, it is unlikely through the fibre mechanisms investigated in this study.

Development of predictive statistical shape models for paediatric lower limb bones

BACKGROUND AND OBJECTIVE

Accurate representation of bone shape is important for subject-specific musculoskeletal models as it may influence modelling of joint kinematics, kinetics, and muscle dynamics. Statistical shape modelling is a method to estimate bone shape from minimal information, such as anatomical landmarks, and to avoid the time and cost associated with reconstructing bone shapes from comprehensive medical imaging. Statistical shape models (SSM) of lower limb bones have been developed and validated for adult populations but are not applicable to paediatric populations. This study aimed to develop SSM for paediatric lower limb bones and evaluate their reconstruction accuracy using sparse anatomical landmarks.

METHODS

We created three-dimensional models of 56 femurs, 29 pelves, 56 tibias, 56 fibulas, and 56 patellae through segmentation of magnetic resonance images taken from 29 typically developing children (15 females; 13 ± 3.5 years). The SSM for femur, pelvis, tibia, fibula, patella, haunch (i.e., combined femur and pelvis), and shank (i.e., combined tibia and fibula) were generated from manual segmentation of comprehensive magnetic resonance images to describe the shape variance of the cohort. We implemented a leave-one-out cross-validation method wherein SSM were used to reconstruct novel bones (i.e., those not included in SSM generation) using full- (i.e., full segmentation) and sparse- (i.e., anatomical landmarks) input, and then compared these reconstructions against bones segmented from magnetic resonance imaging. Reconstruction performance was evaluated using root mean squared errors (RMSE, mm), Jaccard index (0-1), Dice similarity coefficient (DSC) (0-1), and Hausdorff distance (mm). All results reported in this abstract are mean ± standard deviation.

RESULTS

Femurs, pelves, tibias, fibulas, and patellae reconstructed via SSM using full-input had RMSE between 0.89 ± 0.10 mm (patella) and 1.98 ± 0.38 mm (pelvis), Jaccard indices between 0.77 ± 0.03 (pelvis) and 0.90 ± 0.02 (tibia), DSC between 0.87 ± 0.02 (pelvis) and 0.95 ± 0.01 (tibia), and Hausdorff distances between 2.45 ± 0.57 mm (patella) and 9.01 ± 2.36 mm (pelvis). Reconstruction using sparse-input had RMSE ranging from 1.33 ± 0.61 mm (patella) to 3.60 ± 1.05 mm (pelvis), Jaccard indices ranging from 0.59 ± 0.10 (pelvis) to 0.83 ± 0.03 (tibia), DSC ranging from 0.74 ± 0.08 (pelvis) to 0.90 ± 0.02 (tibia), and Hausdorff distances ranging from 3.21 ± 1.19 mm (patella) to 12.85 ± 3.24 mm (pelvis).

CONCLUSIONS

The SSM of paediatric lower limb bones showed reconstruction accuracy consistent with previously developed SSM and outperformed adult-based SSM when used to reconstruct paediatric bones.

Comparison of Walking Biomechanics After Physical Therapist-Led Care or Hip Arthroscopy for Femoroacetabular Impingement Syndrome: A Secondary Analysis From a Randomized Controlled Trial

BACKGROUND

Femoroacetabular impingement syndrome is characterized by chondrolabral damage and hip pain. The specific biomechanics used by people with femoroacetabular impingement syndrome during daily activities may exacerbate their symptoms. Femoroacetabular impingement syndrome can be treated nonoperatively or surgically; however, differential treatment effects on walking biomechanics have not been examined.

PURPOSE

To compare the 12-month effects of physical therapist-led care or arthroscopy on trunk, pelvis, and hip kinematics as well as hip moments during walking.

STUDY DESIGN

Secondary analysis of multi-centre, pragmatic, two-arm superiority randomized controlled trial subsample; Level of evidence, 1.

METHODS

A subsample of 43 participants from the Australian Full randomised controlled trial of Arthroscopic Surgery for Hip Impingement versus best cONventional (FASHIoN trial) underwent gait analysis and completed the International Hip Outcome Tool (iHOT-33) at both baseline and 12 months after random allocation to physical therapist-led care (personalized hip therapy; n = 22; mean age 35; 41% female) or arthroscopy (n = 21; mean age 36; 48% female). Changes in trunk, pelvis, and hip biomechanics were compared between treatment groups across the gait cycle using statistical parametric mapping. Associations between changes in iHOT-33 and changes in hip kinematics across 3 planes of motion were examined.

RESULTS

As compared with the arthroscopy group, the personalized hip therapy group increased its peak hip adduction moments (mean difference = 0.35 N·m/body weight·height [%] [95% CI, 0.05-0.65]; effect size = 0.72; P = .02). Hip adduction moments in the arthroscopy group were unchanged in response to treatment. No other between-group differences were detected. Improvements in iHOT-33 were not associated with changes in hip kinematics.

CONCLUSION

Peak hip adduction moments were increased in the personalized hip therapy group and unchanged in the arthroscopy group. No biomechanical changes favoring arthroscopy were detected, suggesting that personalized hip therapy elicits greater changes in hip moments during walking at 12-month follow-up. Twelve-month changes in hip-related quality of life were not associated with changes in hip kinematics.

Study protocol for double-blind, randomised placebo-controlled trial evaluating semitendinosus function and morbidity following tendon harvesting for anterior cruciate ligament reconstruction augmented by platelet-rich plasma

INTRODUCTION

Anterior cruciate ligament (ACL) rupture is debilitating, often requiring surgical reconstruction. An ACL reconstruction (ACLR) using a tendon autograft harvested from the semitendinosus results in substantial injury to the donor muscle. Following ACLR, patients rarely return to their preinjury level of physical activity, are at elevated risk of secondary lower limb injuries and early onset knee osteoarthritis. To date, no randomised controlled trial has evaluated the efficacy of platelet-rich plasma (PRP) in aiding knee function and semitendinosus morphology of following ALCR.

METHODS AND ANALYSIS

This is a multicentre double-blind randomised placebo-controlled trial. Fifty-four ACLR patients aged 18-50 years will be randomised to receive either a single application of PRP (ACLR+) or placebo saline (ACLR) into the semitendinosus harvest zone at the time of surgery. All patients will undergo normal postoperative rehabilitation recommended by the attending orthopaedic surgeon or physiotherapist. The primary outcome measure is between-limb difference (ACLR compared with intact contralateral) in isometric knee flexor strength at 60^o knee flexion, collected 10-12 months postsurgery. This primary outcome measure will be statistically compared between groups (ACLR+ and standard ACLR). Secondary outcome measures include bilateral assessments of hamstring muscle morphology via MRI, biomechanical and electromyographic parameters during an anticipated 45° running side-step cut and multidirectional hopping task and patient-reported outcomes questionaries. Additionally, patient-reported outcomes questionaries will be collected before (baseline) as well as immediately after surgery, and at 2-6 weeks, 3-4 months, 10-12 months and 22-24 months postsurgery 10-12 months following surgery.

ETHICS AND DISSEMINATION

Ethics approval has been granted by Griffith University Human Research Ethics Committee, Greenslopes Research and Ethics Committee, and Royal Brisbane & Women's Hospital Human Research Ethics Committee. Results will be submitted for publication in a peer-reviewed medical journal.

TRIAL REGISTRATION NUMBER

ACTRN12618000762257p.

Electromyography measurements of the deep hip muscles do not improve estimates of hip contact force

The deep hip muscles are often omitted in studies investigating hip contact forces using neuromusculoskeletal modelling methods. However, recent evidence indicates the deep hip muscles have potential to change the direction of hip contact force and could have relevance for hip contact loading estimates. Further, it is not known whether deep hip muscle excitation patterns can be accurately estimated using neuromusculoskeletal modelling or require measurement (through invasive and time-consuming methods) to inform models used to estimate hip contact forces. We calculated hip contact forces during walking, squatting, and squat-jumping for 17 participants using electromyography (EMG)-informed neuromusculoskeletal modelling with (informed) and without (synthesized) intramuscular EMG for the deep hip muscles (piriformis, obturator internus, quadratus femoris). Hip contact force magnitude and direction, calculated as the angle between hip contact force and vector from femoral head to acetabular center, were compared between configurations using a paired t-test deployed through statistical parametric mapping (P < 0.05). Additionally, root mean square error, correlation coefficients (R²), and timing accuracy between measured and modelled deep hip muscle excitation patterns were computed. No significant between-configuration differences in hip contact force magnitude or direction were found for any task. However, the synthesized method poorly predicted (R²-range 0.02-0.3) deep hip muscle excitation patterns for all tasks. Consequently, intramuscular EMG of the deep hip muscles may be unnecessary when estimating hip contact force magnitude or direction using EMG-informed neuromusculoskeletal modelling, though is likely essential for investigations of deep hip muscle function.

Effect of a valgus brace on medial tibiofemoral joint contact force in knee osteoarthritis with varus malalignment: A within-participant cross-over randomised study with an uncontrolled observational longitudinal follow-up

Background

Previous investigations on valgus knee bracing have mostly used the external knee adduction moment. This is a critical limitation, as the external knee adduction moment does not account for muscle forces that contribute substantially to the medial tibiofemoral contact force (MTCF) during walking. The aims of this pilot study were to: 1) determine the effect of a valgus knee brace on MTCF; 2) determine whether the effect is more pronounced after 8 weeks of brace use; 3) assess the feasibility of an 8-week brace intervention.

Methods

Participants with medial radiographic knee OA and varus malalignment were fitted with an Össur Unloader One^© brace. Participants were instructed to wear the brace for 8 weeks. The MTCF was estimated via an electromyogram-assisted neuromuscular model with and without the knee brace at week 0 and week 8. Feasibility outcomes included change in symptoms, quality of life, confidence, acceptability, adherence and adverse events.

Results

Of the 30 (60% male) participants enrolled, 28 (93%) completed 8-week outcome assessments. There was a main effect of the brace (p<0.001) on peak MTCF and MTCF impulse, but no main effect for time (week 0 and week 8, p = 0.10), and no interaction between brace and time (p = 0.62). Wearing the brace during walking significantly reduced the peak MTCF (-0.05 BW 95%CI [-0.10, -0.01]) and MTCF impulse (-0.07 BW.s 95%CI [-0.09, -0.05]). Symptoms and quality of life improved by clinically relevant magnitudes over the 8-week intervention. Items relating to confidence and acceptability were rated relatively highly. Participants wore the brace on average 6 hrs per day. Seventeen participants reported 30 minor adverse events over an 8-week period.

Conclusion

Although significant, reductions in the peak MTCF and MTCF while wearing the knee brace were small. No effect of time on MTCF was observed. Although there were numerous minor adverse events, feasibility outcomes were generally favourable.

Trial registration

Australian and New Zealand Clinical Trials Registry (12619000622101).

Quadriceps muscle strength at 2 years following anterior cruciate ligament reconstruction is associated with tibiofemoral joint cartilage volume

PURPOSE

Quadriceps strength deficits following anterior cruciate ligament reconstruction (ACLR) are linked to altered lower extremity biomechanics, tibiofemoral joint (TFJ) space narrowing and cartilage composition changes. It is unknown, however, if quadriceps strength is associated with cartilage volume in the early years following ACLR prior to the onset of posttraumatic osteoarthritis (OA) development. The purpose of this cross-sectional study was to examine the relationship between quadriceps muscle strength (peak and across the functional range of knee flexion) and cartilage volume at ~ 2 years following ACLR and determine the influence of concomitant meniscal pathology.

METHODS

The involved limb of 51 ACLR participants (31 isolated ACLR; 20 combined meniscal pathology) aged 18-40 years were tested at 2.4 ± 0.4 years post-surgery. Isokinetic knee extension torque generated in 10° intervals between 60° and 10° knee flexion (i.e. 60°-50°, 50°-40°, 40°-30°, 30°-20°, 20°-10°) together with peak extension torque were measured. Tibial and patellar cartilage volumes were measured using magnetic resonance imaging (MRI). The relationships between peak and angle-specific knee extension torque and MRI-derived cartilage volumes were evaluated using multiple linear regression.

RESULTS

In ACLR participants with and without meniscal pathology, higher knee extension torques at 60°-50° and 50°-40° knee flexion were negatively associated with medial tibial cartilage volume (p < 0.05). No significant associations were identified between peak concentric or angle-specific knee extension torques and patellar cartilage volume.

CONCLUSION

Higher quadriceps strength at knee flexion angles of 60°-40° was associated with lower cartilage volume on the medial tibia ~ 2 years following ACLR with and without concomitant meniscal injury. Regaining quadriceps strength across important functional ranges of knee flexion after ACLR may reduce the likelihood of developing early TFJ cartilage degenerative changes.

LEVEL OF EVIDENCE

III.

Towards Tailored Rehabilitation by Implementation of a Novel Musculoskeletal Finite Element Analysis Pipeline

Tissue-level mechanics (e.g., stress and strain) are important factors governing tissue remodeling and development of knee osteoarthritis (KOA), and hence, the success of physical rehabilitation. To date, no clinically feasible analysis toolbox has been introduced and used to inform clinical decision making with subject-specific in-depth joint mechanics of different activities. Herein, we utilized a rapid state-of-the-art electromyography-assisted musculoskeletal finite element analysis toolbox with fibril-reinforced poro(visco)elastic cartilages and menisci to investigate knee mechanics in different activities. Tissue mechanical responses, believed to govern collagen damage, cell death, and fixed charge density loss of proteoglycans, were characterized within 15 patients with KOA while various daily activities and rehabilitation exercises were performed. Results showed more inter-participant variation in joint mechanics during rehabilitation exercises compared to daily activities. Accordingly, the devised workflow may be used for designing subject-specific rehabilitation protocols. Further, results showed the potential to tailor rehabilitation exercises, or assess capacity for daily activity modifications, to optimally load knee tissue, especially when mechanically-induced cartilage degeneration and adaptation are of interest.

An EMG-assisted muscle-force driven finite element analysis pipeline to investigate joint- and tissue-level mechanical responses in functional activities: towards a rapid assessment toolbox

Joint tissue mechanics (e.g., stress and strain) are believed to have a major involvement in the onset and progression of musculoskeletal disorders, e.g., knee osteoarthritis (KOA). Accordingly, considerable efforts have been made to develop musculoskeletal finite element (MS-FE) models to estimate highly detailed tissue mechanics that predict cartilage degeneration. However, creating such models is time-consuming and requires advanced expertise. This limits these complex, yet promising MS-FE models to research applications with few participants and makes the models impractical for clinical assessments. Also, these previously developed MS-FE models have not been used to assess activities other than gait. This study introduces and verifies a semi-automated rapid state-of-the-art MS-FE modeling and simulation toolbox incorporating an electromyography- (EMG) assisted MS model and a muscle-force driven FE model of the knee with fibril-reinforced poro(visco)elastic cartilages and menisci. To showcase the usability of the pipeline, we estimated joint- and tissue-level knee mechanics in 15 KOA individuals performing different daily activities. The pipeline was verified by comparing the estimated muscle activations and joint mechanics to existing experimental data. To determine the importance of EMG-assisted MS approach, results were compared to those from the same FE models but driven by static-optimization-based MS models. The EMG-assisted MS-FE pipeline bore a closer resemblance to experiments compared to the static-optimization-based MS-FE pipeline. Importantly, the developed pipeline showed great potential as a rapid MS-FE analysis toolbox to investigate multiscale knee mechanics during different activities of individuals with KOA. The template FE model of the study is freely available here.

Is subject-specific musculoskeletal modelling worth the extra effort or is generic modelling worth the shortcut?

The majority of musculoskeletal modelling studies investigating healthy populations use generic models linearly scaled to roughly match an individual’s anthropometry. Generic models disregard the considerable variation in musculoskeletal geometry and tissue properties between individuals. This study investigated the physiological implications of personalizing musculoskeletal model geometry (body segment mass, inertia, joint center, and maximum isometric muscle force). Nine healthy athletes performed ten repetitions of 15 meter sprints at 75–95% of their maximum sprinting speed and ten repetitions of unanticipated sidestep cut trials with a 4.5–5.5 m/s approach running speed. Structural magnetic resonance imaging was collected on the lower extremities, from which subject-specific musculoskeletal models were developed. A one-dimensional statistical parametric mapping paired t-test was used to compare generic and subject-specific musculoskeletal models for: lower-limb kinematics, kinetics, torque matching, as well as hamstrings, adductors, and quadriceps muscle activations and fiber dynamics. Percentage change of geometric parameters between generic and subject-specific models were determined. Compared to generic models, subject-specific models showed significantly lower ankle dorsi/plantar flexion angle during sprinting and several significantly different net joint moments during sprint and cut tasks. Additionally, subject-specific models demonstrated better torque matching, more physiologically plausible fiber lengths, higher fiber velocities, lower muscle forces, and lower simulated activations in a subset of investigated muscles and motor tasks. Furthermore, subject-specific models identified between-limb differences that were not identified with generic models. Use of subject-specific modeling, even in healthy populations, may result in more physiologically plausible muscle fiber mechanics. Implementing subject-specific models may be especially beneficial when investigating populations with substantial geometric between-limb differences, or unilateral musculoskeletal pathologies, as these are not captured by a generic model.

In vivo assessment of the passive stretching response of the bi-compartmental human semitendinosus muscle using shear wave elastography

The semitendinosus muscle contains distinct proximal and distal compartments arranged anatomically in-series but separated by a tendinous inscription, with each compartment innervated by separate nerve branches. Although extensively investigated in other mammals, compartment-specific mechanical properties within the human semitendinosus have scarcely been assessed in vivo. Experimental data obtained during muscle-tendon unit stretching (e.g., slack angle) can also be used to validate and/or improve musculoskeletal model estimates of semitendinosus muscle force. The purpose of this study was to investigate the passive stretching response of proximal and distal humans semitendinosus compartments to distal joint extension. Using two-dimensional shear wave elastography, we bilaterally obtained shear moduli of both semitendinosus compartments from 14 prone-positioned individuals at ten knee flexion angles (from 90° to 0° [full extension] at 10° intervals). Passive muscle mechanical characteristics (slack angle, slack shear modulus, and the slope of the increase in shear modulus) were determined for each semitendinosus compartment by fitting a piecewise exponential model to the shear modulus-joint angle curves. We found no differences between compartments or legs for slack angle, slack shear modulus, or the slope of the increase in shear modulus. We also found the experimentally determined slack angle occurred at ~15-80° higher knee flexion angles compared to estimates from two commonly used musculoskeletal models, depending on participant and model used. Overall, these findings demonstrate that passive shear modulus-joint angle curves do not differ between proximal and distal human semitendinosus compartments, and provide experimental data to improve semitendinosus force estimates derived from musculoskeletal models.

Finite element analysis of the performance of additively manufactured scaffolds for scapholunate ligament reconstruction

Rupture of the scapholunate interosseous ligament can cause the dissociation of scaphoid and lunate bones, resulting in impaired wrist function. Current treatments (e.g., tendon-based surgical reconstruction, screw-based fixation, fusion, or carpectomy) may restore wrist stability, but do not address regeneration of the ruptured ligament, and may result in wrist functional limitations and osteoarthritis. Recently a novel multiphasic bone-ligament-bone scaffold was proposed, which aims to reconstruct the ruptured ligament, and which can be 3D-printed using medical-grade polycaprolactone. This scaffold is composed of a central ligament-scaffold section and features a bone attachment terminal at either end. Since the ligament-scaffold is the primary load bearing structure during physiological wrist motion, its geometry, mechanical properties, and the surgical placement of the scaffold are critical for performance optimisation. This study presents a patient-specific computational biomechanical evaluation of the effect of scaffold length, and positioning of the bone attachment sites. Through segmentation and image processing of medical image data for natural wrist motion, detailed 3D geometries as well as patient-specific physiological wrist motion could be derived. This data formed the input for detailed finite element analysis, enabling computational of scaffold stress and strain distributions, which are key predictors of scaffold structural integrity. The computational analysis demonstrated that longer scaffolds present reduced peak scaffold stresses and a more homogeneous stress state compared to shorter scaffolds. Furthermore, it was found that scaffolds attached at proximal sites experience lower stresses than those attached at distal sites. However, scaffold length, rather than bone terminal location, most strongly influences peak stress. For each scaffold terminal placement configuration, a basic metric was computed indicative of bone fracture risk. This metric was the minimum distance from the bone surface to the internal scaffold bone terminal. Analysis of this minimum bone thickness data confirmed further optimisation of terminal locations is warranted.

Patellar cartilage increase following ACL reconstruction with and without meniscal pathology: a two-year prospective MRI morphological study

Background

Anterior cruciate ligament reconstruction (ACLR) together with concomitant meniscal injury are risk factors for the development of tibiofemoral (TF) osteoarthritis (OA), but the potential effect on the patellofemoral (PF) joint is unclear. The aim of this study was to: (i) investigate change in patellar cartilage morphology in individuals 2.5 to 4.5 years after ACLR with or without concomitant meniscal pathology and in healthy controls, and (ii) examine the association between baseline patellar cartilage defects and patellar cartilage volume change.

Methods

Thirty two isolated ACLR participants, 25 ACLR participants with combined meniscal pathology and nine healthy controls underwent knee magnetic resonance imaging (MRI) with 2-year intervals (baseline = 2.5 years post-ACLR). Patellar cartilage volume and cartilage defects were assessed from MRI using validated methods.

Results

Both ACLR groups showed patellar cartilage volume increased over 2 years (p < 0.05), and isolated ACLR group had greater annual percentage cartilage volume increase compared with controls (mean difference 3.6, 95% confidence interval (CI) 1.0, 6.3%, p = 0.008) and combined ACLR group (mean difference 2.2, 95% CI 0.2, 4.2%, p = 0.028). Patellar cartilage defects regressed in the isolated ACLR group over 2 years (p = 0.02; Z = − 2.33; r = 0.3). Baseline patellar cartilage defect score was positively associated with annual percentage cartilage volume increase (Regression coefficient B = 0.014; 95% CI 0.001, 0.027; p = 0.03) in the pooled ACLR participants.

Conclusions

Hypertrophic response was evident in the patellar cartilage of ACLR participants with and without meniscal pathology. Surprisingly, the increase in patellar cartilage volume was more pronounced in those with isolated ACLR. Although cartilage defects stabilised in the majority of ACLR participants, the severity of patellar cartilage defects at baseline influenced the magnitude of the cartilage hypertrophic response over the subsequent ~ 2 years.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-021-04794-5.

Effects of Pubertal Maturation on ACL Forces During a Landing Task in Females

BACKGROUND

Rates of anterior cruciate ligament (ACL) rupture in young people have increased by >70% over the past two decades. Adolescent and young adult females are at higher risk of ACL injury as compared with their prepubertal counterparts.

PURPOSE

To determine ACL loading during a standardized drop-land-lateral jump in females at different stages of pubertal maturation.

STUDY DESIGN

Controlled laboratory study.

METHODS

On the basis of the Tanner classification system, 19 pre-, 19 early-/mid-, and 24 late-/postpubertal females performed a standardized drop-land-lateral jump while 3-dimensional body motion, ground-reaction forces, and surface electromyography data were acquired. These data were used to model external biomechanics, lower limb muscle forces, and knee contact forces, which were subsequently used in a validated computational model to estimate ACL loading. Statistical parametric mapping analysis of variance was used to compare ACL force and its causal contributors among the 3 pubertal maturation groups during stance phase of the task.

RESULTS

When compared with pre- and early-/midpubertal females, late-/postpubertal females had significantly higher ACL force with mean differences of 471 and 356 N during the first 30% and 48% to 85% of stance, and 343 and 274 N during the first 24% and 59% to 81% of stance, respectively, which overlapped peaks in ACL force. At the point of peak ACL force, contributions from sagittal and transverse plane loading mechanisms to ACL force were higher in late-/postpubertal compared with pre- and early-/midpubertal groups (medium effect sizes from 0.44 to 0.77). No differences were found between pre- and early-/midpubertal groups in ACL force or its contributors.

CONCLUSION

The highest ACL forces were observed in late-/postpubertal females, consistent with recently reported rises of ACL injury rates in females aged 15 to 19 years. It is important to quantify ACL force and its contributors during dynamic tasks to advance our understanding of the loading mechanism and thereby provide guidance to injury prevention.

CLINICAL RELEVANCE

Growth of ACL volume plateaus around 10 years of age, before pubertal maturation, meaning that a late-/postpubertal female could have an ACL of similar size to their less mature counterparts. However, late-/postpubertal females have higher body mass requiring higher muscle forces to accelerate the body during dynamic tasks, which may increase ACL loading. Thus, if greater forces develop in these females, in part because of their increased body mass, these higher forces will be applied to an ACL that is not proportionally larger. This may partially explain the higher rates of ACL injury in late-/postpubertal females.

The deep hip muscles are unlikely to contribute to hip stability in the sagittal plane during walking: a stiffness approach

- Objective: This study determined whether the deep hip muscles could contribute to hip stability.

METHODS

Hip stability was defined as rotational hip stiffness in the sagittal plane, which was calculated for walking trials for 12 participants via electromyography (EMG)informed neuromusculoskeletal modelling that included all 22 hip spanning muscles. Three model configurations were compared that differed in the excitations of the deep hip muscles, but were identical in the excitations of all other muscles: (1) deep hip muscles informed by intramuscular EMG measurements (assisted activation); (2) deep hip muscles with simulated zero activation (no activation); (3) deep hip muscles with simulated maximal activation (maximal activation). Sagittal plane rotational hip stiffness over the gait cycle was compared between model configurations using a within-participant analysis of variance via statistical parametric mapping (p<0.05).

RESULTS

Compared to the assisted activation configuration, hip stiffness (mean (95% confidence interval)) was 0.8% (0.7 to 0.9) lower in the no activation configuration, and 3.2% (2.9 to 3.5) higher in the maximal activation configuration over the gait cycle.

CONCLUSION

Regardless of activation level, the deep hip muscles made little contribution to sagittal plane rotational hip stiffness, which casts uncertainty around their assumed function as hip stabilizers.

SIGNIFICANCE

The merit of targeted deep hip muscle strengthening to improve hip stability in rehabilitation programs for remains unclear.

Reliability of hip muscle strength measured in principal and intermediate planes of movement

BACKGROUND

Muscle strength testing is widely used in clinical and athletic populations. Commercially available dynamometers are designed to assess strength in three principal planes (sagittal, transverse, frontal). However, the anatomy of the hip suggests muscles may only be recruited submaximally during tasks performed in these principal planes.

OBJECTIVE

To evaluate the inter-session reliability of maximal isometric hip strength in the principal planes and three intermediate planes.

METHODS

Twenty participants (26.1 ± 2.7 years, 50% female) attended two testing sessions 6.2 ± 1.8 days apart. Participants completed 3-5 maximal voluntary isometric contractions for hip abduction, adduction, flexion, extension, and internal and external rotation measured using a fixed uniaxial load cell (custom rig) and commercial dynamometer (Biodex). Three intermediate hip actions were also tested using the custom rig: extension with abduction, extension with external rotation, and extension with both abduction and external rotation.

RESULTS

Moderate-to-excellent intraclass correlation coefficients were observed for all principal and intermediate muscle actions using the custom rig (0.72-0.95) and the Biodex (0.85-0.95). The minimum detectable change was also similar between devices (custom rig = 11-31%; Biodex = 9-20%). Bland-Altman analysis revealed poor agreement between devices (range between upper and lower limits of agreement = 77-131%).

CONCLUSIONS

Although the custom rig and Biodex showed similar reliability, both devices may lack the sensitivity to detect small changes in hip strength commonly observed following intervention.

Mechanism of Anterior Cruciate Ligament Loading during Dynamic Motor Tasks

INTRODUCTION

This study determined anterior cruciate ligament (ACL) force and its contributors during a standardized drop-land-lateral jump task using a validated computational model.

METHODS

Three-dimensional whole-body kinematics, ground reaction forces, and muscle activation patterns from eight knee-spanning muscles were collected during dynamic tasks performed by healthy recreationally active females (n = 24). These data were used in a combined neuromusculoskeletal and ACL force model to determine lower limb muscle and ACL forces.

RESULTS

Peak ACL force (2.3 ± 0.5 bodyweight) was observed at ~14% of stance during the drop-land-lateral jump. The ACL force was primarily generated through the sagittal plane, and muscle was the dominant source of ACL loading. The main ACL antagonists (i.e., loaders) were the gastrocnemii and quadriceps, whereas the hamstrings were the main ACL agonists (i.e., supporters).

CONCLUSION

Combining neuromusculoskeletal and ACL force models, the roles of muscle in ACL loading and support were determined during a challenging motor task. Results highlighted the importance of the gastrocnemius in ACL loading, which could be considered more prominently in ACL injury prevention and rehabilitation programs.

Lower Limb Biomechanical Responses During a Standardized Load Carriage Task are Sex Specific

INTRODUCTION

The purpose of this study was to investigate sex-specific lower limb biomechanical adaptations during a standardized load carriage task in response to a targeted physical training program.

MATERIALS AND METHODS

Twenty-five healthy civilians (males [n = 13] and females [n = 12]) completed a load carriage task (5 km at 5.5 km·h-1, wearing a 23 kg vest) before and after a 10-week lower-body-focused training program. Kinematics and ground reaction force data were collected during the task and were used to estimate lower limb joint kinematics and kinetics (i.e., moments and powers). Direct statistical comparisons were not conducted due to different data collection protocols between sexes. A two-way repeated measures ANOVA tested for significant interactions between, and main effects of training and distance marched for male and female data, respectively.

RESULTS

Primary kinematic and kinetic changes were observed at the knee and ankle joints for males and at the hip and knee joints for females. Knee joint moments increased for both sexes over the 5 km distance marched (P > .05), with males demonstrating significant reductions in peak knee joint extension after training. Hip adduction, internal rotation, and knee internal rotation angles significantly increased after the 5 km load carriage task for females but not males.

CONCLUSION

Differences in adaptive gait strategies between sexes indicate that physical training needs to be tailored to sex-specific requirements to meet standardized load carriage task demands. The findings highlighted previously unfound sex-specific responses that could inform military training and facilitate the integration of female soldiers into physically demanding military roles.

Repeated sprints alter mechanical work done by hip and knee, but not ankle, sagittal moments

OBJECTIVES

To quantify the changes in work done by lower limb joint moments during maximal speed running following a sports-specific repeated running protocol.

DESIGN

Observational with repeated-measures.

METHODS

Recreational athletes (n = 18 (9 females), aged = 26.2 ± 6.2 years) performed 12 maximal 30-m sprints on a non-motorised treadmill. Three-dimensional kinematics and ground reaction forces were subsequently recorded during a 10-m maximal overground sprint before and immediately after the repeated running protocol, from which we calculated work done by sagittal plane hip, knee, and ankle moments. Relative work (J/kg) was reported as a percentage of positive and negative work done by the sum of joint moments.

RESULTS

Following the repeated running protocol, maximal sprint speed decreased by 19% and was accompanied by reductions in total positive (-1.47 J/kg) and negative (-0.92 J/kg) work, in addition to work done by hip (-0.43 to -0.82 J/kg) and knee (-0.28 J/kg) moments during swing. Compared to before the repeated running protocol, less relative work was done by hip (-9%) and knee (-3%) extension moments during swing. Reductions in work done by hip and knee joint moments during swing were significantly correlated with reductions in maximum running speed (r = 0.61-0.89, p < 0.05).

CONCLUSIONS

A sports-specific repeated running protocol resulted in reductions in mechanical work done by sagittal plane hip and knee joint moments during maximal overground sprinting. Interventions focused on maintaining positive work done by the hip flexors/extensors and negative work done by knee flexors/extensors during the swing phase of running may help prevent reductions in speed following repeated sprinting.

Effect of exercise on knee joint contact forces in people following medial partial meniscectomy: A secondary analysis of a randomised controlled trial

BACKGROUND

Arthroscopic partial meniscectomy may cause knee osteoarthritis, which may be related to altered joint loading. Previous research has failed to demonstrate that exercise can reduce medial compartment knee loads following meniscectomy but has not considered muscular loading in their estimates.

RESEARCH QUESTION

What is the effect of exercise compared to no intervention on peak medial tibiofemoral joint contact force during walking using an electromyogram-driven neuromusculoskeletal model, following medial arthroscopic partial meniscectomy?

METHODS

This is a secondary analysis of a randomized controlled trial (RCT). 41 participants aged between 30-50 years with medial arthroscopic partial meniscectomy within the past 3-12 months, were randomly allocated to either a 12-week, home-based, physiotherapist-guided exercise program or to no exercise (control group). Three-dimensional lower-body motion, ground reaction forces, and surface electromyograms from eight lower-limb muscles were acquired during self-selected normal- and fast-paced walking at baseline and follow-up. An electromyogram-driven neuromusculoskeletal model estimated medial compartment contact forces (body weight). Linear regression models evaluated between-group differences (mean difference (95% CI)).

RESULTS

There were no significant between-group differences in the change (follow-up minus baseline) in first peak medial contact force during self-selected normal- or fast-paced walking (0.07 (-0.08 to 0.23), P = 0.34 and 0.01 (-0.19 to 0.22), P = 0.89 respectively). No significant between-group difference was found for change in second peak medial contact force during normal- or fast-paced walking (0.09 (-0.09 to 0.28), P = 0.31 and 0.02 (-0.17 to 0.22), P = 0.81 respectively). At the individual level, variability was observed for changes in first (range -26.2% to +31.7%) and second (range -46.5% to +59.9%) peak tibiofemoral contact force.

SIGNIFICANCE

This is the first study to apply electromyogram-driven neuromusculoskeletal modelling to an exercise intervention in a RCT. While our results suggest that a 12-week exercise program does not alter peak medial knee loads after meniscectomy, within-participant variability suggests individual-specific muscle activation patterns that warrant further investigation.

Development and validation of statistical shape models of the primary functional bone segments of the foot

Introduction

Musculoskeletal models are important tools for studying movement patterns, tissue loading, and neuromechanics. Personalising bone anatomy within models improves analysis accuracy. Few studies have focused on personalising foot bone anatomy, potentially incorrectly estimating the foot’s contribution to locomotion. Statistical shape models have been created for a subset of foot-ankle bones, but have not been validated. This study aimed to develop and validate statistical shape models of the functional segments in the foot: first metatarsal, midfoot (second-to-fifth metatarsals, cuneiforms, cuboid, and navicular), calcaneus, and talus; then, to assess reconstruction accuracy of these shape models using sparse anatomical data.

Methods

Magnetic resonance images of 24 individuals feet (age = 28 ± 6 years, 52% female, height = 1.73 ± 0.8 m, mass = 66.6 ± 13.8 kg) were manually segmented to generate three-dimensional point clouds. Point clouds were registered and analysed using principal component analysis. For each bone segment, a statistical shape model and principal components were created, describing population shape variation. Statistical shape models were validated by assessing reconstruction accuracy in a leave-one-out cross validation. Statistical shape models were created by excluding a participant’s bone segment and used to reconstruct that same excluded bone using full segmentations and sparse anatomical data (i.e. three discrete points on each segment), for all combinations in the dataset. Tali were not reconstructed using sparse anatomical data due to a lack of externally accessible landmarks. Reconstruction accuracy was assessed using Jaccard index, root mean square error (mm), and Hausdorff distance (mm).

Results

Reconstructions generated using full segmentations had mean Jaccard indices between 0.77 ± 0.04 and 0.89 ± 0.02, mean root mean square errors between 0.88 ± 0.19 and 1.17 ± 0.18 mm, and mean Hausdorff distances between 2.99 ± 0.98 mm and 6.63 ± 3.68 mm. Reconstructions generated using sparse anatomical data had mean Jaccard indices between 0.67 ± 0.06 and 0.83 ± 0.05, mean root mean square error between 1.21 ± 0.54 mm and 1.66 ± 0.41 mm, and mean Hausdorff distances between 3.21 ± 0.94 mm and 7.19 ± 3.54 mm. Jaccard index was higher (P < 0.01) and root mean square error was lower (P < 0.01) in reconstructions from full segmentations compared to sparse anatomical data. Hausdorff distance was lower (P < 0.01) for midfoot and calcaneus reconstructions using full segmentations compared to sparse anatomical data.

Conclusion

For the first time, statistical shape models of the primary functional segments of the foot were developed and validated. Foot segments can be reconstructed with minimal error using full segmentations and sparse anatomical landmarks. In future, larger training datasets could increase statistical shape model robustness, extending use to paediatric or pathological populations.

Modelling the loading mechanics of anterior cruciate ligament

BACKGROUND AND OBJECTIVES

The anterior cruciate ligament (ACL) plays a crucial role in knee stability and is the most commonly injured knee ligament. Although ACL loading patterns have been investigated previously, the interactions between knee loadings transmitted to ACL remain elusive. Understanding the loading mechanism of ACL during dynamic tasks is essential to prevent ACL injuries. Therefore, we propose a computational model that predicts the force applied to ACL in response to knee loading in three planes of motion.

METHODS

First, a three-dimensional (3D) computational model was developed and validated using available cadaveric experimental data to predict ACL force. This 3D model was then combined with a neuromusculoskeletal model of lower limb and used to estimate in vivo ACL forces during a standardised drop-landing task. The neuromusculoskeletal model utilised movement data collected from female participants during a dynamic task and calculated lower limb joint kinematics and kinetics, as well as muscle forces.

RESULTS

The total ACL force predicted by the 3D computational ACL force model was in good agreement with cadaveric data, as strong correlation (r² = 0.96 and P < 0.001), minimal bias, and narrow limits of agreement were observed. The combined model further illustrated that the ACL is primarily loaded through the sagittal plane, mainly due to muscle loading.

CONCLUSIONS

The proposed computational model is the first validated model that can provide an accessible tool to develop and test knee ACL injury prevention programs for people with normal ACL. This method can be extended to study the abnormal ACL upon the availability of relevant experimental data.

Reliability and Validity of Ultrasonography for Measurement of Hamstring Muscle and Tendon Cross-Sectional Area

The purpose of this study was to determine the reliability and validity of ultrasonography for measurement of hamstring muscle and semitendinosus (ST) tendon cross-sectional area (CSA). On two consecutive days, muscle anatomical CSA (ACSA) and ST tendon CSA were measured at standardized positions (30%-80% of thigh length; half the distance from the distal muscle-tendon junction to the popliteal crease) on 12 legs using ultrasonography and compared with corresponding magnetic resonance imaging measures. Inter-day intraclass correlation coefficients were good-to-excellent (0.882-0.996) for all assessed muscle and tendon sites. The limits of agreement widths were narrowest (range: 17%-52%) when muscle ACSA was large but were wide at sites with relatively small ACSA (≤184%) and for ST tendon CSA (range: 72%). Results suggest ultrasound-based measures of individual hamstring muscle maximal ACSA are reliable and valid and ST tendon CSA measures are reliable but require comparison with cadaveric or intra-operative measurements to verify validity.

Neuromusculoskeletal Modeling-Based Prostheses for Recovery After Spinal Cord Injury

Concurrent stimulation and reinforcement of motor and sensory pathways has been proposed as an effective approach to restoring function after developmental or acquired neurotrauma. This can be achieved by applying multimodal rehabilitation regimens, such as thought-controlled exoskeletons or epidural electrical stimulation to recover motor pattern generation in individuals with spinal cord injury (SCI). However, the human neuromusculoskeletal (NMS) system has often been oversimplified in designing rehabilitative and assistive devices. As a result, the neuromechanics of the muscles is seldom considered when modeling the relationship between electrical stimulation, mechanical assistance from exoskeletons, and final joint movement. A powerful way to enhance current neurorehabilitation is to develop the next generation prostheses incorporating personalized NMS models of patients. This strategy will enable an individual voluntary interfacing with multiple electromechanical rehabilitation devices targeting key afferent and efferent systems for functional improvement. This narrative review discusses how real-time NMS models can be integrated with finite element (FE) of musculoskeletal tissues and interface multiple assistive and robotic devices with individuals with SCI to promote neural restoration. In particular, the utility of NMS models for optimizing muscle stimulation patterns, tracking functional improvement, monitoring safety, and providing augmented feedback during exercise-based rehabilitation are discussed.