Specifically, I want to do this without invoking the Taylor/power series definitions for $e^z$. I know they exist and accept they're valid methods, but I want to show them from the definition of Laurent series, mostly just to see how it might work out - if at all (which goes into my question).
I know from the definition of a Laurent series given by
$$f(z) = \sum_{n = -\infty}^{\infty} a_k (z - z_0)^k$$
we have
$$a_n = \frac{1}{2 \pi i} \int_C \frac{f(\zeta)}{(\zeta - z_0)^{n+1}}d\zeta$$
where $C$ is a closed curve in the annulus in which the series is convergent in oriented counterclockwise.
For $f(z) = e^{1/z}$ have the discontinuity at $z = 0$ but other than that we don't have any "trouble points," so to speak. So we can define our annulus by $0 < |z| < \infty$, in effect.
We begin by trying to calculate the above integral, with a curve $C$ which is a circle of finite, positive radius. It becomes clear that
$$ \int_C \frac{e^{1/\zeta}}{\zeta^{n+1}}d \zeta$$
could be evaluated using the equality
$$\frac{1}{2\pi i} \int_C \frac{e^{1/\zeta}}{\zeta^{n+1}}d \zeta=\frac{1}{n!}f^{(n)}(0)$$
... but it won't work out, since the derivatives all have powers of $z$ in the denominator (thus running into problems with inputs of $0$). I certainly can't think of any other way (with our current skillset in the class) as to how it might work.
Is this just a case where it's like "you just have to use the power series for $e^z$, substitute in $1/z$, and invoke the uniqueness of Laurent series," or is there something I'm overlooking?