Repeatability and precision of a weighted centroid method for estimating dynamicin vivotibiofemoral surface interactions in sheep

A Cadaveric Comparison of the Kinematic and Anatomical Axes and Arthrokinematics of the Metatarsosesamoidal and First Metatarsophalangeal Joints

Presently, developments in weightbearing computed tomography and biplanar fluoroscopy technologies offer exciting avenues for investigating normative and pathologic foot function with increasing precision. Still, data quantifying sesamoid bone and proximal phalange motion are currently sparse. To express joint kinematics and compare various clinical cohorts, future studies of first ray motion will necessitate robust coordinate frames that respect the variations in underlying anatomy while also aligning closely with the functional, physiological axes of motion. These activity-dependent functional axes may be represented by a mean helical axis of the joint motion. Our cadaveric study quantified joint kinematics from weightbearing computed tomography scans during simulated toe lift and heel rise tasks. We compared the spatial orientations of the mean finite helical axes of the metatarsosesamoidal and metatarsophalangeal joints to the primary joint axis of two relevant methods for defining metatarsal coordinate frames: inertial axes and fitting of geometric primitives. The resultant kinematics exhibited less crosstalk when using a metatarsal coordinate system based on fitting cylindrical primitives to the bony anatomy compared to using principal component axes. Respective metatarsophalangeal and metatarsosesamoidal arthrokinematic contact paths and instantaneous centers of rotation were similar between activities and agree well with currently published data. This study outlines a methodology for quantitatively assessing the efficacy and utility of various anatomical joint coordinate system definitions. Improvements in our ability to characterize the shape and motion of foot bones in the context of functional tasks will elucidate their biomechanical roles and aid clinicians in refining treatment strategies.

Does interpolation and intra-user variability affect the accuracy of arthrokinematic measurements in the knee? A dual fluoroscopic imaging and model-based tracking study

Model-based tracking (MBT) is a time-consuming and semiautomatic approach, and thus subject to errors during the tracking process. The present study aimed primarily to quantify the effects that interpolation and intra-user variability associated with MBT have on the kinematic and arthrokinematic measurements in comparison to a gold standard radiostereometric analysis (RSA). Cadaveric knee specimens were imaged at 125 Hz while simulating standing, walking, jogging, and lunging motions. (Arthro)kinematic metrics were calculated via MBT without interpolation, MBT with two interpolation techniques when every fifth or tenth frame was analyzed, and RSA. Tracking the same activity multiple times affected (p-value, largest mean difference) the flexion-extension (FE) joint angle during walking (0.03, 0.6°), and the internal-external joint angle during jogging (0.048, -0.9°). Only during jogging for the FE joint angle was there an effect of interpolation (0.046, 0.3°). Neither tracking multiple times nor interpolation affected arthrokinematic metrics (contact path locations and excursions). The present study is the first to quantify the effects that intra-user variability and interpolation have on the (arthro)kinematic measurement accuracy using MBT. Results suggest interpolation may be used without sacrificing (arthro)kinematic outcome measurement accuracy and the errors associated with intra-user variability, while small, were larger than errors due to interpolation.

Kinematic instability in the joints of flatfoot subjects during walking: A biplanar fluoroscopic study

Abnormal foot kinematics is observed in flatfoot subjects with postural foot deformity. We aimed to investigate joint instability in flatfoot subjects by analyzing the abnormal rotational position and speed of their joints while walking. Five flatfoot subjects participated in our study. Three-dimensional motions of the tibia, talus, calcaneus, navicular, and cuboid were obtained during walking using the biplanar fluoroscopic motion analyses. An anatomical coordinate system was established for each bone. The rotations and ranges of motion (ROMs) of the joints from heel-strike to toe-off were quantified. The relative movements on the articular surfaces were quantified by surface relative velocity vector analysis. The data from flat foot subjects were compared with the data from normal foot subjects in previous studies. The average relative speed on the articular surface of the tibiotalar, subtalar, and calcaneocuboid joints for the flatfoot subjects was significantly higher (p < 0.05) than that for the normal foot subjects. The flatfoot subjects exhibited increased movements toward plantar flexion in the tibiotalar joint, and eversion and external rotations in the talonavicular joint during the stance phase, compared to the normal subjects (p < 0.01). Furthermore, the flatfoot subjects had a significantly larger ROM along with the inversion/eversion rotations (5.6 ± 1.8° vs. 10.7 ± 4.0°) and internal/external rotations (7.1 ± 1.5° vs. 10.5 ± 3.5°) in the tibiotalar joint. The flatfoot subjects demonstrated abnormal kinematics and larger joint movements in multiple joints during the mid-stance and terminal stance phases of walking. This demonstrates their high instability levels.

Three-Dimensional Kinematic Coupling of the Healthy Knee During Treadmill Walking

Accurate joint kinematics plays an important role in estimating joint kinetics in musculoskeletal simulations. Biplanar fluoroscopic (BPF) systems have been introduced to measure skeletal kinematics with six degrees-of-freedom. The purpose of this study was to model knee kinematic coupling using knee kinematics during walking, as measured by the BPF system. Seven healthy individuals (mean age, 23 ± 2 yr) performed treadmill walking trials at 1.2 m/s. Knee kinematics was regressed separately for the swing and stance phases using a generalized mixed effects model. Tibial anterior translation function was y=0.20x−3.09 for the swing phase and y=0.31x−0.54 for the stance phase, where x was the flexion angle and y was the tibial anterior translation. Tibial lateral and inferior translation were also regressed separately for the stance phase and the swing phase. Tibial external rotation was y=−0.002x2+0.19x−0.64 for the swing phase and y=−0.19x−1.22 for the stance phase. The tibial adduction rotation function was also calculated separately for the stance and swing phase. The study presented three-dimensional coupled motion in the knee during the stance and swing phases of walking, and demonstrated the lateral pivoting motion found in previous studies. This expanded understanding of secondary knee motion functions will benefit musculoskeletal simulation and help improve the accuracy of calculated kinetics.

Tibio-femoral joint contact in healthy and osteoarthritic knees during quasi-static squat: A bi-planar X-ray analysis

The aim of this study was to quantify the tibio-femoral contact point (CP) locations in healthy and osteoarthritic (OA) subjects during a weight-bearing squat using stand-alone biplanar X-ray images. Ten healthy and 9 severe OA subjects performed quasi-static squats. Bi-planar X-ray images were recorded at 0°, 15°, 30°, 45°, and 70° of knee flexion. A reconstruction/registration process was used to create 3D models of tibia, fibula, and femur from bi-planar X-rays and to measure their positions at each posture. A weighted centroid of proximity algorithm was used to calculate the tibio-femoral CP locations. The accuracy of the reconstruction/registration process in measuring the quasi-static kinematics and the contact parameters was evaluated in a validation study. The quasi-static kinematics data revealed that in OA knees, adduction angles were greater (p<0.01), and the femur was located more medially relative to the tibia (p<0.01). Similarly, the average CP locations on the medial and lateral tibial plateaus of the OA patients were shifted (6.5±0.7mm; p<0.01) and (9.6±3.1mm; p<0.01) medially compared to the healthy group. From 0° to 70° flexion, CPs moved 8.1±5.3mm and 8.9±5.3mm posteriorly on the medial and lateral plateaus of healthy knees; while in OA joints CPs moved 10.1±8.4mm and 3.6±2.8mm posteriorly. The average minimum tibio-femoral bone-to-bone distances of the OA joints were lower in both compartments (p<0.01). The CPs in the OA joints were located more medially and displayed a higher ratio of medial to lateral posterior translations compared to healthy joints.

Alterations in knee contact forces and centers in stance phase of gait: A detailed lower extremity musculoskeletal model

Evaluation of contact forces-centers of the tibiofemoral joint in gait has crucial biomechanical and pathological consequences. It involves however difficulties and limitations in in vitro cadaver and in vivo imaging studies. The goal is to estimate total contact forces (CF) and location of contact centers (CC) on the medial and lateral plateaus using results computed by a validated finite element model simulating the stance phase of gait for normal as well as osteoarthritis, varus-valgus and posterior tibial slope altered subjects. Using foregoing contact results, six methods commonly used in the literature are also applied to estimate and compare locations of CC at 6 periods of stance phase (0%, 5%, 25%, 50%, 75% and 100%). TF joint contact forces are greater on the lateral plateau very early in stance and on the medial plateau thereafter during 25-100% stance periods. Large excursions in the location of CC (>17mm), especially on the medial plateau in the mediolateral direction, are computed. Various reported models estimate quite different CCs with much greater variations (~15mm) in the mediolateral direction on both plateaus. Compared to our accurately computed CCs taken as the gold standard, the centroid of contact area algorithm yielded least differences (except in the mediolateral direction on the medial plateau at ~5mm) whereas the contact point and weighted center of proximity algorithms resulted overall in greatest differences. Large movements in the location of CC should be considered when attempting to estimate TF compartmental contact forces in gait.

A new measure of tibiofemoral subchondral bone interactions that correlates with early cartilage damage in injured sheep

We have demonstrated previously that chondral damage is associated with increased knee surface velocities following ligament and meniscus injuries in sheep. We tested the hypothesis that cartilage damage scores would correlate with a new bone surface interaction measure that captures complex changes in tibiofemoral alignment, "proximity disturbance" (PD). Six sheep underwent combined anterior cruciate and medial collateral ligament transection (ACL/MCLx), five complete lateral meniscectomy (Mx), and four sham arthrotomy (Sham). Tibiofemoral subchondral bone surfaces were modeled, and the post-operative changes in relative separation of the surfaces (i.e., "proximity") were derived from subject-specific in vivo 3D stifle kinematics. Surface areas of regions of near contact were determined, and PD was calculated as the range of change in tibiofemoral proximity, divided by normalized overlapping proximity surface areas between baseline and post-operative time points. Cartilage morphology was graded at dissection. ACL/MCLx PD was significantly elevated relative to Mx and Shams, and correlated with cartilage damage (r(2)  = 0.88-0.98). Although not statistically significant, Mx PD values tended to be higher than those of Shams, and correlated with cartilage damage. Results from both injury models suggest that increasing change in tibiofemoral surface alignment may be increasingly deleterious to long-term cartilage health in sheep.