Intracellular Vesicle Transport Impairment as a Candidate Systems-Level Bottleneck in Chronic Diabetic Foot Ulcers:...

Background: Growth factor therapy often fails in diabetic foot ulcers (DFUs). The reason remains unclear. Standard differential expression analysis may miss functionally critical genes with modest expression changes. Methods: We performed a secondary computational analysis of a longitudinal DFU transcriptomic dataset (Dryad; 17 patients, 117 serial biopsy samples, 12-week follow-up). Co-expression networks were built separately for healed (n = 37) and non-healed (n = 80) samples. Virtual gene knockout (VGK) was used to rank genes by topological impact on network cohesion. Single-cell analysis (GSE165816) assessed the association between endogenous KIF13A expression and keratinocyte migration-related signatures. A conceptual Hill-equation simulation was used to illustrate the transport-signaling threshold relationship. Drug repurposing used DSigDB enrichment. An independent bulk DFU cohort (GSE134431) was used for external validation. Results: KIF13A showed no differential expression (log₂FC = 0.173, p = 0.263) yet ranked first by VGK topological impact. In keratinocytes, high KIF13A expression correlated with greater migration scores versus zero-detection cells (p = 0.0058). A clear threshold effect emerged: below the 30th expression percentile, EGF, PDGF, and FGF pathway activation scores remained near baseline. In a structural-equation model, transport activity negatively predicted inflammation (standardized β = -0.92, p < 0.001). HIF1A showed the strongest positive correlation with KIF13A in keratinocytes (Spearman ρ = 0.26, p < 0.001), and FOS showed a negative correlation in the single-cell analysis (ρ = -0.16, p < 0.001) and in the bulk longitudinal cohort (ρ = -0.32, p < 0.001, n = 117). Recurrent AKR1B1-related drug signatures nominated the aldose-reductase pathway, and epalrestat was therefore prioritized as a hypothesis-generating candidate compound rather than a direct top-ranked enrichment hit. External validation confirmed consistent upregulation of KIF13A (Fold-Change = 1.58, adj. p = 0.0075), EPN1, and CLIP1 in DFU tissue. Despite population-level upregulation, a subset of cells fell below the functional signaling threshold. Conclusions: These computational findings suggest that KIF13A-associated vesicle transport impairment may represent a candidate systems-level bottleneck for growth-factor responsiveness in DFUs, a network-level pattern not captured by standard differential-expression analysis. Epalrestat, an AKR1B1 inhibitor prioritized through recurrent AKR1B1-related drug signatures, is presented as a candidate compound for further evaluation. As the present analysis is observational and computational, the findings should be interpreted as hypothesis-generating; experimental perturbation studies and prospective clinical validation are required.