The answer is in general no because the zero set's convergence exponent imposes a lower bound on the order of growth of the function, which in turn imposes a growth condition on a subsequence of the Taylor coefficients.

In detail:

Assume $Z$ infinite of course and order all its non-zero elements (obviously $Z$ has finitely many zeroes and for each $r > 0$ there are only finitely many elements in $Z$ bounded by $r$ in absolute value) by the increasing absolute value of its elements $0 < |z_1| \leq |z_2| ... \leq |z_n| ...$, then the convergence exponent $\rho$ of $Z$ is the greatest lower bound of all $\lambda$ s.t. $\Sigma{\frac{1}{|z_n|^{\lambda}}}$ converges, where $\rho$ can be zero or infinity and we can find $Z$ to give us any such $\rho$ (for finite non-zero $\rho$ just pick $z_n = {n^{\frac{1}{\rho}}}$, for zero, $z_n = {2^n}$, for infinity, $z_n = n^{\frac{1}{loglogn}}$)

Assume now $f$ with zero set $Z$ has order $\rho$ as an entire function (finite or infinite, the order being the superior limit of $\frac{\log{\log{M_f(r)}}}{\log{r}}$, when $r$ goes to infinity, where $M_f(r)$ is the maximum of $|f|$ on the disc of radius $r$ centered at the origin), then a basic theorem (Hadamard) about entire functions says that the convergence exponent of its zero set defined above is less or equal than its order (obviously by multiplying $f$ with something like $e^{e^z}$ or $e^{z^q}$, $q$ high integer, we can always increase the order, even make it infinite, and keep same zero set, though in this case the order becomes integral or infinite as for finite non-integral orders there is a converse of the above).

But now there is a formula that gives the order of $f$ in terms of the Taylor coefficients and it essentially says that if $f(z) = \Sigma{b_nz^n}$ has order $\rho$, there is a subsequence of Taylor coefficients $c_n$ s.t $\frac{1}{|c_n}|$ is about $n^{\frac{n}{\rho})}$ - here we need $\rho$ non-zero, while for $\rho$ infinity we just interpret it as $\frac{1}{|c_n}|$ is smaller than $n^{\frac{n}{\rho})}$ for any finite $\rho$. Hence the sequence of Taylor coefficients of an entire function of non-zero order has a subsequence with a definite rate of growth ($\sqrt[n]{|c_n|}$ being about $n^{-\frac{1}{\rho}}$, where again we interpret that as usual for $\rho$ infinite), so in particular we have proved the negative answer to the original question as long as the zero set $Z$ has non-zero convergence exponent.