not included in the original Balitsky-JIMWLK hierarchy and obviously missed in its mean field version, the BK equation. The fluctuations have been exhaustively investigated in recent years and their effects in the high energy evolution will be discussed in the next sections of this contribution. What one can antecipate is that they dramatically influence the high energy scattering, at least in the case when the strong coupling constant as is fixed. Among these investigations, one can cite a recent study [6] of the possible implications of the gluon number fluctuations in the measurements of the F2 proton structure function of deep inelastic scattering (DIS) at HERA, through well known models for the dipole-proton cross section in coordinate space [7,8]. The results show that, although they may be present in the DIS measurements, they seem not to strongly influence the evolution at HERA energies. In these proceedings we use the first model for the dipoleproton cross section developed in momentum space [9] with fluctuation effects included to fit the last HERA data. This model has shown to be successful in the description of DIS data, specially with heavy quark (charm) contributions included, by using a mean-field (BK) inspired parametrization for the scattering amplitude, which enters in the expression for the dipole-proton cross section. Some preliminary results are presented and some future expectations on the effects of fluctuations at HERA are discussed.