Although DFTB3 is a widely used semi-empirical method to study (bio)organic systems and often serves as a viable alternative to DFT methods for QM/MM molecular dynamics simulations, it has significant limitations in accurately describing the energetics of hydride transfers. This limitation manifests as an artificial local minimum along the reaction coordinate instead of a transition state barrier, as demonstrated for Ccr and DHFR. To address this systematic deficiency, we introduced a straightforward, targeted approach based on analyzing C-H and C-C distance distributions, followed by modifications to the repulsive potentials within specific regions using linear combinations of harmonic functions.