Intracellular Vesicle Transport Impairment as a Candidate Systems-Level Bottleneck in Chronic Diabetic Foot Ulcers: Network Medicine Identifies KIF13A as a Potential Therapeutic Vulnerability.

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 (log2FC = 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.