In complete analogy with the bilateral Laplace transform of x(t)=e−at (which doesn't exist), the bilateral Z-transform of an doesn't exist either. The series
∑n=−∞∞anz−n
converges nowhere, simply because an grows without bounds for n→−∞ if |a|<1, or for n→∞ if |a|>1. Of course, for |a|=1 there series doesn't converge either.
EDIT:
As for your computation of the Z-transform of an, the mistake lies in the fact that in addition to the algebraic expression of the transform you also need to consider the region of convergence. If you split an (as you did) as
an=anu[n]+anu[−n−1](1)
you can compute the Z-transform of both right-hand side expressions separately:
Z{anu[n]}=zz−a,|z|>|a|Z{anu[−n−1]}=−zz−a,|z|<|a|(2)
Note that the region of convergence (ROC) for the first part is outside the circle with radius |a|, whereas the ROC of the second part is inside the circle with radius |a|. The ROC of the total expression would be the overlap of the two ROCs, which is zero. Consequently, the sum doesn't converge anywhere and the Z-transform of the total expression doesn't exist.