We derive some bounds for densities of random variables, starting from the Malliavin and Thalmaier representation formula in terms of the Riesz transform. We consider two examples of applications, both under a weak Hormander non-degeneracy hypothesis.