Analysis of the Relationship between Peak Stress and Proteoglycan Loss following Injurious Compression of Human Post-mortem Knee and Ankle Cartilage

While traumatic joint injuries are known to increase the risk of osteoarthritis (OA), the mechanism is not known. Models for injurious compression of cartilage may identify predictors of injury that suggest a clinical mechanism. We investigated the relationship between peak stress during compression and glycosaminoglycan (GAG) loss after injury for knee and ankle cartilages. Human cartilage explant disks were harvested post-mortem from the knee and ankle of three organ donors with no history of OA and subjected to injurious compression to 65% strain in uniaxial unconfined compression at 2 mm/s (400%/s). The GAG content of the conditioned medium was measured 3 days after injury. After injury of knee cartilage disks, damage was visible in 18 of 39 disks (36%). Three days after injury, the increase in GAG loss to the medium (GAG loss from injured disks minus GAG loss from location-matched uncompressed controls) was 1.5±0.3 μg/disk (mean ± SEM). With final strain and compression velocity held constant, we observed that increasing peak stress during injury was associated with less GAG loss after injury (P<0.001). In contrast, ankle cartilage appeared damaged after injury in only 1 of 16 disks (6%), there was no increase in GAG loss (0.0±0.3 μg/disk), and no relationship between peak stress and increase in GAG loss was detected (P=0.51). By itself, increasing peak stress did not appear to be an important cause of GAG loss from human cartilage in our injurious compression model. However, we observed further evidence for differences in the response of knee and ankle cartilages to injury.     

A comparison of patellofemoral cartilage morphology and deformation in anterior cruciate ligament deficient versus uninjured knees

Anterior cruciate ligament (ACL) deficient patients have an increased rate of patellofemoral joint (PFJ) osteoarthritis (OA) as compared to the general population. Although the cause of post-injury OA is multi-factorial, alterations in joint biomechanics may predispose patients to cartilage degeneration. This study aimed to compare in vivo PFJ morphology and mechanics between ACL deficient and intact knees in subjects with unilateral ACL ruptures. Eight male subjects underwent baseline MRI scans of both knees. They then performed a series of 60 single-legged hops, followed by a post-exercise MRI scan. This process was repeated for the contralateral knee. The MR images were converted into three-dimensional surface models of cartilage and bone in order to assess cartilage thickness distributions and strain following exercise. Prior to exercise, patellar cartilage was significantly thicker in intact knees as compared to ACL deficient knees by 1.8%. In response to exercise, we observed average patellar cartilage strains of 5.4 ± 1.1% and 2.5 ± 1.4% in the ACL deficient and intact knees, respectively. Importantly, the magnitude of patellar cartilage strain in the ACL deficient knees was significantly higher than in the intact knees. However, while trochlear cartilage experienced a mean strain of 2.4 ± 1.6%, there was no difference in trochlear cartilage strain between the ACL deficient and uninjured knees. In summary, we found that ACL deficiency was associated with decreased patellar cartilage thickness and increased exercise-induced patellar cartilage strain when compared to the uninjured contralateral knees.     

Influence of kinematic analysis methods on detecting ankle and subtalar joint instability

Patients with subtalar joint instability may be misdiagnosed with ankle instability, which may lead to chronic instability at the subtalar joint. Therefore, it is important to understand the difference in kinematics after ligament sectioning and differentiate the changes in kinematics between ankle and subtalar instability. Three methods may be used to determine the joint kinematics; the Euler angles, the Joint Coordinate System (JCS) and the helical axis (HA). The purpose of this study was to investigate the influence of using either method to detect subtalar and ankle joints instability. 3D kinematics at the ankle and subtalar joint were analyzed on 8 cadaveric specimens while the foot was intact and after sequentially sectioning the anterior talofibular ligament (ATFL), the calcaneofibular ligament (CFL), the cervical ligament and the interosseous talocalcaneal ligament (ITCL). Comparison in kinematics calculated from sensor and anatomical landmarks was conducted as well as the influence of Euler angles and JCS rotation sequence (between ISB recommendation and previous research) on the subtalar joint. All data showed a significant increase in inversion when the ITCL was sectioned. There were differences in the data calculated using sensors coordinate systems vs. anatomic coordinate systems. Anatomic coordinate systems were recommended for these calculations. The Euler angle and JCS gave similar results. Differences in Euler angles and JCS sequence lead to the same conclusion in detecting instability at the ankle and subtalar joint. As expected, the HA detected instability in plantarflexion at the ankle joint and in inversion at the subtalar joint.     

Determination of dynamic ankle ligament strains from a computational model driven by motion analysis based kinematic data

External rotation of the foot has been implicated in high ankle sprains. Recent studies by this laboratory, and others, have suggested that torsional traction characteristics of the shoe-surface interface may play a role in ankle injury. While ankle injuries most often involve damage to ligaments due to excessive strains, the studies conducted by this laboratory and others have largely used surrogate models of the lower extremity to determine shoe-surface interface characteristics based on torque measures alone. The objective of this study was to develop a methodology that would integrate a motion analysis-based kinematic foot model with a computational model of the ankle to determine dynamic ankle ligament strains during external foot rotation. Six subjects performed single-legged, internal rotation of the body with a planted foot while a marker-based motion analysis was conducted to track the hindfoot motion relative to the tibia. These kinematic data were used to drive an established computational ankle model. Ankle ligament strains, as a function of time, were determined. The anterior tibiofibular ligament (ATiFL) experienced the highest strain at 9.2±1.1%, followed by the anterior deltoid ligament (ADL) at 7.8±0.7%, averaged over the six subjects. The peak ATiFL strain occurred prior to peak strain in the ADL in all subjects. This novel methodology may provide new insights into mechanisms of high ankle sprains and offer a basis for future evaluations of shoe-surface interface characteristics using human subjects rather than mechanical surrogate devices.     

Effect of simulated joint instability and bracing on ankle and subtalar joint flexibility

It is clinically challenging to distinguish between ankle and subtalar joints instability in vivo. Understanding the changes in load-displacement at the ankle and subtalar joints after ligament injuries may detect specific changes in joint characteristics that cannot be detected by investigating changes in range of motion alone. The effect of restricting joints end range of motion with ankle braces was already established, but little is known about the effect of an ankle brace on the flexibility of the injured ankle and subtalar joints. Therefore, the purposes of this study were to (1) understand how flexibility is affected at the ankle and subtalar joints after sectioning lateral and intrinsic ligaments during combined sagittal foot position and inversion and during internal rotation and (2) investigate the effect of a semi-rigid ankle brace on the ankle and subtalar joint flexibility. Kinematics and kinetics were collected from nine cadaver feet during inversion through the range of ankle flexion and during internal rotation. Motion was applied with and without a brace on an intact foot and after sequentially sectioning the calcaneofibular ligament (CFL) and the intrinsic ligaments. Segmental flexibility was defined as the slope of the angle-moment curve for each 1 Nm interval. Early flexibility significantly increased at the ankle and subtalar joint after CFL sectioning during inversion. The semi-rigid ankle brace significantly decreased early flexibility at the subtalar joint during inversion and internal rotation for all ligament conditions and at the ankle joint after all ligaments were cut.     

Effect of ACL deficiency on MCL strains and joint kinematics

The knee joint is partially stabilized by the interaction of multiple ligament structures. This study tested the interdependent functions of the anterior cruciate ligament (ACL) and the medial collateral ligament (MCL) by evaluating the effects of ACL deficiency on local MCL strain while simultaneously measuring joint kinematics under specific loading scenarios. A structural testing machine applied anterior translation and valgus rotation (limits 100 N and 10 N m, respectively) to the tibia of ten human cadaveric knees with the ACL intact or severed. A three-dimensional motion analysis system measured joint kinematics and MCL tissue strain in 18 regions of the superficial MCL. ACL deficiency significantly increased MCL strains by 1.8% (p<0.05) during anterior translation, bringing ligament fibers to strain levels characteristic of microtrauma. In contrast, ACL transection had no effect on MCL strains during valgus rotation (increase of only 0.1%). Therefore, isolated valgus rotation in the ACL-deficient knee was nondetrimental to the MCL. The ACL was also found to promote internal tibial rotation during anterior translation, which in turn decreased strains near the femoral insertion of the MCL. These data advance the basic structure-function understanding of the MCL, and may benefit the treatment of ACL injuries by improving the knowledge of ACL function and clarifying motions that are potentially harmful to secondary stabilizers.     

Tibiotalar cartilage stress corresponds to T2 mapping: application to barefoot running in novice and marathon-experienced runners

Abstract

A 3D anatomically-based finite-element foot model was adopted for predicting von Mises stresses within tibiotalar cartilage following 5?km barefoot running. To compare this predicted stress with T2 maps, magnetic resonance scans of the right ankle and plantar pressure were obtained from ten novices and ten marathon-experienced runners before and after running. Following running, tibiotalar cartilage stress was decreased in experienced runners. This corresponded with T2 values that did not change between pre- and post-running suggesting no increase in cartilage fluid levels. In contrast, novices maintained the same level of von Mises stress and this corresponded with a significant T2 increase in tibiotalar cartilage.     

Determination of real-time in-vivo cartilage contact deformation in the ankle joint

The knowledge of real-time in-vivo cartilage deformation is important for understanding of cartilage function and biomechanical factors that may relate to cartilage degeneration. This study investigated cartilage contact area and peak contact compressive strain of four healthy human ankle joints as a function of time using a combined magnetic resonance (MR) and dual-orthogonal fluoroscopic imaging technique. Each ankle was subjected to a different constant loading (between 700 and 820 N). The cartilage contact deformation was obtained from the first second to 300 s after the load was applied. In all ankle joints studied in this paper, contact strains increased to 24-38% at first 20 s after loading. Beyond 20 s, the change of cartilage contact deformation was relatively small and varied in a rate close to zero beyond 50 s. These data indicated that the cartilage contact areas and contact strain could raise dramatically right after loading and reach a relatively stable condition within 1 min after constant loading. The history of cartilage deformation determined in this study may provide a real-time boundary condition for 3D finite element simulation of in vivo cartilage contact stress in the joint as a function of time.     

A computational algorithm to simulate disorganization of collagen network in injured articular cartilage

Cartilage defects are a known risk factor for osteoarthritis. Estimation of structural changes in these defects could help us to identify high risk defects and thus to identify patients that are susceptible for the onset and progression of osteoarthritis. Here, we present an algorithm combined with computational modeling to simulate the disorganization of collagen fibril network in injured cartilage. Several potential triggers for collagen disorganization were tested in the algorithm following the assumption that disorganization is dependent on the mechanical stimulus of the tissue. We found that tensile tissue stimulus alone was unable to preserve collagen architecture in intact cartilage as collagen network reoriented throughout the cartilage thickness. However, when collagen reorientation was based on both tensile tissue stimulus and tensile collagen fibril strains or stresses, the collagen network architecture was preserved in intact cartilage. Using the same approach, substantial collagen reorientation was predicted locally near the cartilage defect and particularly at the cartilage–bone interface. The developed algorithm was able to predict similar structural findings reported in the literature that are associated with experimentally observed remodeling in articular cartilage. The proposed algorithm, if further validated, could help to predict structural changes in articular cartilage following post-traumatic injury potentially advancing to impaired cartilage function.     

A real-time computational model for estimating kinematics of ankle ligaments

Background: An accurate assessment of ankle ligament kinematics is crucial in understanding the injury mechanisms and can help to improve the treatment of an injured ankle, especially when used in conjunction with robot-assisted therapy. A number of computational models have been developed and validated for assessing the kinematics of ankle ligaments. However, few of them can do real-time assessment to allow for an input into robotic rehabilitation programs. Method: An ankle computational model was proposed and validated to quantify the kinematics of ankle ligaments as the foot moves in real-time. This model consists of three bone segments with three rotational degrees of freedom (DOFs) and 12 ankle ligaments. This model uses inputs for three position variables that can be measured from sensors in many ankle robotic devices that detect postures within the foot–ankle environment and outputs the kinematics of ankle ligaments. Validation of this model in terms of ligament length and strain was conducted by comparing it with published data on cadaver anatomy and magnetic resonance imaging. Results: The model based on ligament lengths and strains is in concurrence with those from the published studies but is sensitive to ligament attachment positions. Conclusions: This ankle computational model has the potential to be used in robot-assisted therapy for real-time assessment of ligament kinematics. The results provide information regarding the quantification of kinematics associated with ankle ligaments related to the disability level and can be used for optimizing the robotic training trajectory.     

Full field strain measurements of collagenous tissue by tracking fiber alignment through vector correlation

Full field strain measurements of biological tissue during loading are often limited to the quantification of fiduciary marker displacements on the tissue surface. These marker measurements can lack the necessary spatial resolution to characterize non-uniform deformation and may not represent the deformation of the load-bearing collagen microstructure. To overcome these potential limitations, a method was developed to track the deformation of the collagen fiber microstructure in ligament tissue. Using quantitative polarized light imaging, fiber alignment maps incorporating both direction and alignment strength at each pixel were generated during facet capsular ligament loading. A grid of virtual markers was superimposed over the tissue in the alignment maps, and the maximization of a vector correlation calculation between fiber alignment maps was used to track marker displacement. Tracking error was quantified through comparisons to the displacements of excised ligament tissue (n=3); separate studies applied uniaxial tension to isolated facet capsular ligament tissue (n=4) to evaluate tracking capabilities during large tissue deformations. The average difference between virtual marker and tissue displacements was 0.07±0.06 pixels. This error in marker location produced principal strain measurements of 1.2±1.6% when markers were spaced 4 pixels apart. During tensile tissue loading, substantial inhomogeneity was detected in the strain field using vector correlation tracking, and the location of maximum strain differed from that produced by standard tracking techniques using coarser meshes. These findings provide a method to directly measure fiber network strains using quantitative fiber alignment data, enabling a better understanding of structure–function relationships in tissues at different length scales.     

Can altered neuromuscular coordination restore soft tissue loading patterns in anterior cruciate ligament and menisci deficient knees during walking?

Injuries to the anterior cruciate ligament (ACL) and menisci commonly lead to early onset osteoarthritis. Treatments that can restore normative cartilage loading patterns may mitigate the risk of osteoarthritis, though it is unclear whether such a goal is achievable through conservative rehabilitation. We used musculoskeletal simulation to predict cartilage and ligament loading patterns during walking in intact, ACL deficient, menisci deficient, and ACL-menisci deficient knees. Stochastic simulations with varying coordination strategies were then used to test whether neuromuscular control could be modulated to restore normative knee mechanics in the pathologic conditions. During early stance, a 3 mm increase in anterior tibial translation was predicted in the ACL deficient knee. Mean cartilage contact pressure increased by 18% and 24% on the medial and lateral plateaus, respectively, in the menisci deficient knee. Variations in neuromuscular coordination were insufficient to restore normative cartilage contact patterns in either the ACL or menisci deficient knees. Elevated cartilage contact pressures in the pathologic knees were observed in regions where cartilage wear patterns have previously been reported. These results suggest that altered cartilage tissue loading during gait may contribute to region-specific degeneration patterns, and that varying neuromuscular coordination in isolation is unlikely to restore normative knee mechanics.     

Adsorption and Distribution of Fluorescent Solutes near the Articular Surface of Mechanically Injured Cartilage

The development of cartilage-specific imaging agents supports the improvement of tissue assessment by minimally invasive means. Techniques for highlighting cartilage surface damage in clinical images could provide for sensitive indications of posttraumatic injury and early stage osteoarthritis. Previous studies in our laboratory have demonstrated that fluorescent solutes interact with cartilage surfaces strongly enough to affect measurement of their partition coefficients within the tissue bulk. In this study, these findings were extended by examining solute adsorption and distribution near the articular surface of mechanically injured cartilage. Using viable cartilage explants injured by an established protocol, solute distributions near the articular surface of three commonly used fluorophores (fluorescein isothiocyanate (FITC), tetramethylrhodamine isothiocyanate (TRITC), and carboxytetramethylrhodamine (TAMRA)) were observed after absorption and subsequent desorption to assess solute-specific matrix interactions and reversibility. Both absorption and desorption processes demonstrated a trend of significantly less solute adsorption at surfaces of fissures compared to adjacent intact surfaces of damaged explants or surfaces of uninjured explants. After adsorption, normalized mean surface intensities of fissured surfaces of injured explants were 6%, 40%, and 32% for FITC, TRITC, and TAMRA, respectively, compared to uninjured surfaces. Similar values were found for sliced explants and after a desorption process. After desorption, a trend of increased solute adsorption at the site of intact damaged surfaces was noted (316% and 238% for injured and sliced explants exposed to FITC). Surface adsorption of solute was strongest for FITC and weakest for TAMRA; no solutes negatively affected cell viability. Results support the development of imaging agents that highlight distinct differences between fissured and intact cartilage surfaces.     

T2 relaxation time measurements in tibiotalar cartilage after barefoot running and its relationship to ankle biomechanics

The influence of ankle kinematics and plantar pressure from mid-range barefoot running on T2 relaxation times of tibiotalar cartilage is unknown. This study aimed to quantitatively evaluate the T2 relaxation time of tibiotalar cartilage and ankle biomechanics following 5 km barefoot running. Twenty healthy runners (who had no 5 km barefoot running experience) underwent 3.0-Tesla magnetic resonance (MR) scans and assessment of running gait before and after 5 km barefoot running. Participants were divided into two groups consisting of marathon-experienced (n = 10) and novice (n = 10) with equal number of males and females in each group. Three musculoskeletal radiologists measured T2 relaxation times in 18 regions of the ankle cartilage: anterior zone, central zone, and posterior zone, or lateral, middle, and medial sections in the sagittal plane. Three-dimensional ankle kinetics, kinematics, and plantar pressure were all also assessed during barefoot running. In the novice group, the T2 relaxation time in the posterior zone of tibial cartilage (p = 0.001) and lateral section in both tibial (p = 0.02) and talar (p = 0.02) cartilage were significantly increased after barefoot running. Ankle kinematics exhibited significant changes in females. Plantar loading was shifted from the medial to lateral aspect after running. This included a significant reduction in the loading under the toes and the 1st, 2nd and 3rd metatarsals, with a significant increase under the 4th and 5th metatarsals and lateral midfoot. The results suggest that plantar pressure may directly lead to local increases in cartilage T2 signal, which was not associated with changes in ankle kinematics.     

3D characterization of bone strains in the rat tibia loading model

Bone strain is considered one of the factors inducing bone tissue response to loading. Nevertheless, where animal studies can provide detailed data on bone response, they only offer limited information on experimental bone strains. Including micro-CT-based finite element (micro FE) models in the analysis represents a potent methodology for quantifying strains in bone. Therefore, the main objective of this study was to develop and validate specimen-specific micro FE models for the assessment of bone strains in the rat tibia compression model. Eight rat limbs were subjected to axial compression loading; strain at the medio-proximal site of the tibiae was measured by means of strain gauges. Specimen-specific micro FE models were created and analyzed. Repeated measurements on each limb indicated that the effect of limb positioning was small (COV?= 6.45 ± 2.27 %). Instead, the difference in the measured strains between the animals was high (54.2%). The computational strains calculated at the strain gauge site highly correlated to the measured strains (R ²?=?0.95). Maximum peak strains calculated at exactly 25% of the tibia length for all specimens were equal to 435.11 ± 77.88 microstrains (COV?=?17.19%). In conclusion, we showed that strain gauge measurements are very sensitive to the exact strain gauge location on the bone; hence, the use of strain gauge data only is not recommended for studies that address at identifying reliable relationships between tissue response and local strains. Instead, specimen-specific micro FE models of rat tibiae provide accurate estimates of tissue-level strains.     

Ligament fibre recruitment and forces for the anterior drawer test at the human ankle joint

Although the anterior drawer test at the ankle joint is commonly used in routine clinical practice, very little is known about the sharing of load between the individual passive structures and the joint response at different flexion angles.A mathematical model of the ankle joint was devised to calculate ligament fibre recruitment and load/displacement curves at different flexion angles. Ligaments were modelled as three-dimensional arrays of fibres, and their orientations at different flexion angles were taken from a previously validated four-bar-linkage model in the sagittal plane. A non-linear stress/strain relationship was assumed for ligament fibres and relevant mechanical parameters were taken from two reports in the literature. Talus and calcaneus were assumed to move as a single rigid body. Antero/distal motion of the talus relative to the tibia was analysed.The ankle joint was found to be stiffer at the two extremes of the flexion range, and the highest laxity was found around the neutral position, confirming previous experimental works. With a first dataset, a 20N anterior force produced 4.3, 5.5, and 4.4mm displacement respectively at 20 degrees plantarflexion, at neutral, and at 20 degrees dorsiflexion. At 10 degrees plantarflexion, for a 6mm displacement, 65% of the external force was supported by the anterior talofibular, 11% by the deep anterior tibiotalar and 5.5% by the tibionavicular ligament. Corresponding results from a second dataset were 1.4, 2.4 and 1.8mm at 40N force, and 80%, 0% and 2% for a 3mm displacement. A component of the contact force supported the remainder.     

Partial medial meniscectomy and rotational differences at the knee during walking

Loss of meniscal function due to injury or partial meniscectomy is common and represents a significant risk factor for premature osteoarthritis. The menisci can influence the transverse plane movements (anterior–posterior (AP) translation and internal–external (IE) rotation) of the knee during walking. While walking is the most frequent activity of daily living, the kinematic differences at the knee during walking associated with the meniscal injury are not well understood. This study examined the influence of partial medial meniscectomy (PMM) on the kinematics and kinetics of the knee during the stance phase of gait by testing the differences in anterior–posterior translation, internal–external rotation, knee flexion range of movement, peak flexion/extension moments, and adduction moments between the PMM and healthy contralateral limbs. Ten patients (45±9 years old, height 1.75±0.06 m, weight 76.7±13.5 kg) who had undergone partial medial meniscectomy (33±100 months post-op) in one limb with a healthy contralateral limb were tested during normal walking. The contralateral limb was compared to a matched control group and no differences were found. The primary kinematic difference was a significantly greater external rotation (3.2°) of the tibia that existed through stance phase, with 8 of 10 subjects demonstrating the same pattern. The PMM subjects also exhibited significantly lower peak flexion and extension moments in their PMM limbs. The altered rotational position found likely results in changes of tibio-femoral contact during walking and could cause the type of degenerative changes found in the articular cartilage following meniscal injury.     

Viscoelastic and failure properties of spine ligament collagen fascicles

The microstructural volume fractions, orientations, and interactions among components vary widely for different ligament types. If these variations are understood, however, it is conceivable to develop a general ligament model that is based on microstructural properties. This paper presents a part of a much larger effort needed to develop such a model. Viscoelastic and failure properties of porcine posterior longitudinal ligament (PLL) collagen fascicles were determined. A series of subfailure and failure tests were performed at fast and slow strain rates on isolated collagen fascicles from porcine lumbar spine PLLs. A finite strain quasi-linear viscoelastic model was used to fit the fascicle experimental data. There was a significant strain rate effect in fascicle failure strain (P < 0.05), but not in failure force or failure stress. The corresponding average fast-rate and slow-rate failure strains were 0.098 ± 0.062 and 0.209 ± 0.081. The average failure force for combined fast and slow rates was 2.25 ± 1.17 N. The viscoelastic and failure properties in this paper were used to develop a microstructural ligament failure model that will be published in a subsequent paper.     

Rotational stiffness of football shoes influences talus motion during external rotation of the foot

Shoe-surface interface characteristics have been implicated in the high incidence of ankle injuries suffered by athletes. Yet, the differences in rotational stiffness among shoes may also influence injury risk. It was hypothesized that shoes with different rotational stiffness will generate different patterns of ankle ligament strain. Four football shoe designs were tested and compared in terms of rotational stiffness. Twelve (six pairs) male cadaveric lower extremity limbs were externally rotated 30 deg using two selected football shoe designs, i.e., a flexible shoe and a rigid shoe. Motion capture was performed to track the movement of the talus with a reflective marker array screwed into the bone. A computational ankle model was utilized to input talus motions for the estimation of ankle ligament strains. At 30 deg of rotation, the rigid shoe generated higher ankle joint torque at 46.2?±?9.3 Nm than the flexible shoe at 35.4?±?5.7 Nm. While talus rotation was greater in the rigid shoe (15.9?±?1.6 deg versus 12.1?±?1.0 deg), the flexible shoe generated more talus eversion (5.6?±?1.5 deg versus 1.2± 0.8 deg). While these talus motions resulted in the same level of anterior deltoid ligament strain (approxiamtely 5%) between shoes, there was a significant increase of anterior tibiofibular ligament strain (4.5± 0.4% versus 2.3?±?0.3%) for the flexible versus more rigid shoe design. The flexible shoe may provide less restraint to the subtalar and transverse tarsal joints, resulting in more eversion but less axial rotation of the talus during foot∕shoe rotation. The increase of strain in the anterior tibiofibular ligament may have been largely due to the increased level of talus eversion documented for the flexible shoe. There may be a direct correlation of ankle joint torque with axial talus rotation, and an inverse relationship between torque and talus eversion. The study may provide some insight into relationships between shoe design and ankle ligament strain patterns. In future studies, these data may be useful in characterizing shoe design parameters and balancing potential ankle injury risks with player performance.