Abstract Objective. Four-dimensional CT (4DCT) enables visualization of dynamic wrist joint pathology but may be compromised by motion artifacts. Deep learning (DL)-assisted motion correction has the potential to enhance wrist 4DCT quality, but model development is hindered by scarce training data. To address this, we developed a model-based motion simulation framework that generated realistic 4DCT data with paired static ground truth images to demonstrate wrist motion correction using a physics-informed DL model trained with simulated data. Approach. Wrist motion was simulated by synchronizing rotational and translational movements of carpal bones and connective tissues. Elastic deformation guided by analytically-derived motion vector fields was applied to static wrist scans at neutral positions to create artifact-free images at multiple angular positions. Virtual motion-corrupted 4DCT images were produced via forward projection and reconstruction with synchronized gantry and wrist motion. Experimental data from 4DCT cadaver scans were collected to validate simulation quality. Low-intensity region score (LIRS), positivity, and entropy were used to assess simulated motion artifacts. A physics-informed motion correction network was adapted exclusively using simulated training data and subsequently deployed on simulated, real cadaveric, and in vivo 4DCT scans. LIRS and structural similarity index (SSIM) were used to evaluate motion correction quality. Wilcoxon signed-rank tests were conducted ( α = 0.05). Results. The simulation method created realistic motion artifacts compared to real scan reference ( p > 0.05). Motion artifacts were effectively mitigated across testing cases demonstrating generalizability. In simulated scans, LIRS improved from 0.55 ± 0.12–0.75 ± 0.04 and SSIM from 0.87 ± 0.08–0.95 ± 0.03 at bone regions (each p < 0.01). In experimental cadaveric scans, LIRS improved from 0.05 ± 0.01–0.31 ± 0.17 in marker regions and from 0.58 ± 0.19–0.77 ± 0.14 in bone regions (each p < 0.01). Significance. Our simulation method facilitates training DL-assisted motion correction algorithms by providing high-quality training data, effectively supporting 4DCT image quality enhancement.
Rapid personalized computational modeling of the wrist
Abstract Background/Purpose Diagnosing wrist ligament injuries is challenging; early detection and treatment are important to prevent osteoarthritis progression. Interosseous proximity maps, a proxy measure for joint space, can be generated from volumetric imaging data and may provide important information about wrist health. Artificial intelligence (AI) could enhance accuracy of noninvasive diagnosis based on imaging-derived metrics. This work demonstrates feasibility of AI training using synthetic proximity map data generated from finite element models (FEMs). Methods Personalized wrist FEMs for two asymptomatic participants were created from four-dimensional computed tomography-derived anatomic and kinematic data. Monte Carlo sampling varied 22 ligament material properties and simulated 7,500 unique injury scenarios generating 9,000,000 labeled red, green, and blue (RGB) images of interosseous proximity vector fields from FEM-derived motions. Images were associated with 17 descriptive metrics, including gross wrist angles and bone surface pairs, and used to develop mixed-input convolutional neural networks (CNNs). Model performance was evaluated for identifying specific ligament injuries. Results Average area under receiver operating characteristic curve (AUROC) for CNNs was 0.757 across all injury types and kinematics. In a subset with clinically-relevant functional angles, the average AUROC was 0.824. Best-performing individual ligament AUROCs ranged from 0.807 to 0.999. Sensitivities and specificities exceeded 0.99 for some ligament injury simulations under specific wrist angles and bone surface pairs. Conclusion This study demonstrates the feasibility of using synthetic data from FEMs to train AI models for classifying wrist ligament injuries. Proximity-based RGB images may be a relevant biomarker of ligamentous injury.
Factors affecting trainees’ preferred timing of clinical clerkships during MD-PhD training
Abstract M.D.-Ph.D. programs in the United States have traditionally followed a “2- n -2” curricular model, with the graduate phase occurring between the pre-clerkship and clerkship portions of medical training. While well established, this format can limit trainee autonomy in shaping their physician-scientist training trajectories. In response, some programs have introduced a “3- n -1” model, allowing students to complete clerkships before beginning Ph.D. training. Our institution implemented multiple flexible curricula in 2017. Understanding why trainees choose one pathway is important as programs consider implementing more adaptable training structures. To investigate these factors, we surveyed M.D.-Ph.D. students at our institution, which offers multiple flexible curricular alternatives, about contributions to their curricular decisions. Responses supported that trainees weigh considerations across medical education, scientific development, and integrated physician-scientist training domains. Although the traditional pathway was a popular option, most students pursued one of the flexible pathways. Our findings suggest that introducing flexibility in the timing of undergraduate medical education and graduate training can support diverse educational, logistical, and personal priorities while maintaining the rigor of physician-scientist training. Offering multiple pathways empowers trainees to design trajectories that best fit their needs. Continued longitudinal studies are needed to assess the long-term impacts of flexible curricula on physician-scientist career outcomes.
Finite Element Modeling of the Scaphoid Shift Maneuver: Implications for Scapholunate Ligament injuries
Abstract Background Computational modeling is a tool being deployed for orthopaedic solutions but its use in the hand and wrist remains limited. This work used a model to simulate a clinically relevant provocative scaphoid shift maneuver (SSM) with different levels of scapholunate interosseous ligament (SLIL) injuries to observe the effect on different metrics. Methods A personalized model simulated the full SSM motion cycle from ulnar deviation with extension to radial deviation with flexion informed by the participant’s motion obtained from dynamic computed tomography. Models repeated the SSM under different levels of SLIL injury and reported changes in joint kinematics, contact mechanics, and ligament forces. Results The fully injured model increased scaphoid dorsal translation, flexion, and radial deviation compared to the intact condition and caused a subluxation of the scaphoid. Radioscaphoid contact areas were approximately 200% greater in the fully injured model compared with all others and the fully injured model was the only condition where contact force decreased across the motion cycle. Ligament forces in the intact condition were on average 33.0 N and 54.2 N for the volar and dorsal SLIL, respectively. Lastly, the long radiolunate, an extrinsic stabilizer, had forces that increased following SLIL injury. Conclusions Computational models can successfully recreate clinically observed behaviors of an SSM, including scaphoid subluxation, while providing new insights via quantification of contact mechanics and ligament forces. Contact mechanics metrics may be important for understanding the long-term progression of untreated SLIL injuries to osteoarthritis. Additionally, ligament force metrics may explain the progression of SLIL injuries from volar SLIL to dorsal SLIL and highlight the importance of repairing extrinsic stabilizers of the joint, due to increased force sharing following SLIL injury. This work provides a pathway to future studies investigating the effects of SLIL injury and repair, both acutely and chronically.
Clinical Reasoning: A 30-Year-Old Female Patient With Multiple Sclerosis Presenting With Rapidly Progressive Cranial Neuropathies, Weakness, and Ataxia
Background: Flexor pollicis longus (FPL) palsy following both-bone forearm fracture (BBFF) is a rare complication. Methods: A retrospective review of acute BBFF treated with open reduction internal fixation by a single surgeon from 2005 to 2023 was performed. Injury and surgical characteristics of patients with documented FPL palsy were reviewed. In addition, 10 cadaveric dissections were performed to evaluate the anatomy of the anterior interosseous nerve (AIN) and its branches. The distance of these branches from palpable elbow landmarks and variability in branching pattern were evaluated. Results: Twenty-nine patients underwent surgery for acute BBFF. Of these, 5 (17%) had evidence of an FPL palsy either at the time of injury presentation (n = 2) or immediately following surgery (n = 3). All patients with FPL palsy sustained fractures in the middle one-third of the radius. All palsies resolved after an average of 33 days of observation. In cadaveric dissections, the average distance from the lateral epicondyle to the AIN takeoff and branch to the FPL was 5.5 and 7.6 cm, respectively. The AIN takeoff and branch to the FPL were never less than 4 and 7 cm from the lateral epicondyle, respectively. Conclusion: Flexor pollicis longus palsy following BBFF can occur at the time of injury or following surgery. All FPL palsies involved midshaft radial fractures and were likely neurapraxia. The etiology of FPL palsy remains unclear, but cadaveric dissection suggests the FPL motor branch may be at risk from mid-to-proximal radius fracture fragments or excessive traction during surgery.
Four-Dimensional Computed Tomography-Derived Radiolunate Arthrokinematics With a Case Study in Four-Corner Arthrodesis
Background: The radiolunate joint is conventionally thought to be spared from radiocarpal osteoarthritis, but proximal lunate degenerative changes have been noted. An investigation of radiolunate arthrokinematics in ostensibly healthy wrists may improve understanding of posttraumatic and postoperative outcomes. We hypothesize that four-dimensional computed tomography (4DCT)-derived radiolunate arthrokinematics during motion will differ between wrists with type 1 and type 2 lunates. The clinical-translational relevance of this analysis is highlighted using a 4DCT case study of scaphoid excision with 4-corner arthrodesis (SEFCA). Methods: Four-dimensional computed tomography data were collected during flexion-extension and radioulnar deviation in 19 healthy wrists (mean [SD] age: 43.6 [11.7] years, 21.1% female) and 1 wrist 3 years after SEFCA (male, 52 years). Carpal osteokinematics and joint arthrokinematics were calculated to derive median radiolunate interosseous proximities at each 4DCT-captured timepoint. Results: Median radiolunate interosseous proximities during motion are well approximated by second-order regressions in uninjured wrists. In the arthrodesed wrist, closest proximities occurred at the extremes of motion except maximum radial deviation. Conclusions: This study quantifies a normative range of median radiolunate interosseous proximities during wrist motion. The normative 4DCT datasets provide a reference range for evaluating postoperative arthrokinematics, suggesting the clinical-translational potential of dynamic imaging. The case study highlights applications for patient-level postoperative evaluation, which can be expanded to other pathologies, surgeries, or conditions in the future.
Morphology of the scaphotrapeziotrapezoid joint: A multi‐domain statistical shape modeling approach
Abstract The scaphotrapeziotrapezoid (STT) joint is involved in load transmission between the wrist and thumb. A quantitative description of baseline STT joint morphometrics is needed to capture the variation of normal anatomy as well as to guide staging of osteoarthritis. Statistical shape modeling (SSM) techniques quantify variations in three‐dimensional shapes and relative positions. The objectives of this study are to describe the morphology of the STT joint using a multi‐domain SSM. We asked: (1) What are the dominant modes of variation that impact bone and articulation morphology at the STT joint, and (2) what are the morphometrics of SSM‐generated STT joints? Thirty adult participants were recruited to a computed tomography study of normal wrist imaging and biomechanics. Segmentations of the carpus were converted to three‐dimensional triangular surface meshes. A multi‐domain, particle‐based entropy system SSM was used to quantify variation in carpal bone shape and position as well as articulation morphology. Articular surface areas and interosseous proximity distributions were calculated between mesh vertex pairs on adjacent bones within distance (2.0 mm) and surface‐normal angular (35°) thresholds. In the SSM, the first five modes of variation captured 76.2% of shape variation and contributed to factors such as bone scale, articular geometries, and carpal tilt. Median interosseous proximities—a proxy for joint space—were 1.39 mm (scaphotrapezium), 1.42 mm (scaphotrapezoid), and 0.61 mm (trapeziotrapezoid). This study quantifies morphological and articular variations at the STT joint, presenting a range of normative anatomy. The range of estimated interosseous proximities may guide interpretation of imaging‐derived STT joint space.
Fourier Analysis of the Vertical Ground Reaction Force During Walking: Applications for Quantifying Differences in Gait Strategies
Time series biomechanical data inform our understanding of normal gait mechanics and pathomechanics. This study examines the utility of different quantitative methods to distinguish vertical ground reaction forces (VGRFs) from experimentally distinct gait strategies. The goals of this study are to compare measures of VGRF data—using the shape factor method and a Fourier series-based analysis—to (1) describe how these methods reflect and distinguish gait patterns and (2) determine which Fourier series coefficients discriminate normal walking, with a relatively stiff-legged gait, from compliant walking, using deep knee flexion and limited vertical oscillation. This study includes a reanalysis of previously presented VGRF data. We applied the shape factor method and fit 3- to 8-term Fourier series to zero-padded VGRF data. We compared VGRF renderings using Euclidean L 2 distances and correlations stratified by gait strategy. Euclidean L 2 distances improved with additional harmonics, with limited improvement after the seventh term. Euclidean L 2 distances were greater in shape factor versus Fourier series renderings. In the 8 harmonic model, amplitudes of 9 Fourier coefficients—which contribute to VGRF features including peak and local minimum amplitudes and limb loading rates—were different between normal and compliant walking. The results suggest that Fourier series-based methods distinguish between gait strategies.