A rigorous approach for constructing of energy density functionals is developed using the local-scale point transformation method. The exact kinetic energy density τ[ρ] contains two terms: the original Weizsäcker term τw[ρ] and a Thomas-Fermi-like term τw[ρ]. An energy density functionalE[ρ] is built up within a Slater determinant approximation using the most general Skyrme-type forces. It is demonstrated that all Hartree-Fock nuclear ground state results are reproduced quite accurately and consequentlyE[ρ] includes the shell effects accounted mostly by the particular form of τ[ρ]. The approach presented is compared with both Extended Thomas-Fermi and Expectation Value methods. Further possible applications are briefly considered.