Investigation into the limits of perturbation theory at low $Q^2$ using HERA deep inelastic scattering data