Relevance of the parallel nonlinearity in gyrokinetic simulations of tokamak plasmas