On gradient estimates for heat kernels