Stefan van Drongelen
Mariska Wesseling2†
Jana Holder
Felix Stief- 1Dr. Rolf M. Schwiete Research Unit for Osteoarthritis, Orthopaedic University Hospital Friedrichsheim gGmbH, Frankfurt, Germany
- 2Human Movement Biomechanics Research Group, Department of Movement Sciences, KU Leuven, Leuven, Belgium
- 3Orthopaedic University Hospital Friedrichsheim gGmbH, Frankfurt, Germany
- 4Goethe University Frankfurt, Frankfurt, Germany
Reduced external knee adduction moments in the second half of stance after total hip replacement have been reported in hip osteoarthritis patients. This reduction is thought to shift the load from the medial to the lateral knee compartment and as such increase the risk for knee osteoarthritis. The knee adduction moment is a surrogate for the load distribution between the medial and lateral compartments of the knee and not a valid measure for the tibiofemoral contact forces which are the result of externally applied forces and muscle forces. The purpose of this study was to investigate whether the distribution of the tibiofemoral contact forces over the knee compartments in unilateral hip osteoarthritis patients 1 year after receiving a primary total hip replacement differs from healthy controls. Musculoskeletal modeling on gait was performed in OpenSim using the detailed knee model of Lerner et al. (2015) for 19 patients as well as for 15 healthy controls of similar age. Knee adduction moments were calculated by the inverse dynamics analysis, medial and lateral tibiofemoral contact forces with the joint reaction force analysis. Moments and contact forces of patients and controls were compared using Statistical Parametric Mapping two-sample t-tests. Knee adduction moments and medial tibiofemoral contact forces of both the ipsi- and contralateral leg were not significantly different compared to healthy controls. The contralateral leg showed 14% higher medial tibiofemoral contact forces compared to the ipsilateral (operated) leg during the second half of stance. During the first half of stance, the lateral tibiofemoral contact force of the contralateral leg was 39% lower and the ratio 32% lower compared to healthy controls. In contrast, during the second half of stance the forces were significantly higher (39 and 26%, respectively) compared to healthy controls. The higher ratio indicates a changed distribution whereas the increased lateral tibiofemoral contact forces indicate a higher lateral knee joint loading in the contralateral leg in OA patients after total hip replacement (THR). Musculoskeletal modeling using a detailed knee model can be useful to detect differences in the load distribution between the medial and lateral knee compartment which cannot be verified with the knee adduction moment.
Introduction
Hip osteoarthritis (OA) is one of the most common degenerative diseases of the musculoskeletal system (Fuchs et al., 2013). Patients, who already had a unilateral total hip replacement (THR) for end-stage hip OA, have an increased risk of developing OA in other joints of their lower extremities (Shakoor et al., 2002).
Stief et al. (2018) showed that after THR, patients reduce the external knee adduction moment (KAM) in the second half of stance of both legs compared to healthy controls whereas Shakoor et al. (2003) found a significantly higher KAM in the contralateral knee. The reduction in KAM is thought to shift the load from the medial to the lateral knee compartment (Schmidt et al., 2017). An increased probability of lateral knee OA would support the findings of Weidow et al. (2005), who showed that lateral OA of the knee was more frequently associated with hip OA on the ipsilateral side than medial knee OA.
KAM is considered to be only a surrogate for the load distribution between the medial and lateral compartments of the knee and not a valid measure for the tibiofemoral contact forces (CF) which are the result of externally applied forces and muscle forces (Winby et al., 2009). Joint CF, calculated by musculoskeletal modeling, have been validated by direct (in vivo) measurements in patients with instrumented implants (Fregly et al., 2012). An increased KAM is related to higher medial tibiofemoral contact forces (MKCF) during normal gait (Kutzner et al., 2013). However, the relation is not always straightforward as various gait retraining programs, with the intention to reduce KAM, found similar MKCF as in normal gait (Walter et al., 2010; Pizzolato et al., 2017). In addition, Pizzolato et al. (2017) reported that reduced MKCF are not inevitably associated with increased lateral tibiofemoral contact forces (LKCF).
Most studies investigating joint CF of patients receiving THR and the progression of these forces after THR, reported on hip CF (Meyer et al., 2018; O’Connor et al., 2018), whereas only few studies also examined knee CF in these patients (Shakoor et al., 2003; Wesseling et al., 2018). Shakoor et al. (2003) reported increased reaction forces (not including the muscle activation), in the medial knee compartment of the contralateral leg compared to the operated leg. When muscle forces are not included in the calculation of joint load, this may lead to an underestimation of the actual load in the joint (Lloyd and Buchanan, 2001). Several studies already showed a contribution of muscular co-contraction to higher joint CF (Trepczynski et al., 2018; Hoang et al., 2019). Wesseling et al. (2018) reported that 12 months after THR surgery the total CF on the ipsilateral knee were still lower compared to healthy controls. Since Wesseling et al. (2018) focused on the total knee CF, they could not provide information on the load on the medial and lateral compartments. It is known that after THR, hip OA patients have a higher risk for the development of OA in the contralateral knee joint (Jungmann et al., 2015; Joseph et al., 2016), however, these studies did not specify the medial or lateral knee compartment. Information on a possible asymmetry or imbalance across the medial and lateral compartments is important with regard to the development of OA: a non-physiological knee loading can precipitate OA in a healthy joint (Andriacchi et al., 2004, 2009).
In the present study the MKCF and LKCF are studied independently, using a more detailed knee model introduced by Lerner et al. (2015). Consequently, the purpose of the present study was to investigate tibiofemoral CF in the knees of unilateral hip THR patients after THR surgery separately for the medial and lateral compartment. This study hypothesized that MKCF and LKCF are different in patients with unilateral hip OA after THR surgery compared to healthy controls in such a way that the lateral knee compartment is loaded more and the medial knee compartment is loaded less.
Materials and Methods
Participants
Nineteen patients, scheduled for THR due to unilateral hip OA, were included in this study (Table 1). Exclusion criteria were: the inability to walk without walking aids, Body Mass Index (BMI) above 30, inflammatory arthritis, orthopedic surgeries within the past 6 months, OA in lower limb joints other than the affected hip, and previous lower extremity joint replacement. All surgeries were performed, using a lateral approach, by experienced orthopedic surgeons.
Table 1. Anthropometric data and walking speed of patients and healthy controls.
Patient gait data were compared to data of 15 healthy controls with similar age distribution (Stief et al., 2018). Healthy subjects were included if they had no history of orthopedic surgeries or chronic and neuromuscular disease. Our institution’s medical ethics committee approved the study (123/13 and 497/15) and all patients and healthy subjects gave informed consent prior to participation.
Gait Analysis
For the patients, three-dimensional gait analysis was performed in the week before surgery and on average 16 (6) months postoperatively. All patients were asked to walk barefoot at a self-selected speed in the level gait laboratory. Kinematic data were collected using 8 Vicon MX T10 cameras (VICON Motion Systems, Oxford, United Kingdom) operating at 200 Hz. Additionally, two AMTI force plates (Advanced Mechanical Technology Inc., Watertown, MA, United States) were used to synchronously collect kinetic data at 1000 Hz.
A marker set (called MA) was used which includes, in addition to the standardized Plug-in-Gait marker set (Kadaba et al., 1990), reflective markers on the medial malleolus, medial femoral condyle and greater trochanter (Stief et al., 2013). An overview of all markers can be found in Supplementary Figure 1. A static standing trial was performed to be able to scale the musculoskeletal model. During the static upright standing trial, participants stood barefoot, feet shoulder width apart with the knees fully extended. Of all performed trials, three postoperative trials were processed per subject.
For the control subjects, the three-dimensional gait analysis was performed only once. First, the healthy controls walked at a self-selected speed, to be able to compare the patients’ postoperative data (Table 1). Further, to compare preoperative patients to healthy controls, healthy controls also walked at a slow walking speed of approximately 1.0 ms–1, comparable to patients shortly before surgery (Schmidt et al., 2017). In the control subjects, only one side (the left side was chosen randomly) was used for further analysis and for comparison with the ipsi- and contralateral side of the patient group.
Musculoskeletal Modeling
Figure 1 shows the full-body musculoskeletal model with a detailed knee joint (Lerner et al., 2015) that was used to analyze the gait data in OpenSim 3.3 (Delp et al., 2007). The model is based on a previously described full body musculoskeletal model (Delp et al., 1990; Demers et al., 2014) and includes 18 body segments and 92 muscle-tendon actuators. Metatarsophalangeal and subtalar joints were fixed in anatomical neutral positions for all analyses (O’Connor et al., 2018). A fourth order zero-lag low-pass Butterworth filter with a cut-off frequency of 10 Hz was applied to the ground reaction forces, whereas a Woltring filter with a mean square error value of 10 was used to smooth the marker data (Woltring, 1991). Input for the model was created with the MOtoNMS toolbox (Mantoan et al., 2015) which processes the experimental data (3D marker positions and ground reaction forces) from C3D files. The joint centers of the hip, knee and ankle joint were calculated within the MOtoNMS toolbox (Mantoan et al., 2015). The hip joint center was calculated according to Harrington et al. (2007), while the knee and ankle joint centers were computed as the mid points between femoral condyle markers and the medial and lateral malleolus, respectively.
Figure 1. Model with the used marker set, combined with the workflow on implementing the motion capture data into OpenSim to calculate the joint contact forces. The inset graphic shows the medial and lateral compartment of the knee of the used musculoskeletal model according to Lerner et al. (2015).
The model was scaled for each subject and for the preoperative and postoperative measurements separately, using the marker positions of the static standing trial and the measured body mass. Inverse kinematics and inverse dynamics were performed to calculate joint angles and moments, respectively. The knee adduction moment around the medial condyle (MKAM) was normalized to body mass and expressed in Nmkg–1. Muscle forces were calculated using a static optimization approach (Anderson and Pandy, 2001). The objective function for static optimization was to minimize the sum of squared muscle activations. In a final step, joint reaction analysis was performed to calculate the MKCF and LKCF (Figure 1). CF were normalized to body weight. All data were time-normalized to the stance phase of gait (101 time points, from heel strike to foot-off of the same foot, detected using the ground reaction forces with a threshold of 20N).
Statistics
Statistical data analyses for the anthropometrics and gait speed were performed with SPSS (version 26, IBM Corporation, New York, NY, United States). Shapiro-Wilk tests were used to test for normal distribution of the above mentioned parameters. Unpaired Student’s t-tests were used to determine statistical differences between the anthropometrics of controls and patients in age, height, body mass and speed, whereas Mann-Whitney tests were used to determine statistical differences in BMI. Further, a χ2 test was used to compare gender distribution in the two groups.
Joint angles, MKAM, MKCF, LKCF and the ratio between the lateral and total CF were evaluated using Statistical Parametric Mapping (SPM). SPM is based on Random Field Theory (Adler and Taylor, 2007) and has been validated for 1D data (Pataky et al., 2013, 2016). All SPM analyses were implemented using the open-source spm1d code (version M.0.4.31) in MATLAB. Controls and patients were compared using SPM two-sample t-tests, data of the ipsilateral and contralateral leg were compared with SPM paired sample t-tests. A critical threshold of α = 0.05 was used. When the waveform exceeded the critical threshold, the data were considered significantly different in that part of stance. Differences were considered significant when differences were found between more than four successive time points, i.e., at least 4% of the stance phase of the gait cycle, similar to Wesseling et al. (2018).
Post hoc power calculations and the effect size Cohens’s d were determined for the main findings using Gpower (Faul et al., 2007) and according to Cohen (2013). To express the percentage differences between groups as well as between legs, and as input for the power calculations the values and standard deviations at the time of the maximum difference between the curves (within the band of significant differences from the SPM) were extracted.
Results
Participants
Significant differences between postoperative patients and healthy controls were found in body mass and BMI (Table 1). No significant differences were found in gender distribution, age, height, and walking speed.
Kinetics
Although slightly lower in the second half of stance, MKAM of the ipsilateral and the contralateral leg after THR were not significantly different compared to MKAM of healthy controls (Figure 2). No differences were found between the ipsilateral and contralateral leg.
Figure 2. Curves of the external knee adduction moment (bands represent standard deviation) on the left side and the SPM t-test results on the right side; comparing controls (norm in black) to ipsilateral side of postoperative patients (post-op ipsi in red) at the top (A,B), ipsilateral side of postoperative patients (post-op ipsi in red) to the contralateral side (post-op contra in blue) in the middle (C,D) and controls (norm in black) to the contralateral side of the postoperative patients (post-op contra in blue) at the bottom (E,F). When the SPM t-values exceed the critical threshold (dashed horizontal line), differences are significant. Where significant differences are found for more than 4 successive time points (4% of the stance phase of the gait cycle), areas are shaded gray and p-values are reported. The knee adduction moments are normalized to body mass and expressed in Nmkg–1.
The MKCF of the ipsilateral leg showed slightly lower values compared to the healthy controls in both the first and second half of stance, however, these differences were not significant (Figure 3). In contrast to the ipsilateral leg, the MKCF of the contralateral leg showed only lower values during the first half of stance, but also these differences did not reach a significant level. MKCF of the contralateral leg was higher (max 14%) compared to the ipsilateral leg in the second half of stance (p < 0.001, 67–81% of stance, power = 100%, effect size d = 3.66).
Figure 3. Curves of the tibiofemoral contact force on the medial epicondyle (bands represent standard deviation) on the left side and the SPM t-test results on the right side; comparing controls (norm in black) to ipsilateral side of postoperative patients (post-op ipsi in red) at the top (A,B), ipsilateral side of postoperative patients (post-op ipsi in red) to the contralateral side (post-op contra in blue) in the middle (C,D) and controls (norm in black) to the contralateral side of the postoperative patients (post-op contra in blue) at the bottom (E,F). When the SPM t-values exceed the critical threshold (dashed horizontal line), differences are significant. Where significant differences are found for more than 4 successive time points (4% of the stance phase of the gait cycle), areas are shaded gray and p-values are reported. Contact forces are normalized to body weight (BW).
Figure 4 shows that the LKCF of the ipsilateral leg was not significantly different compared to healthy controls. The contralateral leg did not differ from the ipsilateral leg; however, some differences were detected compared to the healthy controls. The contralateral leg of the patients showed a lower LKCF (39%) between 6 and 16% of stance (p = 0.002, power = 67%, effect size d = 0.73) in relation to the healthy controls whereas the contralateral leg of the patients showed a significantly higher LKCF (39%) between 82 and 91% of stance (p = 0.003, power = 64%, effect size d = 0.70).
Figure 4. Curves of the tibiofemoral contact force on the lateral epicondyle (bands represent standard deviation) on the left side and the SPM t-test results on the right side; comparing controls (norm in black) to ipsilateral side of postoperative patients (post-op ipsi in red) at the top (A,B), ipsilateral side of postoperative patients (post-op ipsi in red) to the contralateral side (post-op contra in blue) in the middle (C,D) and controls (norm in black) to the contralateral side of the postoperative patients (post-op contra in blue) at the bottom (E,F). When the SPM t-values exceed the critical threshold (dashed horizontal line), differences are significant. Where significant differences are found for more than 4 successive time points (4% of the stance phase of the gait cycle), areas are shaded gray and p-values are reported. Contact forces are normalized to body weight (BW).
The ratio between the CF on the lateral condyle and the total CF for the ipsilateral leg was not significantly different compared to the healthy controls (Figure 5). Also, no differences were found compared to the contralateral leg. Short significant differences were detected when comparing the contralateral leg to the healthy controls. The healthy controls showed a higher ratio (max 32%) between 6 and 13% of stance (p = 0.009, power = 61%, effect size d = 0.68) and a lower ratio (max 26%) compared to the contralateral leg of the patients between 85 and 91% of stance (p = 0.010, power = 75%, effect size d = 0.82).
Figure 5. Curves of the ratio between the lateral tibiofemoral contact force and the total tibiofemoral contact force (bands represent standard deviation) on the left side and the SPM t-test results on the right side; comparing controls (norm in black) to ipsilateral side of postoperative patients (post-op ipsi in red) at the top (A,B), ipsilateral side of postoperative patients (post-op ipsi in red) to the contralateral side (post-op contra in blue) in the middle (C,D) and controls (norm in black) to the contralateral side of the postoperative patients (post-op contra in blue) at the bottom (E,F). When the SPM t-values exceed the critical threshold (dashed horizontal line), differences are significant. Where significant differences are found for more than 4 successive time points (4% of the stance phase of the gait cycle), areas are shaded gray and p-values are reported.
Figures with the results of the SPM analyses of the preoperative data were reported in the Supplementary Material (Supplementary Figures 2–4).
Kinematics
The kinematics for the hip (flexion-extension and adduction-abduction) and the knee (flexion-extension) after THR were presented in Supplementary Material (Supplementary Figures 5–7).
For hip extension no differences were found compared to the healthy controls, but the ipsilateral leg showed significantly less extension in the second phase of stance compared to the contralateral leg (p < 0001, 54–100% of stance, Supplementary Figure 5).
Hip adduction was similar for the ipsilateral and contralateral leg, but significantly less compared to healthy controls during the first half of stance (18–40 and 16–58% of stance for the ipsilateral and contralateral leg, respectively, Supplementary Figure 6).
Knee flexion during the first half of stance equaled that of healthy controls and of the contralateral leg. The extent of knee extension in the ipsilateral leg was significantly lower compared to healthy controls (p < 0001, 61–100% of stance) and compared to the contralateral leg (p < 0001, 52–100% of stance, Supplementary Figure 7).
Discussion
The goal of this research was to study tibiofemoral CF during gait in unilateral hip OA patients after THR. It has been proposed that a reduced KAM in the second half of the stance phase shifts the distribution of the knee joint CF from medial to lateral and therefore, increases the load on the lateral compartment (Schmidt et al., 2017). In the present study we found no differences in the MKCF compared to healthy controls for both the ipsi-and contralateral leg. In the ipsilateral leg no significant shift of the tibiofemoral CF from the medial to the lateral condyle after THR was detected, as the LKCF were nearly identical to healthy controls. Only in the contralateral leg a higher lateral knee joint loading in the second half of the stance phase was found.
The MKAM reported here agree with the KAM reported in the study of Stief et al. (2018) in which KAM in the second half of the stance was still lower in patients after THR compared to healthy controls. Stief et al. reported significant differences whereas in the present study the values were only slightly lower. An explanation for these discrepancies might be that different methods were used to calculate the external joint moments. Stief et al. (2013) used an adapted lower body protocol and calculated the joint moments around the joint centers directly from the gait analysis (Davis et al., 1991). In the present study, the external joint moments were calculated with the inverse dynamics tool of OpenSim, where KAM is calculated at the medial condyle, which means it is shifted medially relative to the joint center. Further, Stief et al. compared peak values whereas in the present study, the whole stance phase was compared by SPM.
Like MKAM, MKCF of the ipsi- and contralateral leg did not deviate significantly from healthy controls. All in all, our MKCF agree with that of Wesseling et al. (2018) who studied longitudinal total knee CF after THR. Specifically, the second peak of the ipsilateral leg remained lower than controls 1 year after surgery, and they concluded that no overloading was present 1 year after THR, neither in the ipsilateral nor the contralateral knee. In contrast to Wesseling et al. (2018), who used the generic musculoskeletal model of OpenSim, we used the model by Lerner et al. (2015) which allows detailed analysis of medial and lateral tibiofemoral CF.
For the ipsilateral leg, lateralization of the load was not found, as the LKCF were similar to healthy controls. For the contralateral leg the LKCF were lower during the first half of stance and higher during the second half of stance compared to healthy controls. The latter might indicate a shift of the tibiofemoral CF to the lateral compartment. It must be noted that the absolute differences are small and only for a short period (between 82 and 91% of stance). The differences as such might not be clinically relevant.
The higher LKCF for the contralateral leg resulted in a higher ratio of the lateral to total CF for a period of 6% in the second half of stance. In general, the distribution of the forces over the medial and lateral compartments (a clear higher medial to lateral loading) agrees with previous data (Lerner et al., 2015). A finite element study of Yang et al. (2010) showed that, as in the present study (Figure 4), approximately 80% of the tibiofemoral forces are distributed to the medial knee compartment and 20% to the lateral compartment. Yang et al. (2010) were even able to distinguish leg alignments as the force distributed to the lateral compartment was 19, 21, and 22% for the varus, normal and valgus aligned knee, respectively. The higher ratio at the beginning and at the end of the stance phase is due to the external knee abduction moment at heel strike and toe-off (Yang et al., 2010).
Changes in external hip and knee joint moments and joint CF are caused by changes in the kinematics (Wesseling et al., 2015). Deviating kinematics in the ipsilateral leg can persist even 1 year after THR (Foucher et al., 2007). In the present study, it was shown that hip adduction (in the first half of stance) and knee extension (in the second half of stance) were still less in postoperative patients compared to the healthy controls. Recently, it has been shown that differences in sagittal plane gait kinematics between patients with knee OA and asymptomatic controls appear to be mainly a result of reduced walking speed (Ismailidis et al., 2020). Walking speed can also influence joint CF (O’Connor et al., 2018). O’Connor et al. (2018) reported higher hip CF in THR patients during walking at higher speeds. However, for the knee CF less information is available, especially when considering MKCF and LKCF. Lenton et al. (2018) published on the effect of speed on medial and total tibiofemoral CF, whereas no effect on the lateral CF was found. However, the walking speeds were much higher (1.53 and 1.81 ms–1) so it must be questioned if the results can be transferred to the slower walking speed of patients, as measured in the present study. The strength of the present study is that CFs were compared between patients after THR and a healthy control group with an almost identical walking speed (no significant differences). Accordingly, our results cannot be related to differences in walking speed.
In the present study a significant difference in body mass and BMI was found between healthy controls and patients. However, we normalized the joint moments and joint CF to body mass/body weight, thus eliminating the effect of BMI and body mass on our main findings. Also, Harding et al. (2016) found significant OA effects in peak tibiofemoral compression forces when the forces were normalized to body mass. In addition, they found significant BMI effects only in the case of absolute, non-normalized forces. In the present study, and in the above cited study by Harding et al. (2016) subjects with nominal BMI were included: therefore it is possible that differences in tibiofemoral CF, greater than those which can be attributed to higher body mass are only valid for grossly obese individuals (Browning and Kram, 2007; Messier et al., 2014).
Higher CF on the contralateral (non-operated) side in the second half of stance compared to the ipsilateral side confirm the results of other studies (Shakoor et al., 2003; Foucher and Wimmer, 2012; Wesseling et al., 2018), although not all studies reported significant differences between patients (ipsilateral or contralateral) and healthy controls (Wesseling et al., 2018) so that a surplus load cannot be confirmed. The present results do support studies reporting a higher risk for the development of OA in the contralateral knee joint after THR (Shakoor et al., 2002; Umeda et al., 2009; Gillam et al., 2013; Jungmann et al., 2015; Joseph et al., 2016) since patients are walking with an asymmetrical limb load. However, these studies only report more structural damage and progression of degenerative findings in the contralateral knee and do not specify the medial or lateral knee compartment. Gillam et al. (2013) stated that further studies are needed to investigate whether an increased risk of receiving an arthroplasty in the contralateral knee is only related to having a hip arthroplasty or if other factors such as pain play a role as well. One factor which can be decisive is the leg alignment. Previous studies found that a valgus malalignment increases LKCF (Holder et al., 2020) and therefore might increase the risk to suffer lateral knee OA and lateral cartilage damage (Felson et al., 2013). The missing information on the leg alignment is a limitation of the present study. With that information the musculoskeletal model could have been optimized for the patient-specific tibiofemoral alignment. Previous research by Lerner et al. (2015) used data from anterior-posterior full-leg radiographs to show that MKCF and LKCF can accurately predict in vivo measurements. In the present study, these data were not available as in the clinical standard only a pelvic overview X-Ray was made to check for loosening of the prosthesis. We expect that the leg alignment has only a marginal effect on the results as previous work showed that the implantation of a hip prosthesis only led to a slight increased varus alignment (1°) of the operated leg (van Drongelen et al., 2019).
Another limitation of the present study is the relatively small sample size. Post hoc power analysis and the determination of the effect size for the primary findings, however, revealed that the sample size of 19 patients provided 100% power for the difference in the MKCF between the ipsilateral and contralateral leg. Further, the differences in the LKCF and the ratio between the LKCF and total CF between our study groups (19 patients, 15 controls) had a power between 61 and 75%. The calculated effect sizes between 0.68 and 3.66 observed for the main findings in the current study reinforces the relevance of the results. Nevertheless, some results of the present study may be underpowered and the differences could be significant with greater sample sizes.
Conclusion
The hypothesis that the distribution of tibiofemoral CF is different in patients with unilateral hip OA patients after THR compared to healthy controls cannot be fully confirmed by this study. The CF on the medial tibiofemoral compartment for both the ipsi- and contralateral leg did not differ from healthy controls. Also, the CF on the lateral compartment of the ipsilateral leg was similar to controls. Only a higher LKCF and higher ratio between the LKCF and total CF in the contralateral leg indicates a higher lateral knee joint loading in hip OA patients after receiving a THR. Nevertheless, it must be noted that the differences found are only for a short period during the second half of stance. Though, musculoskeletal modeling can be useful to detect differences in the load distribution between the medial and lateral knee compartment, which cannot be verified with KAM. The present study cannot finally clarify whether this difference is sufficient to cause permanent damage to the lateral knee compartment in patients with unilateral hip OA after receiving a THR.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Medical Ethics Committee of the Faculty of Medicine of the Goethe University Frankfurt. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
SD, AM, and FS were involved in the design of the study. SD, MW, and JH carried out all data analysis. SD, MW, JH, AM, and FS were involved in the interpretation and discussion of the results. SD took the lead in writing the manuscript. All authors provided critical feedback and contributed to the final manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank Hanna Kaldowski, André Schmidt, and Anna Schäfer for their help with the data collection and Zoe Feja and Sandra Kmiec for their help with the data analysis.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2020.578030/full#supplementary-material
Abbreviations
BMI, body mass index; CF, contact forces; KAM, knee adduction moment; LKCF, lateral tibiofemoral contact force; MKAM, knee adduction moment around the medial condyle; MKCF, medial tibiofemoral contact force; OA, osteoarthritis; SPM, statistical parametric mapping; THR, total hip replacement.
Footnotes
References
Anderson, F. C., and Pandy, M. G. (2001). Static and dynamic optimization solutions for gait are practically equivalent. J. Biomech. 34, 153–161. doi: 10.1016/s0021-9290(00)00155-x
Andriacchi, T. P., Koo, S., and Scanlan, S. F. (2009). Gait mechanics influence healthy cartilage morphology and osteoarthritis of the knee. J. Bone Joint Surg. Am. 91(Suppl. 1), 95–101. doi: 10.2106/jbjs.h.01408
Andriacchi, T. P., Mündermann, A., Smith, R. L., Alexander, E. J., Dyrby, C. O., and Koo, S. (2004). A framework for the in vivo pathomechanics of osteoarthritis at the knee. Ann. Biomed. Eng. 32, 447–457. doi: 10.1023/b:abme.0000017541.82498.37
Browning, R. C., and Kram, R. (2007). Effects of Obesity on the Biomechanics of Walking at Different Speeds. Med. Sci. Sports Exerc. 39, 1632–1641. doi: 10.1249/mss.0b013e318076b54b
Davis, R. B., Õunpuu, S., Tyburski, D., and Gage, J. R. (1991). A gait analysis data collection and reduction technique. Hum. Mov. Sci. 10, 575–587. doi: 10.1016/0167-9457(91)90046-Z
Delp, S. L., Anderson, F. C., Arnold, A. S., Loan, P., Habib, A., John, C. T., et al. (2007). OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans. Biomed. Eng. 54, 1940–1950. doi: 10.1109/tbme.2007.901024
Delp, S. L., Loan, J. P., Hoy, M. G., Zajac, F. E., Topp, E. L., and Rosen, J. M. (1990). An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures. IEEE Trans. Biomed. Eng. 37, 757–767. doi: 10.1109/10.102791
Demers, M. S., Pal, S., and Delp, S. L. (2014). Changes in tibiofemoral forces due to variations in muscle activity during walking. J. Orthop. Res. 32, 769–776. doi: 10.1002/jor.22601
Faul, F., Erdfelder, E., Lang, A. G., and Buchner, A. (2007). G∗Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav. Res. Methods 39, 175–191. doi: 10.3758/BF03193146
Felson, D. T., Niu, J., Gross, K. D., Englund, M., Sharma, L., Cooke, T. D. V., et al. (2013). Valgus malalignment is a risk factor for lateral knee osteoarthritis incidence and progression: findings from MOST and the osteoarthritis initiative. Arthritis Rheum. 65, 355–362. doi: 10.1002/art.37726
Foucher, K. C., Hurwitz, D. E., and Wimmer, M. A. (2007). Preoperative gait adaptations persist one year after surgery in clinically well-functioning total hip replacement patients. J. Biomech. 40, 3432–3437. doi: 10.1016/j.jbiomech.2007.05.020
Foucher, K. C., and Wimmer, M. A. (2012). Contralateral hip and knee gait biomechanics are unchanged by total hip replacement for unilateral hip osteoarthritis. Gait Posture 35, 61–65. doi: 10.1016/j.gaitpost.2011.08.006
Fregly, B. J., Besier, T. F., Lloyd, D. G., Delp, S. L., Banks, S. A., Pandy, M. G., et al. (2012). Grand challenge competition to predict in vivo knee loads. J. Orthop. Res. 30, 503–513. doi: 10.1002/jor.22023
Fuchs, J., Rabenberg, M., and Scheidt-Nave, C. (2013). Prevalence of selected musculoskeletal conditions in Germany - results of the german health interview and examination survey for adults (DEGS1). Bundesgesundheitsblatt 56, 723–732. doi: 10.1007/s00103-013-1684-7
Gillam, M. H., Lie, S. A., Salter, A., Furnes, O., Graves, S. E., Havelin, L. I., et al. (2013). The progression of end-stage osteoarthritis: analysis of data from the Australian and Norwegian joint replacement registries using a multi-state model. Osteoarthr. Cartil. 21, 405–412. doi: 10.1016/j.joca.2012.12.008
Harding, G. T., Dunbar, M. J., Hubley-Kozey, C. L., Stanish, W. D., and Astephen Wilson, J. L. (2016). Obesity is associated with higher absolute tibiofemoral contact and muscle forces during gait with and without knee osteoarthritis. Clin. Biomech. 31, 79–86. doi: 10.1016/j.clinbiomech.2015.09.017
Harrington, M. E., Zavatsky, A. B., Lawson, S. E., Yuan, Z., and Theologis, T. N. (2007). Prediction of the hip joint centre in adults, children, and patients with cerebral palsy based on magnetic resonance imaging. J. Biomech. 40, 595–602. doi: 10.1016/j.jbiomech.2006.02.003
Hoang, H. X., Diamond, L. E., Lloyd, D. G., and Pizzolato, C. (2019). A calibrated EMG-informed neuromusculoskeletal model can appropriately account for muscle co-contraction in the estimation of hip joint contact forces in people with hip osteoarthritis. J. Biomech. 83, 134–142. doi: 10.1016/j.jbiomech.2018.11.042
Holder, J., Feja, Z., van Drongelen, S., Adolf, S., Böhm, H., Meurer, A., et al. (2020). Effect of guided growth intervention on static leg alignment and dynamic knee contact forces during gait. Gait Posture 78, 80–88. doi: 10.1016/j.gaitpost.2020.03.012
Ismailidis, P., Egloff, C., Hegglin, L., Pagenstert, G., Kernen, R., Eckardt, A., et al. (2020). Kinematic changes in patients with severe knee osteoarthritis are a result of reduced walking speed rather than disease severity. Gait Posture 79, 256–261. doi: 10.1016/j.gaitpost.2020.05.008
Joseph, G. B., Hilton, J. F., Jungmann, P. M., Lynch, J. A., Lane, N. E., Liu, F., et al. (2016). Do persons with asymmetric hip pain or radiographic hip OA have worse pain and structure outcomes in the knee opposite the more affected hip? Data from the osteoarthritis initiative. Osteoarthr. Cartil. 24, 427–435. doi: 10.1016/j.joca.2015.10.001
Jungmann, P. M., Nevitt, M. C., Baum, T., Liebl, H., Nardo, L., Liu, F., et al. (2015). Relationship of unilateral total hip arthroplasty (THA) to contralateral and ipsilateral knee joint degeneration - a longitudinal 3T MRI study from the osteoarthritis initiative (OAI). Osteoarthr. Cartil. 23, 1144–1153. doi: 10.1016/j.joca.2015.03.022
Kadaba, M. P., Ramakrishnan, H. K., and Wootten, M. E. (1990). Measurement of lower extremity kinematics during level walking. J. Orthop. Res. 8, 383–392. doi: 10.1002/jor.1100080310
Kutzner, I., Trepczynski, A., Heller, M. O., and Bergmann, G. (2013). Knee adduction moment and medial contact force–facts about their correlation during gait. PLoS One 8:e81036. doi: 10.1371/journal.pone.0081036
Lenton, G. K., Bishop, P. J., Saxby, D. J., Doyle, T. L. A., Pizzolato, C., Billing, D., et al. (2018). Tibiofemoral joint contact forces increase with load magnitude and walking speed but remain almost unchanged with different types of carried load. PLoS One 13:e0206859. doi: 10.1371/journal.pone.0206859
Lerner, Z. F., DeMers, M. S., Delp, S. L., and Browning, R. C. (2015). How tibiofemoral alignment and contact locations affect predictions of medial and lateral tibiofemoral contact forces. J. Biomech. 48, 644–650. doi: 10.1016/j.jbiomech.2014.12.049
Lloyd, D. G., and Buchanan, T. S. (2001). Strategies of muscular support of varus and valgus isometric loads at the human knee. J. Biomech. 34, 1257–1267. doi: 10.1016/s0021-9290(01)00095-1
Mantoan, A., Pizzolato, C., Sartori, M., Sawacha, Z., Cobelli, C., and Reggiani, M. (2015). MOtoNMS: a MATLAB toolbox to process motion data for neuromusculoskeletal modeling and simulation. Source Code Biol. Med. 10:12. doi: 10.1186/s13029-015-0044-4
Messier, S. P., Pater, M., Beavers, D. P., Legault, C., Loeser, R. F., Hunter, D. J., et al. (2014). Influences of alignment and obesity on knee joint loading in osteoarthritic gait. Osteoarthr. Cartil. 22, 912–917. doi: 10.1016/j.joca.2014.05.013
Meyer, C. A. G., Wesseling, M., Corten, K., Nieuwenhuys, A., Monari, D., Simon, J. P., et al. (2018). Hip movement pathomechanics of patients with hip osteoarthritis aim at reducing hip joint loading on the osteoarthritic side. Gait Posture 59, 11–17. doi: 10.1016/j.gaitpost.2017.09.020
O’Connor, J. D., Rutherford, M., Bennett, D., Hill, J. C., Beverland, D. E., Dunne, N. J., et al. (2018). Long-term hip loading in unilateral total hip replacement patients is no different between limbs or compared to healthy controls at similar walking speeds. J. Biomech. 80, 8–15. doi: 10.1016/j.jbiomech.2018.07.033
Pataky, T. C., Robinson, M. A., and Vanrenterghem, J. (2013). Vector field statistical analysis of kinematic and force trajectories. J. Biomech. 46, 2394–2401. doi: 10.1016/j.jbiomech.2013.07.031
Pataky, T. C., Robinson, M. A., and Vanrenterghem, J. (2016). Region-of-interest analyses of one-dimensional biomechanical trajectories: bridging 0D and 1D theory, augmenting statistical power. PeerJ 4:e2652. doi: 10.7717/peerj.2652
Pizzolato, C., Reggiani, M., Saxby, D. J., Ceseracciu, E., Modenese, L., and Lloyd, D. G. (2017). Biofeedback for gait retraining based on real-time estimation of tibiofemoral joint contact forces. IEEE Trans. Neural Syst. Rehabil. Eng. 25, 1612–1621. doi: 10.1109/tnsre.2017.2683488
Schmidt, A., Meurer, A., Lenarz, K., Vogt, L., Froemel, D., Lutz, F., et al. (2017). Unilateral hip osteoarthritis: the effect of compensation strategies and anatomic measurements on frontal plane joint loading. J. Orthop. Res. 35, 1764–1773. doi: 10.1002/jor.23444
Shakoor, N., Block, J. A., Shott, S., and Case, J. P. (2002). Nonrandom evolution of end-stage osteoarthritis of the lower limbs. Arthritis Rheum. 46, 3185–3189. doi: 10.1002/art.10649
Shakoor, N., Hurwitz, D. E., Block, J. A., Shott, S., and Case, J. P. (2003). Asymmetric knee loading in advanced unilateral hip osteoarthritis. Arthritis Rheum. 48, 1556–1561. doi: 10.1002/art.11034
Stief, F., Böhm, H., Michel, K., Schwirtz, A., and Döderlein, L. (2013). Reliability and accuracy in three-dimensional gait analysis: a comparison of two lower body protocols. J. Appl. Biomech. 29, 105–111. doi: 10.1123/jab.29.1.105
Stief, F., Schmidt, A., van Drongelen, S., Lenarz, K., Froemel, D., Tarhan, T., et al. (2018). Abnormal loading of the hip and knee joints in unilateral hip osteoarthritis persists two years after total hip replacement. J. Orthop. Res. 36, 2167–2177. doi: 10.1002/JOR.23886
Trepczynski, A., Kutzner, I., Schwachmeyer, V., Heller, M. O., Pfitzner, T., and Duda, G. N. (2018). Impact of antagonistic muscle co-contraction on in vivo knee contact forces. J. Neuroeng. Rehabil. 15:101. doi: 10.1186/s12984-018-0434-3
Umeda, N., Miki, H., Nishii, T., Yoshikawa, H., and Sugano, N. (2009). Progression of osteoarthritis of the knee after unilateral total hip arthroplasty: minimum 10-year follow-up study. Arch. Orthop. Trauma Surg. 129, 149–154. doi: 10.1007/s00402-008-0577-y
van Drongelen, S., Kaldowski, H., Tarhan, T., Assi, A., Meurer, A., and Stief, F. (2019). Are changes in radiological leg alignment and femoral parameters after total hip replacement responsible for joint loading during gait? BMC Musculoskelet. Disord. 20:526. doi: 10.1186/s12891-019-2832-5
Walter, J. P., D’Lima, D. D., Colwell, C. W. Jr., and Fregly, B. J. (2010). Decreased knee adduction moment does not guarantee decreased medial contact force during gait. J. Orthop. Res. 28, 1348–1354. doi: 10.1002/jor.21142
Weidow, J., Mars, I., and Karrholm, J. (2005). Medial and lateral osteoarthritis of the knee is related to variations of hip and pelvic anatomy. Osteoarthr. Cartil. 13, 471–477. doi: 10.1016/j.joca.2005.01.009
Wesseling, M., de Groote, F., Meyer, C., Corten, K., Simon, J. P., Desloovere, K., et al. (2015). Gait alterations to effectively reduce hip contact forces. J. Orthop. Res. 33, 1094–1102. doi: 10.1002/jor.22852
Wesseling, M., Meyer, C., Corten, K., Desloovere, K., and Jonkers, I. (2018). Longitudinal joint loading in patients before and up to one year after unilateral total hip arthroplasty. Gait Posture 61, 117–124. doi: 10.1016/j.gaitpost.2018.01.002
Winby, C. R., Lloyd, D. G., Besier, T. F., and Kirk, T. B. (2009). Muscle and external load contribution to knee joint contact loads during normal gait. J. Biomech. 42, 2294–2300. doi: 10.1016/j.jbiomech.2009.06.019
Woltring, H. J. (1991). Representation and calculation of 3-D joint movement. Hum. Mov. Sci. 10, 603–616. doi: 10.1016/0167-9457(91)90048-3
Keywords: joint contact forces, musculoskeletal modeling, total hip replacement, knee adduction moment, walking, osteoarthritis
Citation: van Drongelen S, Wesseling M, Holder J, Meurer A and Stief F (2020) Knee Load Distribution in Hip Osteoarthritis Patients After Total Hip Replacement. Front. Bioeng. Biotechnol. 8:578030. doi: 10.3389/fbioe.2020.578030
Received: 30 June 2020; Accepted: 28 August 2020;
Published: 21 September 2020.
Edited by:
Fabio Galbusera, Galeazzi Orthopedic Institute (IRCCS), ItalyReviewed by:
Nicola Francesco Lopomo, University of Brescia, ItalyAbdelwahed Barkaoui, International University of Rabat, Morocco
Copyright © 2020 van Drongelen, Wesseling, Holder, Meurer and Stief. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Stefan van Drongelen, cy52YW5kcm9uZ2VsZW5AZnJpZWRyaWNoc2hlaW0uZGU=; Felix Stief, Zi5zdGllZkBmcmllZHJpY2hzaGVpbS5kZQ==
†Present address: Mariska Wesseling, Materialise NV, Leuven, Belgium