Medial and lateral knee contact forces during walking, stair ascent and stair descent are more affected by contact locations than tibiofemoral alignment in knee osteoarthritis patients with varus malalignment.

Introduction: Knee OA progression is related to medial knee contact forces, which can be altered by anatomical parameters of tibiofemoral alignment and contact point locations. There is limited and controversial literature on medial-lateral force distribution and the effect of anatomical parameters, especially in motor activities different from walking. We analyzed the effect of tibiofemoral alignment and contact point locations on knee contact forces, and the medial-lateral force distribution in knee OA subjects with varus malalignment during walking, stair ascending and stair descending. Methods: Fifty-one knee OA subjects with varus malalignment underwent weight-bearing radiographs and motion capture during walking, stair ascending and stair descending. We created a set of four musculoskeletal models per subject with increasing level of personalization, and calculated medial and lateral knee contact forces. To analyze the effect of the anatomical parameters, statistically-significant differences in knee contact forces among models were evaluated. Then, to analyze the force distribution, the medial-to-total contact force ratios were calculated from the fully-informed models. In addition, a multiple regression analysis was performed to evaluate correlations between forces and anatomical parameters. Results: The anatomical parameters significantly affected the knee contact forces. However, the contact points decreased medial forces and increased lateral forces and led to more marked variations compared to tibiofemoral alignment, which produced an opposite effect. The forces were less medially-distributed during stair negotiation, with medial-to-total ratios below 50% at force peaks. The anatomical parameters explained 30%-67% of the variability in the knee forces, where the medial contact points were the best predictors of medial contact forces. Discussion: Including personalized locations of contact points is crucial when analyzing knee contact forces in subjects with varus malalignment, and especially the medial contact points have a major effect on the forces rather than tibiofemoral alignment. Remarkably, the medial-lateral force distribution depends on the motor activity, where stair ascending and descending show increased lateral forces that lead to less medially-distributed loads compared to walking.

Differences in hip musculoskeletal loads between limbs during daily activities in patients with 3D-printed hemipelvic reconstructions following tumor surgery.

BACKGROUND: Anatomical custom-made prostheses, thanks to computer-aided design and 3D-printing technology, help improve osseointegration and reduce mechanical complications in bone reconstructions following bone tumors. A recent quantitative analysis of long-term recovery in patients with 3D-printed reconstructions following pelvic tumor surgery showed asymmetries in ground reaction forces between limbs during different motor activities, while standing very good motor performance and quality of life. RESEARCH QUESTION: We analyzed hip contact forces and muscle forces in that cohort of six patients with an innovative custom-made reconstruction of the hemipelvis, and we tested the hypothesis that asymmetries in ground reaction forces would result in more marked differences in musculoskeletal forces. METHODS: State-of-the-art musculoskeletal modeling in an optimization-based inverse-dynamics workflow was used to calculate hip contact forces and muscle forces during five motor activities, and the differences between limbs were statistically evaluated across the motor activity cycles and on the force peaks. RESULTS: The musculoskeletal loads were found to be not symmetric, as hip loads were generally higher in the contralateral limb. We found significant differences in considerable portions of the motor activities cycles except squat, load symmetry indices indicating a load increase (median up to 25%) on the contralateral limb, especially during stair descent and chair rise/sit, and significantly higher values in the contralateral limb at force peaks. SIGNIFICANCE: We confirmed the hypothesis that residual asymmetries found in ground reaction forces were amplified when hip musculoskeletal loads were investigated, reflecting a shift of the loads toward the intact limb. Despite the general trend of higher loads found in the contralateral hip, this cannot be considered a risk of overloading, as both hips supported loads in a physiological range or lower, indicating a likely optimal recovery.

Long-term functional recovery in patients with custom-made 3D-printed anatomical pelvic prostheses following bone tumor excision.

BACKGROUND: Anatomical custom-made prostheses make it possible to reconstruct complicated bone defects following excision of bone tumors, thanks to 3D-printed technology. To date, clinical measures have been used to report clinical-functional outcome and provide evidence for the effectiveness of this new surgical approach. However, there are no studies that quantified the achievable recovery during common activities by using instrumental clinical-functional evaluation in these patients. RESEARCH QUESTION: What is the motor performance, functional outcome and quality of life in patients with custom-made 3D-printed pelvic prostheses following bone tumor? METHODS: To analyze motor performance, six patients performed motion analysis during five motor activities at follow-up of 32 ± 18 months. Joint angles, ground reaction forces and joint moments of the operated and contralateral limbs were compared. On-off activity of lower-limb muscles were calculated from electromyography and compared to a healthy matched population. To analyze functional outcome and quality of life, differences in measured hip abductor strength between limbs were evaluated, as well as clinical-functional scores (Harris Hip Score, Barthel Index, Musculoskeletal Tumor Society score), and quality of life (SF-36 health survey). RESULTS: We found only slight differences in joint kinematics when comparing operated and contralateral limb. The activity of gluteal muscles was normal, while hamstrings showed out-of-phase activities. Ground reaction forces and hip moments showed asymmetries between limbs, particularly in more demanding motor activities. We found a mean difference in hip abductor strength of 48 ± 82 N between limbs, good clinical-functional scores, and quality of life scores within normative. SIGNIFICANCE: Our study showed optimal long-term results in functional recovery, mainly achieved through recovery of the gluteal function, although minor impairments were found, which may be considered for future improvement of this innovative surgery. The effect of a more loaded contralateral limb on internal loads and long-term performance of the implant remains unknown and deserves further investigation.

Abductor muscle strengthening in THA patients operated with minimally-invasive anterolateral approach for developmental hip dysplasia.

OBJECTIVE: In developmental hip dysplasia (DDH) patients, the chronic dislocation of the femoral head with respect to the true acetabulum determines muscle contracture and atrophy, particularly of the abductor muscles, and leads to secondary osteoarthritis (OA) with severe motor dysfunction, pain and disability. The correct positioning when a total hip replacement (THR) is performed is fundamental to the recovery of gait function. Also, inadequate rehabilitation of the abductor muscles for pelvic stabilisation can be responsible for residual dysfunction. Consensus on a programme for abductor muscle strengthening in these patients is not currently available. The aim of this study was to evaluate the effectiveness of a specific program of exercises for strengthening the abductor muscles in these patients. METHODS: A multicentre, prospective, randomised clinical trial was carried out in an outpatient rehabilitation setting on 103 patients given THR for DDH through a minimally-invasive anterolateral approach. Patients were randomly divided into a Study Group, including 46 patients, and a Control Group, including 57 patients. All patients underwent standard early postoperative rehabilitation. In addition, the Study Group were given an extra 2-week rehabilitation once full weight-bearing on the operated limb was allowed, aimed at strengthening the abductor muscles. All patients were evaluated preoperatively, and at about 3 and 6 months postoperatively. Clinical measures (lower limb-length differences, hip range of motion, abductor muscle strength), and functional measures (WOMAC, HHS, 10mt WT, SF-12) were taken. RESULTS: Hip range of motion and functional outcome measures showed a progressive improvement at follow ups in both groups, significantly higher in the Study Group. In particular, abductor strength at 6 months post-op improved by 92.5% with respect to 35.7% in the Control Group. CONCLUSION: In addition to standard rehabilitation, a rehabilitation programme for strengthening the gluteal muscles in DDH patients who underwent THR determined an increase in muscle strength that improved functional performance and patient satisfaction.

Component positioning and ceramic damage in cementless ceramic-on-ceramic total hip arthroplasty.

BACKGROUND: In ceramic-on-ceramic (CoC) total hip arthroplasty (THA), component positioning demonstrated to influence the bearing damage: however the connection between angles and clinical outcomes at long-term follow-ups is currently lacking. Aims of this study were: the computer tomography (CT) assessment of component positioning in CoC THAs; the correlation analysis between positioning and ceramic damage; the identification of safe zones. METHODS: 91 consecutive post-operative CT scans including two types of CoC implants, with a mean follow-up of 12 ± 4.4 years, were evaluated. III generation (74.2%) and IV generation (25.8%) CoC surfaces were included. The angle measurements (cup abduction, anteversion, cup tilt, stem antetorsion, sacral slope) were automated using a CT-based software. The combined anteversion was assessed as well as the cup-neck position at -15°, 0°, 45° and 90° of flexion. Ceramic damage was diagnosed using synovial fluid analyses and radiological criteria. RESULTS: 63.7% of THAs was inside the cup abduction target 30°-45° and 68.1% was inside the cup anteversion target 5°-25°. 19 patients (20.9%) showed signs of ceramic damage. High cup abduction and high cup-neck 45° minimum angle (which stood for high abduction and extreme combined version) significantly correlated with ceramic damage. No demographical features apart III generation ceramic bearings influenced the results. No safe zones could be detected. CONCLUSIONS: In CoC THA, no safe zones can be described. However it is important to avoid cup inclination over 45° and a combination of steep cup and extreme combined version.

nmsBuilder: Freeware to create subject-specific musculoskeletal models for OpenSim.

BACKGROUND AND OBJECTIVE: Musculoskeletal modeling and simulations of movement have been increasingly used in orthopedic and neurological scenarios, with increased attention to subject-specific applications. In general, musculoskeletal modeling applications have been facilitated by the development of dedicated software tools; however, subject-specific studies have been limited also by time-consuming modeling workflows and high skilled expertise required. In addition, no reference tools exist to standardize the process of musculoskeletal model creation and make it more efficient. Here we present a freely available software application, nmsBuilder 2.0, to create musculoskeletal models in the file format of OpenSim, a widely-used open-source platform for musculoskeletal modeling and simulation. nmsBuilder 2.0 is the result of a major refactoring of a previous implementation that moved a first step toward an efficient workflow for subject-specific model creation. METHODS: nmsBuilder includes a graphical user interface that provides access to all functionalities, based on a framework for computer-aided medicine written in C++. The operations implemented can be used in a workflow to create OpenSim musculoskeletal models from 3D surfaces. A first step includes data processing to create supporting objects necessary to create models, e.g. surfaces, anatomical landmarks, reference systems; and a second step includes the creation of OpenSim objects, e.g. bodies, joints, muscles, and the corresponding model. RESULTS: We present a case study using nmsBuilder 2.0: the creation of an MRI-based musculoskeletal model of the lower limb. The model included four rigid bodies, five degrees of freedom and 43 musculotendon actuators, and was created from 3D surfaces of the segmented images of a healthy subject through the modeling workflow implemented in the software application. CONCLUSIONS: We have presented nmsBuilder 2.0 for the creation of musculoskeletal OpenSim models from image-based data, and made it freely available via nmsbuilder.org. This application provides an efficient workflow for model creation and helps standardize the process. We hope this would help promote personalized applications in musculoskeletal biomechanics, including larger sample size studies, and might also represent a basis for future developments for specific applications.

Bone adaptation of a biologically reconstructed femur after Ewing sarcoma: Long-term morphological and densitometric evolution.

Combining bone allografts and vascularized fibular autografts in intercalary reconstructions after resection of bone sarcomas is of particular interest in young patients as it facilitates bone healing and union and helps reduce fractures. However, adverse events related to bone adaptation still occur. Bone adaptation is driven by mechanical loading, but no quantitative biomechanical studies exist that would help surgical planning and rehabilitation. We analyzed the bone adaptation of a successful femoral reconstruction after Ewing sarcoma during 76-month follow-up using a novel methodology that allows CT-based quantification of morphology and density. The results indicated that the vital allograft promoted bone adaptation in the reconstruction. However, an overall negative balance of bone remodeling and a progressive mineral density decrease in the femoral neck might threaten long-term bone safety. These concerns seem related to both surgical technique and mechanical stimuli, where a stiff metal implant may determine load sharing, which negatively affects bone remodeling.

Relationship between bone adaptation and in-vivo mechanical stimulus in biological reconstructions after bone tumor: A biomechanical modeling analysis.

BACKGROUND: Biomechanical interpretations of bone adaptation in biological reconstructions following bone tumors would be crucial for orthopedic oncologists, particularly if based on quantitative observations. This would help plan for surgical treatments, rehabilitative programs and communication with the patients. We aimed to analyze the biomechanical adaptation of a femoral reconstruction after Ewing sarcoma according to an increasingly-used surgical technique, and to relate in-progress bone resorption to the mechanical stimulus induced by different motor activities. METHODS: We created a multiscale musculoskeletal and finite element model from CT scans and motion analysis data at a 76-month follow-up of a patient, to analyze muscle and joint loads, and to compare the mechanical competence of the reconstructed bone with the contralateral limb, in the current real condition and in a possible revision surgery that removed proximal screws. FINDINGS: Our results showed strategies of muscle coordination that led to differences in joint loads between limbs more marked in more demanding motor activities, and generally larger in the contralateral limb. The operated femur presented a markedly low ratio of physiological strain due to load-sharing with the metal implant, particularly in the lateral aspect. The possible revision surgery would help restore a physiological strain configuration, while the safety of the reconstruction would not be threatened. INTERPRETATION: We suggest that bone resorption is related to load-sharing and to the internal forces exerted during movement, and the mechanical stimulus should be improved by adopting modifications in the surgical treatment and by promoting physical therapy aimed at specific muscle strengthening.

Effect of lower-limb joint models on subject-specific musculoskeletal models and simulations of daily motor activities.

Understanding the validity of using musculoskeletal models is critical, making important to assess how model parameters affect predictions. In particular, assumptions on joint models can affect predictions from simulations of movement, and the identification of image-based joints is unavoidably affected by uncertainty that can decrease the benefits of increasing model complexity. We evaluated the effect of different lower-limb joint models on muscle and joint contact forces during four motor tasks, and assessed the sensitivity to the uncertainties in the identification of anatomical four-bar-linkage joints. Three MRI-based musculoskeletal models having different knee and ankle joint models were created and used for the purpose. Model predictions were compared against a baseline model including simpler and widely-adopted joints. In addition, a probabilistic analysis was performed by perturbing four-bar-linkage joint parameters according to their uncertainty. The differences between models depended on the motor task analyzed, and there could be marked differences at peak loading (up to 2.40 BW at the knee and 1.54 BW at the ankle), although they were rather small over the motor task cycles (up to 0.59 BW at the knee and 0.31 BW at the ankle). The model including more degrees of freedom showed more discrepancies in predicted muscle activations compared to measured muscle activity. Further, including image-based four-bar-linkages was robust to simulate walking, chair rise and stair ascent, but not stair descent (peak standard deviation of 2.66 BW), suggesting that joint model complexity should be set according to the imaging dataset available and the intended application, performing sensitivity analyses.

Sensitivity of predicted muscle forces during gait to anatomical variability in musculotendon geometry.

Scaled generic musculoskeletal models are commonly used to drive dynamic simulations of motions. It is however, acknowledged that not accounting for variability in musculoskeletal geometry and musculotendon parameters may confound the simulation results, even when analysing control subjects. This study documents the three-dimensional anatomical variability of musculotendon origins and insertions of 33 lower limb muscles determined based on magnetic resonance imaging in six subjects. This anatomical variability was compared to the musculotendon point location in a generic musculoskeletal model. Furthermore, the sensitivity of muscle forces during gait, calculated using static optimization, to perturbations of the musculotendon point location was analyzed with a generic model. More specific, a probabilistic approach was used: for each analyzed musculotendon point, the three-dimensional location was re-sampled with a uniform Latin hypercube method within the anatomical variability and the static optimization problem was then re-solved for all perturbations. We found that musculotendon point locations in the generic model showed only variable correspondences with the anatomical variability. The anatomical variability of musculotendon point location did affect the calculated muscle forces: muscles most sensitive to perturbations within the anatomical variability are iliacus and psoas. Perturbation of the gluteus medius anterior, iliacus and psoas induces the largest concomitant changes in muscle forces of the unperturbed muscles. Therefore, when creating subject-specific musculoskeletal models, these attachment points should be defined accurately. In addition, the size of the anatomical variability of the musculotendon point location was not related to the sensitivity of the calculated muscle forces.

Sensitivity of a subject-specific musculoskeletal model to the uncertainties on the joint axes location.

Subject-specific musculoskeletal models have become key tools in the clinical decision-making process. However, the sensitivity of the calculated solution to the unavoidable errors committed while deriving the model parameters from the available information is not fully understood. The aim of this study was to calculate the sensitivity of all the kinematics and kinetics variables to the inter-examiner uncertainty in the identification of the lower limb joint models. The study was based on the computer tomography of the entire lower-limb from a single donor and the motion capture from a body-matched volunteer. The hip, the knee and the ankle joint models were defined following the International Society of Biomechanics recommendations. Using a software interface, five expert anatomists identified on the donor's images the necessary bony locations five times with a three-day time interval. A detailed subject-specific musculoskeletal model was taken from an earlier study, and re-formulated to define the joint axes by inputting the necessary bony locations. Gait simulations were run using OpenSim within a Monte Carlo stochastic scheme, where the locations of the bony landmarks were varied randomly according to the estimated distributions. Trends for the joint angles, moments, and the muscle and joint forces did not substantially change after parameter perturbations. The highest variations were as follows: (a) 11° calculated for the hip rotation angle, (b) 1% BW × H calculated for the knee moment and (c) 0.33 BW calculated for the ankle plantarflexor muscles and the ankle joint forces. In conclusion, the identification of the joint axes from clinical images is a robust procedure for human movement modelling and simulation.

Are subject-specific musculoskeletal models robust to the uncertainties in parameter identification?

Subject-specific musculoskeletal modeling can be applied to study musculoskeletal disorders, allowing inclusion of personalized anatomy and properties. Independent of the tools used for model creation, there are unavoidable uncertainties associated with parameter identification, whose effect on model predictions is still not fully understood. The aim of the present study was to analyze the sensitivity of subject-specific model predictions (i.e., joint angles, joint moments, muscle and joint contact forces) during walking to the uncertainties in the identification of body landmark positions, maximum muscle tension and musculotendon geometry. To this aim, we created an MRI-based musculoskeletal model of the lower limbs, defined as a 7-segment, 10-degree-of-freedom articulated linkage, actuated by 84 musculotendon units. We then performed a Monte-Carlo probabilistic analysis perturbing model parameters according to their uncertainty, and solving a typical inverse dynamics and static optimization problem using 500 models that included the different sets of perturbed variable values. Model creation and gait simulations were performed by using freely available software that we developed to standardize the process of model creation, integrate with OpenSim and create probabilistic simulations of movement. The uncertainties in input variables had a moderate effect on model predictions, as muscle and joint contact forces showed maximum standard deviation of 0.3 times body-weight and maximum range of 2.1 times body-weight. In addition, the output variables significantly correlated with few input variables (up to 7 out of 312) across the gait cycle, including the geometry definition of larger muscles and the maximum muscle tension in limited gait portions. Although we found subject-specific models not markedly sensitive to parameter identification, researchers should be aware of the model precision in relation to the intended application. In fact, force predictions could be affected by an uncertainty in the same order of magnitude of its value, although this condition has low probability to occur.

Influence of weak hip abductor muscles on joint contact forces during normal walking: probabilistic modeling analysis.

The weakness of hip abductor muscles is related to lower-limb joint osteoarthritis, and joint overloading may increase the risk for disease progression. The relationship between muscle strength, structural joint deterioration and joint loading makes the latter an important parameter in the study of onset and follow-up of the disease. Since the relationship between hip abductor weakness and joint loading still remains an open question, the purpose of this study was to adopt a probabilistic modeling approach to give insights into how the weakness of hip abductor muscles, in the extent to which normal gait could be unaltered, affects ipsilateral joint contact forces. A generic musculoskeletal model was scaled to each healthy subject included in the study, and the maximum force-generating capacity of each hip abductor muscle in the model was perturbed to evaluate how all physiologically possible configurations of hip abductor weakness affected the joint contact forces during walking. In general, the muscular system was able to compensate for abductor weakness. The reduced force-generating capacity of the abductor muscles affected joint contact forces to a mild extent, with 50th percentile mean differences up to 0.5 BW (maximum 1.7 BW). There were greater increases in the peak knee joint loads than in loads at the hip or ankle. Gluteus medius, particularly the anterior compartment, was the abductor muscle with the most influence on hip and knee loads. Further studies should assess if these increases in joint loading may affect initiation and progression of osteoarthritis.

Muscle discretization affects the loading transferred to bones in lower-limb musculoskeletal models.

Modelling the mechanical effect of muscles is important in several research and clinical contexts. However, few studies have investigated the effect of different muscle discretizations from a mechanical standpoint. The present study evaluated the errors of a reduced discretization of the lower-limb muscles in reproducing the muscle loading transferred to bones. Skeletal geometries and a muscle data collection were derived from clinical images and dissection studies of two cadaver specimens. The guidelines of a general method previously proposed for a different anatomical district were followed. The data collection was used to calculate the mechanical effect of muscles, i.e. the generalized force vectors, and the errors between a large and a reduced discretization, in a reference skeletal pose and in the extreme poses of the range of motion of joints. The results showed that the errors committed using a reduced representation of muscles could be significant and higher than those reported for a different anatomical region. In particular, the calculated errors were found to be dependent on the individual anatomy and on the skeletal pose. Since different biomechanical applications may require different discretization levels, care is suggested in identifying the number of muscle lines of action to be used in musculoskeletal models.

Femoral loads during gait in a patient with massive skeletal reconstruction.

BACKGROUND: Biological massive skeletal reconstructions in tumours adopt a long rehabilitation protocol aimed at minimising the fracture risk. To improve rehabilitation and surgical procedures it is important to fully understand the biomechanics of the reconstructed limb. The aim of the present study was to develop a subject-specific musculoskeletal model of a patient with a massive biological skeletal reconstruction, to investigate the loads acting on the femur during gait, once the rehabilitation protocol was completed. METHODS: A personalised musculoskeletal model of the patient's lower limbs was built from a CT exam and registered with the kinematics recorded in a gait analysis session. Predicted activations for major muscles were compared to EMG signals to assess the model's predictive accuracy. FINDINGS: Gait kinematics showed only minor discrepancies between the two legs and was compatible with normality data. External moments showed slightly higher differences and were almost always lower on the operated leg exhibiting a lower variability. In the beginning of the stance phase, the joint moments were, conversely, higher on the operated side and showed a higher variability. This pattern was reflected and amplified on the femoral forces where the differences became important: on the hip, a maximum difference of 1.6 BW was predicted. The variability of the forces seemed, generally, lower on the operated leg than on the contralateral one. INTERPRETATION: Small asymmetries in kinematic patterns might be associated, in massive skeletal reconstruction, to significant difference in the skeletal loads (up to 1.6 BW for the hip joint reaction) during gait.