It is shown that in the leading twist approximation of the Wilson operator product expansion with "frozen" and analytic strong coupling constants, Bessel-inspired behavior of the structure functions F 2 and F cc 2 and also the derivative ∂ ln F 2 /∂ ln(1/x) at small x values, obtained for a flat initial condition in the DGLAP evolution equations, leads to good agreement with the deep inelastic scattering experimental data from HERA.