Potential analysis for a class of diffusion equations: A Gaussian bounds approach