Self-similar properties of proton structure in the kinematic region of low Bjorken x are introduced and studied numerically. A description of the proton structure function F 2 (x, Q 2 ) reflecting self-similarity is proposed with a few parameters which are fitted to recent HERA data. The specific parameterisation provides an excellent description of the data which cover a region of four momentum transfer squared, 0.045 ≤ Q 2 ≤ 120 GeV 2 , and of Bjorken x, 6.2 · 10 −7 ≤ x ≤ 0.01 .Recent measurements of the H1 [1] and ZEUS [2] collaborations at HERA enable to study the proton structure in the region of low Q 2 1 GeV 2 where perturbative QCD has to face computation difficulties arising from the increase of the strong coupling constant α s (Q 2 ). Nevertheless, there is a number of approaches to describe the transition to low Q 2 at small x together with the region, Q 2 > 1 GeV 2 , which is described by perturbative QCD very well. Such attempts involve Reggeon exchange ideas [3], dipole interactions [4], vector meson dominance (VMD) [5], efficient parametrisations [6] and others.This letter presents a different point of view, based on the idea that the proton structure at low x is of fractal nature. Using the fractal dimension concept [7], a simple parametrisation of the proton structure function F 2 (x, Q 2 ) is obtained with a few well defined parameters. A numeric study is made using recent small x HERA data, for Q 2 between 0.045 GeV 2 and 120 GeV 2 .