We have extended our supercomputer-enabled Monte Carlo simulations of hopping transport in completely disordered 2D conductors to the case of substantial electron-electron Coulomb interaction. Such interaction may not only suppress the average value of hopping current, but also affect its fluctuations rather substantially. In particular, the spectral density S(I)(f) of current fluctuations exhibits, at sufficiently low frequencies, a 1/f-like increase which approximately follows the Hooge scaling, even at vanishing temperature. At higher f, there is a crossover to a broad range of frequencies in which S(I)(f) is nearly constant, hence allowing characterization of the current noise by the effective Fano factor [Formula: see text]. For sufficiently large conductor samples and low temperatures, the Fano factor is suppressed below the Schottky value (F = 1), scaling with the length L of the conductor as F = (L(c)/L)(α). The exponent α is significantly affected by the Coulomb interaction effects, changing from α = 0.76 ± 0.08 when such effects are negligible to virtually unity when they are substantial. The scaling parameter L(c), interpreted as the average percolation cluster length along the electric field direction, scales as [Formula: see text] when Coulomb interaction effects are negligible and [Formula: see text] when such effects are substantial, in good agreement with estimates based on the theory of directed percolation.