Deep-inelastic heavy-quark production in electron-proton interactions at the HERA collider energies is considered within the QCD semihard (kT-factorization) approach. The dependence of the total and differential cross sections for the production of b quarks (and for muons arising in the subsequent semileptonic decay of B mesons) on the choice of unintegrated distributions of gluons in the proton is investigated. In particular, the gluon distributions obtained by solving the Catani-Ciafaloni-Fiorani-Marchesini evolution equation and the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi-Balitsky-Fadin-Kuraev-Lipatov combined equation are used in numerical calculations. The results of the theoretical calculations are compared with latest experimental data obtained by the H1 and ZEUS Collaborations at the HERA collider.