Introduction. Adolescent Idiopathic Scoliosis (AIS) is a three-dimensional deformity of the spine with unclear etiology. Due to the asymmetry of lateral curves, there are differences in the muscle activation between the convex and concave sides. This study utilized a comprehensive thoracic spine and ribcage musculoskeletal model to improve the biomechanical understanding of the development of AIS deformity and approach an explanation of the condition. Methods. In this study, we implemented a motion capture model using a generic rigid-body thoracic spine and ribcage model, which is kinematically determinate and controlled by spine posture obtained, for instance, from radiographs. This model is publicly accessible via a GitHub repository. We simulated gait and standing models of two AIS (averaging 15 years old, both with left lumbar curve and right thoracic curve averaging 25 degrees) and one control subject. The marker set included extra markers on the sternum and the thoracic and lumbar spine. The study was approved by the regional Research Ethics Committee (Journal number: H17034237). Results. We investigated the difference between the muscle activation on the right and left sides including erector spinae (ES), psoas major (PS), and multifidus (MF). Results of the AIS simulations indicated that, on average throughout the gait cycle, the right ES, left PS and left MF had 46%, 44%, and 23% higher activities compared to the other side, respectively. In standing, the ratios were 28%, 40%, and 19%, respectively. However, for the control subject, the differences were under 7%, except ES throughout the gait, which was 17%. Conclusion. The musculoskeletal model revealed distinct differences in force patterns of the right and left sides of the spine, indicating an instability phenomenon, where larger curves lead to higher muscle activations for stabilization. Acknowledgement. The project is funded by the European Union's Horizon 2020 program through Marie Skłodowska-Curie grant No. [764644]

For a proper rehabilitation of the knee following knee arthroplasty, a comprehensive understanding of bony and soft tissue structures and their effects on biomechanics of the individual patient is essential. Musculoskeletal models have the potential, however, to predict dynamic interactions of the knee joint and provide knowledge to the understanding of knee biomechanics. Our goal was to develop a generic musculoskeletal knee model which is adaptable to subject-specific situations and to use in-vivo kinematic measurements obtained under full-weight bearing condition in a previous Upright-MRI study of our group for a proper validation of the simulation results. The simulation model has been developed and adapted to subject-specific cases in the multi-body simulation software AnyBody. For the implementation of the knee model a reference model from the AnyBody Repository was adapted for the present issue. The standard hinge joint was replaced with a new complex knee joint with 6DoF. The 3D bone geometries were obtained from an optimized MRI scan and then post-processed in the mesh processing software MeshLab. A homogenous dilation of 3 mm was generated for each bone and used as articulating surfaces. The anatomical locations of viscoelastic ligaments and muscle attachments were determined based on literature data. Ligament parameters, such as elongation and slack length, were adjusted in a calibration study in two leg stance as reference position. For the subject-specific adaptation a general scaling law, taking segment length, mass and fat into account, was used for a global scaling. The scaling law was further modified to allow a detailed adaption of the knee region, e.g. align the subject-specific knee morphology (including ligament and muscle attachments) in the reference model. The boundary conditions were solely described by analytical methods since body motion (apart from the knee region) or force data were not recorded in the Upright-MRI study. Ground reaction forces have been predicted and a single leg deep knee bend was simulated by kinematic constraints, such as that the centre of mass is positioned above the ankle joint. The contact forces in the knee joint were computed using the force dependent kinematic algorithm. Finally, the simulation model was adapted to three subjects, a single leg deep knee bend was simulated, subject-specific kinematics were recorded and then compared to their corresponding in-vivo kinematic measurements data. We were able to simulate the whole group of subjects over the complete range of motion. The tibiofemoral kinematics of three subjects could be simulated showing the overall trend correctly, whereas absolute values partially differ. In conclusion, the presented simulation model is highly adaptable to an individual situation and seems to be suitable to approximate subject-specific knee kinematics without consideration of cartilage and menisci. The model enables sensitivity analyses regarding changes in patient specific knee kinematics following e.g. surgical interventions on bone or soft tissue as well as related to the design and placement of partial or total knee joint replacement. However, model optimisation, a higher case number, sensitivity analyses of selected parameters and a semi-automation of the workflow are parts of our ongoing work

Aims. The surgical target for optimal implant positioning in robotic-assisted total knee arthroplasty remains the subject of ongoing discussion. One of the proposed targets is to recreate the knee’s functional behaviour as per its pre-diseased state. The aim of this study was to optimize implant positioning, starting from mechanical alignment (MA), toward restoring the pre-diseased status, including ligament strain and kinematic patterns, in a patient population. Methods. We used an active appearance model-based approach to segment the preoperative CT of 21 osteoarthritic patients, which identified the osteophyte-free surfaces and estimated cartilage from the segmented bones; these geometries were used to construct patient-specific musculoskeletal models of the pre-diseased knee. Subsequently, implantations were simulated using the MA method, and a previously developed optimization technique was employed to find the optimal implant position that minimized the root mean square deviation between pre-diseased and postoperative ligament strains and kinematics. Results. There were evident biomechanical differences between the simulated patient models, but also trends that appeared reproducible at the population level. Optimizing the implant position significantly reduced the maximum observed strain root mean square deviations within the cohort from 36.5% to below 5.3% for all but the anterolateral ligament; and concomitantly reduced the kinematic deviations from 3.8 mm (SD 1.7) and 4.7° (SD 1.9°) with MA to 2.7 mm (SD 1.4) and 3.7° (SD 1.9°) relative to the pre-diseased state. To achieve this, the femoral component consistently required translational adjustments in the anterior, lateral, and proximal directions, while the tibial component required a more posterior slope and varus rotation in most cases. Conclusion. These findings confirm that MA-induced biomechanical alterations relative to the pre-diseased state can be reduced by optimizing the implant position, and may have implications to further advance pre-planning in robotic-assisted surgery in order to restore pre-diseased knee function. Cite this article: Bone Joint J 2024;106-B(11):1231–1239

Head collisions in sport can result in catastrophic cervical spine injuries. Musculo-skeletal (MSK) modelling can help analyse the relationship between players' motion, external loading and internal stresses that lead to injury. However, the literature lacks sport specific MSK models. In automotive research the intervertebral disc behaviour has been represented as viscoelastic elements (“bushing”), whose stiffness and damping parameters are often estimated under quasi-static conditions and may lack validity in dynamic impacts. The aim of this study was to develop a validated cervical spine model for axial impacts for future use in the analysis of head-first rugby collisions.

A drop test rig was used to replicate a sub-catastrophic axial head impact. A load of 80 N from 0.5 m was applied to the cranial aspect of a C2-C6 porcine spinal specimen mounted in the neutral position. The 3D motion of C3-C5 vertebras (4 kHz) and the cranial axial load (1 MHz) were measured via motion capture (Qualysis, Sweden) and a uniaxial load cell (RDP Electronics Ltd, UK). Specimen specific models were created in NMSBuilder and OpenSim after the vertebrae geometries were obtained from the segmentation of micro-CT images of the specimens. The compressive viscoelastic properties of four vertebral joints (C2-C3 through to C5-C6) were optimised via a Genetic Algorithm (MATLAB v2016b, The Mathworks Inc) to minimise tracking errors.

The optimisation converged to a solution of 140–49000 kN/m and 2000–8000 Ns/m for stiffness and damping respectively (RMSE=5.1 mm). Simulated joint displacements ranged between 0.09 – 1.75 mm compared to experimental 0.1 – 0.8 mm.

Optimal bushing parameters were higher than previously reported values measured through quasi-static testing. Higher stiffness and damping values could be explained by the higher-dynamics nature of the event analysed related to a different part of the non-linear intervertebral disc load-displacement curve.

Purpose: Previous work demonstrated that negative pressure wound therapy (NPWT) resulted in less Pseudomonas aeruginosa than standard wet-to-dry (WTD) dressings in a complex orthopaedic wound model. Staphylococcus aureus is more clinically relevant in open fractures, and is the most prevalent bacteria in osteomyelitis. The purpose of this study is to determine if S. aureus responds similarly to P. aeruginosa when treated with NPWT.

Methods: A complex musculoskeletal wound was created on the hindlimb of 20 goats and contaminated with S. aureus (lux) bacteria. The bacteria are genetically engineered to emit photons, allowing for quantification with a photon-counting camera system. The wounds were débrided and irrigated with 9 L of normal saline using gravity flow irrigation 6 hours after inoculation. Goats were assigned to two different treatment groups: a control group using WTD dressing changes and an experimental group using NPWT. The wounds were débrided and irrigated every other day for 6 days. Bacteria within the wounds were quantified both before and after each débridement.

Results: There was no difference between treatment groups in amounts of bacteria in the wound at all time points (p≥0.37).

Conclusion: Previous work demonstrated that NPWT resulted in a significant and clinically relevant reduction of P. aeruginosa at all time points in a similar model. We presume that NPWT was effective because it created an environment that allowed the body to ward off this “opportunistic” gram negative. However, as shown in this study, S. aureus is less affected by NPWT and persists within the wound.

Introduction. This research aims to enhance the control of intricate musculoskeletal spine models, a critical tool for comprehending both healthy and pathological spinal conditions. State-of-the-art musculoskeletal spine models incorporate segments for all vertebra, each possessing 3 degrees-of-freedom (DOF). Manually defining the posture with this amount of DOFs presents a significant challenge. The prevalent method of equally distributing the spine's overall rotation among the vertebrae often proves to be an inadequate assumption, particularly when dealing with the entire spine. Method. We have engineered a comprehensive non-linear spine rhythm and the requisite tools for its implementation in widely utilized musculoskeletal modelling software (1). The rhythm controls lateral bending, axial rotation, and flexion/extension. The mathematical and implementation details of the rhythm are beyond this abstract, but it's noteworthy that the implementation accommodates non-linear rhythms. This means, for example, that one set of rhythm coefficients is used for flexion and another for extension. The rhythm coefficients, which distinguish the movement along the spine, were derived from a review of spine literature. The values for spine and vertebra range-of-motion (ROM) vary significantly in published studies, and no complete dataset was found in any single study. Consequently, the rhythm presented here is a composite, designed to provide the most consistent and average set of rhythm coefficients. Result. The novel spine rhythm simplifies the control of detailed spine models, accommodating varying amounts of input data. It allows for the specification of only the overall motion or the posture at a more detailed level (i.e., lumbar, thoracic, neck). The tools and rhythm coefficients are publicly available on GitHub. Conclusion. The innovative spine rhythm enhances the usability of cutting-edge spine models. For flexion/extension of the spine, it introduces a non-linear rhythm, exhibiting distinct behaviour between flexion and extension - a feature not previously observed in musculoskeletal spine models. 1) The AnyBody Modeling System

Abstract. OBJECTIVE. Knee varus malalignment increases medial knee compartment loading and is associated with knee osteoarthritis (OA) progression and severity. 1. Altered biomechanical loading and dysregulation of joint tissue biology drive OA progression, but mechanistic links between these factors are lacking. Subchondral bone structural changes are biomechanically driven, involve bone resorption, immune cell influx, angiogenesis, and sensory nerve invasion, and contribute to joint destruction and pain. 2. We have investigated mechanisms underlying this involving RANKL and alkaline phosphatase (ALP), which reflect bone resorption and mineralisation respectively. 3. and the axonal guidance factor Sema3A. Sema3A is osteotropic, expressed by mechanically sensitive osteocytes, and an inhibitor of sensory nerve, blood vessel and immune cell invasion. 4. Sema3A is also differentially expressed in human OA bone. 5. HYPOTHESIS: Medial knee compartment overloading in varus knee malalignment patients causes dysregulation of bone derived Sema3A signalling directly linking joint biomechanics to pathology and pain. METHODS. Synovial fluid obtained from 30 subjects with medial knee OA (KL grade II-IV) undergoing high tibial osteotomy surgery (HTO) was analysed by mesoscale discovery and ELISA analysis for inflammatory, neural and bone turnover markers. 11 of these patients had been previously analysed in a published patient-specific musculoskeletal model. 6. of gait estimating joint contact location, pressure, forces, and medial-lateral condyle load distribution in a published data set included in analyses. Data analysis was performed using Pearson's correlation matrices and principal component analyses. Principal Components (PCs) with eigenvalues greater than 1 were analysed. RESULTS. PC1 (32.94% of variation) and PC2 (25.79% of variation) from PCA analysis and correlation matrices separated patients according to correlated clusters of established inflammatory markers of OA pain and progression (IL6/IL8, r=0.754, p<0.001) and anti-inflammatory mediators (IL4/IL10, r=0.469, p=0.005). Bone turnover marker ALP was positively associated with KL grade (r=0.815, p=0.002) and negatively associated with IL10 (r=−0.402, p=0.018) and first peak knee loading pressures (r=−0.688, p=0.019). RANKL was positively associated with IL4 (r=0.489, p=0.003). Synovial fluid Sema3A concentrations showed separate clustering from all OA progression markers and was inversely correlated with TNF-α (r=−0.423, p=0.022) in HTO patients. Sema3A was significantly inversely correlated with total predicted force in the medial joint compartment (r=−0.621, p=0.041), mean (r=−0.63, p=0.038) and maximum (r=−0.613, p=0.045) calculated medial compartment joint pressures during the first phase and mean (r=−0.618, p=0.043) and maximum (r=−0.641, p=0.034) medial compartment joint pressures during midstance outputs of patient-specific musculoskeletal model. CONCLUSIONS. This study shows joint inflammatory status and mechanical overloading influence subchondral bone-remodelling. Synovial Sema3A concentrations are inversely correlated to patient-specific musculoskeletal model estimations of pathological medial overloading. This study reveals Sema3A as a biological mediator with capacity to induce OA pain and disease progression that is directly regulated by gait mechanical loading. Declaration of Interest. (b) declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported:I declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research project

Abstract. Objectives. Osteoporosis of the pelvis and femur is diagnosed in a high proportion of lower-limb amputees which carries an increased fracture risk and subsequently serious implications on mobility, physical dependency and morbidity. Through the development of biofidelic musculoskeletal and finite element (FE) models, we aim to determine the effect of lower-limb amputation on long-term bone remodelling in the hip and to understand the potential underpinning mechanisms for bone degradation in the younger amputee population. Methods. Our models are patient specific and anatomically accurate. Geometries are derived from MRI-scans of one bilateral, above-knee, amputee and one body-matched control subject. Musculoskeletal modelling enables comparison of muscle and joint reaction-forces throughout gait. This provides the loading scenario implemented in FE. FE modelling demonstrates the effect of loading on the amputated limb via a prosthetic socket by comparing bone mechanical stimulation in amputee and control cases. Results. Musculoskeletal modelling shows that the bilateral amputee has 25% higher peak hip-reaction force than controls but a 54% lower peak knee-reaction force. Compensation for missing muscles and joints cause large-scale changes to the muscle loading patterns of the residual limb. FE analysis shows a 32% reduction in bone stimulation within the proximal femur and an 81% reduction in the distal femoral shaft when compared to the healthy control. A shielding effect from weight-bearing through a prosthetic socket was observed that may offset any increases in joint and muscle loading at the amputated hip. Conclusions. Bone loss in the young amputee population could be driven by unloading osteopenia where altered joint and muscle loads cause altered mechanical stimulus in the femur. Over many cycles of remodelling, a net bone loss occurs. Importantly, this suggests that the issue is preventable, or even reversible, with the implementation of targeted loading regimes or changes to the design of the prosthetic socket. Declaration of Interest. (b) declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported:I declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research project

Introduction. Fretting corrosion at the taper interface of modular connections can be studied using Finite Element (FE) analyses. However, the loading conditions in FE studies are often simplified, or based on generic activity patterns. Using musculoskeletal modeling, subject-specific muscle and joint forces can be calculated, which can then be applied to a FE model for wear predictions. The objective of the current study was to investigate the effect of incorporating more detailed activity patterns on fretting simulations of modular connections. Methods. Using a six-camera motion capture system, synchronized force plates, and 45 optical markers placed on 6 different subjects, data was recorded for three different activities: walking at a comfortable speed, chair rise, and stair climbing. Musculoskeletal models, using the Twente Lower Extremity Model 2.0 implemented in the AnyBody modeling System™ (AnyBody Technology A/S, Aalborg, Denmark; figure1), were used to determine the hip joint forces. Hip forces for the subject with the lowest and highest peak force, as well as averaged hip forces were then applied to an FE model of a modular taper connection (Biomet Type-1 taper with a Ti6Al4V Magnum +9 mm adaptor; Figure 2). During the FE simulations, the taper geometry was updated iteratively to account for material removal due to wear. The wear depth was calculated based on Archard's Law, using contact pressures, micromotions, and a wear factor, which was determined from accelerated fretting experiments. Results. The forces for the comfortable walking speed had the highest peak forces for the maximum peak subject, with a maximum peak force of 3644 N, followed by walking up stairs, with a similar maximum peak force of 3626 N. The chair rise had a lower maximum peak force of 2240 N (−38.5%). The simulated volumetric wear followed the trends seen in the peaks of the predicted hip joint forces, with the largest wear volumes predicted for a comfortable walking speed, followed by the stairs up activity and the chair rise (Figure 3). The subjects with the highest peak forces produced the most volumetric wear in all cases. However, the lowest peak subject had a higher volumetric wear for the stairs up case than the average subject. Discussion. This study explored the effect of subject-specific variations in hip joint loads on taper fretting. The results indicate that taper wear was predominantly affected by the magnitudes of the peak forces, rather than by the orientation of the force. A more comprehensive study, capturing the full spectrum of patient variability, can help identifying parameters that accelerate fretting corrosion. Such a study should also incorporate other sources of variability, including surgical factors such as implant orientation, sizing, and offset. These factors also affect hip joint forces, and can be evaluated in musculoskeletal models such as presented here

Occlusal loading and muscle forces during mastication aids in assessment of dental restorations and implants and jaw implant design; however, three-dimensional bite forces cannot be measured with conventional transducers, which obstruct the native occlusion. The aim of this study was to combine accurate jaw kinematics measurements, together with subject-specific computational modelling, to estimate subject-specific occlusal loading and muscle forces during mastication. Motion experiments were performed on one male participant (age: 39yrs, weight: 82kg) with healthy dentition. Two low-profile magnetic sensors were fixed to the participant's teeth and the two dental arches digitised using an intra-oral scanner. The participant performed ten continuous of chewing on a polyurethane rubber sample of known material properties, followed by maximal compression (clenching). This was repeated at the molars, premolars of both the left and right sides, and central incisors. Jaw motion was simultaneously recorded from the sensors, and finite element modelling used to estimate bite force. Specifically, simulations of chewing and biting were performed by driving the model using the measured kinematics, and bite force magnitude and direction quantified. Muscle forces were then evaluated using a rigid-body musculoskeletal model of the patient's jaw. The first molars generated the largest bite forces during chewing (left: 309 N, right: 311 N) and maximum-force biting (left: 496 N, right: 495 N). The incisors generated the smallest bite forces during chewing (75 N) and maximum-force biting (114 N). The anterior temporalis and superficial masseter muscles had the largest contribution to maximum bite force, followed by the posterior temporalis and medial pterygoid muscles. This study presents a new method for estimating dynamic occlusal loading and muscle forces during mastication. These techniques provide new knowledge of jaw biomechanics, including muscle and occlusal loading, which will be useful in surgical planning and jaw implant design

The aim of this research was to determine biomechanical markers which differentiate medial knee osteoarthritis (OA) patients who do and do not show structural progression over a 2-year period. A cohort of 36 subjects was selected from a longitudinal study (Meireles et al 2017) using Kellgren-Lawrence (KL) scores at baseline and 2-year follow-up. The cohort consisted of 10 healthy controls (HC) (KL=0 at both time points), 15 medial knee OA non-progressors (NPKOA) (KL≥1 at baseline and no change over 2 years), and 11 medial knee OA progressors (PKOA) (KL≥1 at baseline and increase of ≥1 over 2 years). 3D integrated motion capture data from three walking trials were processed through a musculoskeletal modelling framework (Smith et al 2016) to estimate knee joint loading parameters (i.e., magnitude of mean contact pressure, and centre of pressure (COP)). Parameters at first and second peak were extracted and compared between groups using Kruskal-Wallis and Mann-Whitney tests. Higher magnitudes were observed in PKOA vs NPKOA, and PKOA vs HC groups at both time points. Additionally, a posterior (1st and 2nd peak), and lateral (2nd peak) shift in medial compartment COP was shown between PKOA and NPKOA, and PKOA and HC subjects. Interestingly, in the studied parameters, no differences were observed between NPKOA and HC groups. Significantly higher magnitude, and a more posterior and lateral COP was observed between PKOA and NPKOA patients. These differences, combined with an absence of difference between NPKOA and HC suggest structural OA progression is driven by a combination of altered loading magnitude and location. These results may serve as guidelines for targeted gait retraining rehabilitation to slow or stop knee OA progression whereby shifting COP anterior and medial and reducing magnitude by ~22% may shift patients from a PKOA to a NPKOA trajectory

Altered mechanical loading is a widely suggested, but poorly understood potential cause of cartilage degeneration in osteoarthritis. In rodents, osteoarthritis is induced following destabilization of the medial meniscus (DMM). This study estimates knee kinematics and contact forces in rats with DMM to gain better insight into the specific mechanisms underlying disease development in this widely-used model. Unilateral knee surgery was performed in adult male Sprague-Dawley rats (n=5 with DMM, n=5 with sham surgery). Radio-opaque beads were implanted on their femur and tibia. 8 weeks following knee surgery, rat gait was recorded using the 3D²YMOX setup (Sanctorum et al. 2019, simultaneous acquisition of biplanar XRay videos and ground reaction forces). 10 trials (1 per rat) were calibrated and processed in XMALab (Knörlein et al. 2016). Hindlimb bony landmarks were labeled on the XRay videos using transfer learning (Deeplabcut, Mathis et al. 2019; Laurence-Chasen et al. 2020). A generic OpenSim musculoskeletal model of the rat hindlimb (Johnson et al. 2008) was adapted to include a 3-degree-of-freedom knee. Inverse kinematics, inverse dynamics, static optimization of muscle forces, and joint reaction analysis were performed. In rats with DMM, knee adduction was lower compared to sham surgery. Ground reaction forces were less variable with DMM, resulting in less variability in joint external moments. The mediolateral ground reaction force was lower, resulting in lower hip adduction moment, thus less force was produced by the rectus femoris. Rats with DMM tended to break rather than propel, resulting in lower hip flexion moment, thus less force was produced by the semimembranosus. These results are consistent with lower knee contact forces in the anteroposterior and axial directions. These preliminary data indicate no overloading of the knee joint in rats with DMM, compared with sham surgery. We are currently expanding our workflow to finite element analysis, to examine mechanical cues in the cartilage of these rats (Fig1G)

Musculoskeletal models of the lower limb lend insight into muscle forces and joint mechanics during dynamic activities. However, traditional musculoskeletal modeling is based on rigid body assumptions, and frequently represents the knee as a hinge joint, neglecting the complex interactions between the patella, femur, and tibia. Implementation of the musculoskeletal modeling framework in an explicit finite element environment allows joint contact to be easily incorporated, as well as representation of any structure as rigid or fully deformable in order to evaluate, for example, implant stresses or bone strain. Prediction of these values is particularly valuable when evaluating implant mechanics after total knee replacement. A finite element, musculoskeletal model of an implanted right lower limb was constructed, including thirteen muscles crossing the knee joint. A Hill-type muscle model was developed to allow muscle activation within the explicit FE framework. Muscle forces were predicted by optimization of muscle activation patterns during flexion-extension and chair-rise activities. The effect of muscle path representation was investigated using two approaches: lines of action directly between the origin and insertion sites of the muscles, and lines of action along the centroid of the muscle bodies. Incorporating anatomic muscle paths into the model reduced the predicted peak quadriceps force during the chair-rise activity by 46%, and reduced the peak tibio-femoral contact pressure by 14%. In addition, bone strain was predicted during the activity for the implanted patella, and showed peak bone strain at the edge of the implant near the inferior pole. The muscle-activated models demonstrated the advantages of an explicit finite element framework, and allow rapid, rigid body simulation in addition to the full contact, deformable analyses when greater resolution is required

This study aimed to quantify self-reported outcomes and walking gait biomechanics in patients following primary and revision THA. The specific goals of this study were to investigate: (i) if primary and revision THA patients have comparable preoperative outcomes; and (2) if revision THA patients have worse postoperative outcomes than primary THA patients. Forty-three patients undergoing primary THA for osteoarthritis and 23 patients undergoing revision THA were recruited and followed longitudinally for their first 12 postoperative months. Reasons for revision were loosening (73%), dislocation (9%), and infection (18%). Patients completed the Hip dysfunction and Osteoarthritis Outcome Score (HOOS), and underwent gait analysis preoperatively, and at 3 and 12 months postoperatively. A 10 camera motion analysis system (V5 Vantage, Vicon, UK) recorded marker trajectories (100 Hz) during walking at self- selected speeds. A generic lower-body musculoskeletal model (Gait2392) was scaled using principal component analysis [1] and the inverse kinematics tool in Opensim 3.3 was used to compute joint angles for the lower limbs in the sagittal plane. Independent samples t-test were used to compare patient reported outcomes between the primary and revision groups at each timepoint. Statistical parametric mapping was used to compare gait patterns between the two groups at each timepoint. Preoperatively, patients undergoing primary THA reported significantly worse pain (p<0.001), symptoms (p<0.001), function (p<0.001), and quality of life (p=0.004). No differences were observed at 3 and 12 months postoperatively between patients who had received a primary or revision THA. The only observed difference in gait pattern was that patients with a revision THA had reduced hip extension at 3 months, but no differences were observed preoperatively and 12 months. Despite the suggestions in the literature that revision THA is bound to have worse outcomes compared to primary THA, we found no differences in in patient-reported outcomes and gait patterns at 12 months postoperatively. This suggests that it may be possible, in some circumstances, for patients following revision THA to achieve similar outcomes to their peers undergoing primary THA

Introduction. Total knee replacement (TKR) is an established and effective surgical procedure in case of advanced osteoarthritis. However, the rate of satisfied patients amounts only to about 75 %. One common cause for unsatisfied patients is the anterior knee pain, which is partially caused by an increase in patellofemoral contact force and abnormal patellar kinematics. Since the malpositioning of the tibial and the femoral component affects the interplay in the patellofemoral joint and therefore contributes to anterior knee pain, we conducted a computational study on a cruciate-retaining (CR) TKR and analysed the effect of isolated femoral and tibial component malalignments on patellofemoral dynamics during a squat motion. Methods. To analyse different implant configurations, a musculoskeletal multibody model was implemented in the software Simpack V9.7 (Simpack AG, Gilching, Germany) from the SimTK data set (Fregly et al.). The musculoskeletal model comprised relevant ligaments with nonlinear force-strain relation according to Wismans and Hill-type muscles spanning the lower extremity. The experimental data were obtained from one male subject, who received an instrumented CR TKR. Muscle forces were calculated using a variant of the computed muscle control algorithm. To enable roll-glide kinematics, both tibio- and patellofemoral joint compartments were modelled with six degrees of freedom by implementing a polygon-contact-model representing the detailed implant surfaces. Tibiofemoral contact forces were predicted and validated using data from experimental squat trials (SimTK). The validated simulation model has been used as reference configuration corresponding to the optimal surgical technique. In the following, implant configurations, i.e. numerous combinations of relative femoral and tibial component alignment were analysed: malposition of the femoral/tibial component in mediolateral (±3 mm) and anterior-posterior (±3 mm) direction. Results. Mediolateral translation/malposition of the tibial component did not show high influence on the maximal patellofemoral contact force. Regarding the mediolateral translation of the femoral component, similar tendencies were observed. However, lateralisation of the femoral component (3 mm) clearly increased the lateral patella shift and medialisation of the tibial component (3 mm) led to a slightly increased lateral patella shift. Compared to the reference model, pronounced posterior translation of the tibial and femoral component resulted in a lower patellofemoral contact force, further increasing with higher anterior translation of the components. The translation of the tibial component showed smaller influence on the patellofemoral contact force than the translation of the femoral component. Discussion. In our present study, the mediolateral malposition of the femoral and tibial component showed no major impact on patellofemoral contact force and contribution to anterior knee pain in patients with CR TKR. However, the influence of implant component positioning in anterior-posterior direction on patellofemoral contact force is evident, especially for the femoral component. Our generated musculoskeletal model can contribute to computer-assisted preclinical testing of TKR and may support clinical decision-making in preoperative planning

The tendency towards using inertial sensors for remote monitoring of the patients at home is increasing. One of the most important characteristics of the sensors is sampling rate. Higher sampling rate results in higher resolution of the sampled signal and lower amount of noise. However, higher sampling frequency comes with a cost. The main aim of our study was to determine the validity of measurements performed by low sampling frequency (12.5 Hz) accelerometers (SENS) in patients with knee osteoarthritis compared to standard sensor-based motion capture system (Xsens). We also determined the test-retest reliability of SENS accelerometers. Participants were patients with unilateral knee osteoarthritis. Gait analysis was performed simultaneously by using Xsens and SENS sensors during two repetitions of over-ground walking at a self-selected speed. Gait data from Xsens were used as an input for AnyBody musculoskeletal modeling software to measure the accelerations at the exact location of two defined virtual sensors in the model (VirtualSENS). After preprocessing, the signals from SENS and VirtualSENS were compared in different coordinate axes in time and frequency domains. ICC for SENS data from first and second trials were calculated to assess the repeatability of the measurements. We included 32 patients (18 females) with median age 70.1[48.1 – 85.4]. Mean height and weight of the patients were 173.2 ± 9.6 cm and 84.2 ± 14.7 kg respectively. The correlation between accelerations in time domain measured by SENS and VirtualSENS in different axes was r = 0.94 in y-axis (anteroposterior), r = 0.91 in x-axis (vertical), r = 0.83 in z-axis (mediolateral), and r = 0.89 for the magnitude vector. In frequency domain, the value and the power of fundamental frequencies (F. 0. ) of SENS and VirtualSENS signals demonstrated strong correlation (r = 0.98 and r = 0.99 respectively). The result of test-retest evaluation showed excellent repeatability for acceleration measurement by SENS sensors. ICC was between 0.89 to 0.94 for different coordinate axes. Low sampling frequency accelerometers can provide valid and reliable measurements especially for home monitoring of the patients, in which handling big data and sensors cost and battery lifetime are among important issues

After total hip replacement, force generating capacity of gluteal muscles is an impotant parameter on joint contact forces and primary fixation of total hip replacement. Femoral offset is an option to optimize muscle moment arms, especially main abductor Gluteus Medius and Minimus. To investigate relationship with weak gluteal muscles (Gluteus Medius and Minimus) and increased femoral offset, we build a musculoskeletal model. Creating of three-dimensional femur geometry and scaling of the musculoskeletal model according to the subject were performed with computed tomography data. Obtained gait kinematic and kinetic data were applied and to mimic gluteal muscle weakness, the force generating capacities of Gluteus Medius and Minimus reduced (%20-%80). Analysis were done for both anatomical and +10mm offset. Then, muscle and joint reaction forces obtained from musculoskeletal analysis transfered to CT based finite element model to evaluate changes in maximum principle stresses on femur. According to the results of the musculoskeletal analysis, the weakness of the gluteal muscles caused an increase in the activation of Gluteus Maximus, Rectus Femoris and Tensor Fasciae Latae. Effects of +10 mm femoral offset on total abductor muscle activity increased with reduced muscle strength. As a result of the finite element analysis, no significant difference was observed for maximum principle stresses on femur with varying muscle activites. The results of these analyses are important to understand weakness of gluteal muscles and for planning hip surgery

Objectives. The aim of this study was to investigate the biomechanical effect of the anterolateral ligament (ALL), anterior cruciate ligament (ACL), or both ALL and ACL on kinematics under dynamic loading conditions using dynamic simulation subject-specific knee models. Methods. Five subject-specific musculoskeletal models were validated with computationally predicted muscle activation, electromyography data, and previous experimental data to analyze effects of the ALL and ACL on knee kinematics under gait and squat loading conditions. Results. Anterior translation (AT) significantly increased with deficiency of the ACL, ALL, or both structures under gait cycle loading. Internal rotation (IR) significantly increased with deficiency of both the ACL and ALL under gait and squat loading conditions. However, the deficiency of ALL was not significant in the increase of AT, but it was significant in the increase of IR under the squat loading condition. Conclusion. The results of this study confirm that the ALL is an important lateral knee structure for knee joint stability. The ALL is a secondary stabilizer relative to the ACL under simulated gait and squat loading conditions. Cite this article: Bone Joint Res 2019;8:509–517

Summary. The quantification of T1Rho relaxation times is not related with internal loading. Improvements in modeling and imaging techniques might lead to better understanding of the pathomechanics of the knee. Introduction. The onset and progression of knee osteoarthritis has been associated with an increased external knee adduction moment (EKAM). However, this external measure has no direct relationship with internal loading of the knee. For a better understanding of the pathomechanics of the knee musculoskeletal models could be used to relate external and internal knee loading. Consequently, high internal loading might cause cartilage degeneration in patients with OA. T1RhoMRI can detect changes in proteoglycan content and is therefore a non-invasive measure of cartilage degeneration in knee OA. The purpose of this study was to relate internal loading of the knee simulated by musculoskeletal models with cartilage health using T1rhoMRI. Patients & Methods. Preliminary results showed data of seven women (50–65yrs), four healthy and three OA. Subjects underwent 3D gait analysis (VICON Nexus) at comfortable walking speed, EKAM was calculated. Simulations of multi-body musculoskeletal models were driven based on the motion capture data, in order to calculate internal medial-lateral knee forces (MLforce). Besides a T1RhoMRI scan of the knee (Phillips 3T) provided cartilage health of the midsection of the medial condyle according to Pedersen et al, 2011 [4]. Differences between healthy and OA were tested with a one sided T-test, correlations between EKAM and MLforce were calculated. Results. Anthropometrics and walking speed showed no significantly different between OA patients and healthy controls. OA patients had significant larger EKAM and MLforce (p<0.05). T1Rho values were not significantly different between the groups. EKAM was positively correlated with MLforce (R. 2. =0.91, p<0.05) in healthy subjects, no association was found in knee OA patients (R. 2. < 0.01). Discussion / Conclusion. The current study demonstrates that external loading of the knee does not predict internal loading in knee OA patients. We did not find a significant effect of knee OA on cartilage quality assessed by T1Rho MRI. However a non-significant increase was visible at the posterior region of the femoral condyle in OA patients. This elevated T1Rho relaxation is in line with expectations and could be related to an increased cartilage degeneration

Abstract. Objectives. Neonatal motor development transitions from initially spontaneous to later increasingly complex voluntary movements. A delay in transitioning may indicate cerebral palsy (CP). The general movement optimality score (GMOS) evaluates infant movement variety and is used to diagnose CP, but depends on specialized physiotherapists, is time-consuming, and is subject to inter-observer differences. We hypothesised that an objective means of quantifying movements in young infants using motion tracking data may provide a more consistent early diagnosis of CP and reduce the burden on healthcare systems. This study assessed lower limb kinematic and muscle force variances during neonatal infant kicking movements, and determined that movement variances were associated with GMOS scores, and therefore CP. Methods. Electromagnetic motion tracking data (Polhemus) was collected from neonatal infants performing kicking movements (min 50° knee extension-flexion, <2 seconds) in the supine position over 7 minutes. Tracking data from lower limb anatomical landmarks (midfoot inferior, lateral malleolus, lateral knee epicondyle, ASIS, sacrum) were applied to subject-scaled musculoskeletal models (Gait2354_simbody, OpenSim). Inverse kinematics and static optimisation were applied to estimate lower limb kinematics (knee flexion, hip flexion, hip adduction) and muscle forces (quadriceps femoris, biceps femoris) for isolated kicks. Functional principal component analysis (fPCA) was carried out to reduce kicking kinematic and muscle force waveforms to PC scores capturing ‘modes’ of variance. GMOS scores (lower scores = reduced variety of movement) were collected in parallel with motion capture by a trained operator and specialised physiotherapist. Pearson's correlations were performed to assess if the standard deviation (SD) of kinematic and muscle force waveform PC scores, representing the intra-subject variance of movement or muscle activation, were associated with the GMOS scores. Results. The study compared GMOS scores, kinematics, and muscle force variances from a total of 26 infants with a mean corrected gestational age of 39.7 (±3.34) weeks and GMOS scores between 21 and 40. There was a significant association between the SD of the PC scores for knee flexion and the GMOS scores (PC1: R = 0.59, p = 0.002; PC2: R = 0.49, p = 0.011; PC3: R = 0.56, p = 0.003). The three PCs captured variances of the overall flexion magnitude (66% variance explained), early-to-late kick knee extension (20%), and continual to biphasic kicking (6%). For hip flexion, only the SD of PC1 correlated with GMOS scores (PC1: R = 0.52, p = 0.0068), which captured the variance of the overall flexion magnitude (81%). For the biceps femoris, the SD of PC1 and PC3 associated with GMOS scores (PC1: R = 0.50, p = 0.002; PC3: R = 0.45, p = 0.03), which captured the variance of the overall bicep force magnitude (79%) and early-to-late kick bicep activation (8%). Conclusions. Infants with reduced motor development as scored in the GMOS displayed reduced variances of knee and hip flexion and biceps femoris activation across kicking cycles. These findings suggest that combining objectively measured movement variances with existing classification methods could facilitate the development of more consistent and accurate diagnostic tools for early detection of CP. Declaration of Interest. (b) declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported:I declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research project