At present, optical anapole resonances in nanostructures have attracted increasing attention due to the strong field confinement and substantially suppressed scattering. This study provides the results of three-dimensional finite-difference time-domain simulations exhibiting the possibility of the anapole state in high refractive index dielectric nanohole array metasurfaces having different profiles of the holes (square, hexagon, and octagon). Behavior of the effective optical Kerr nonlinearity of the metasurfaces in the vicinity of the anapole state is investigated. Depending on the geometry, the absolute value of the effective nonlinear Kerr coefficient of the metasurface may be up to three orders of magnitude greater than that of the unstructured film. A square transverse section of the nanohole is preferable for the optical Kerr effect in the holey metasurfaces. The effect of the random rotation of the square holes representing the metasurface irregularity on the optical nonlinearity is examined. As a result, the dielectric nanohole array metasurfaces display a concrete possibility to have the anapole state with large enhancement of the optical nonlinearity.