Great to see this question already answered! I was about to post a question / comment on exactly the same topic. Silly me, thought my question was very specific, and here it is, beautifully answered already, could've saved myself time thinking about it! I am coming from the physics side and was only vaguely aware of all the maths that deals with remainders of power series, having forgotten most of it.
My solution is as follows. We are talking about the Taylor series for the exponential,
$$
\exp(z)=\sum_{n=0}^{\infty}\frac{z^{n}}{n!},
$$
and I am interested in estimating the remainder (using the original notation) $e^z-P_n(z)$
$$
e^z-P_n(z) = \sum_{k=n+1}^{\infty}\frac{z^{k}}{k!}.
$$
It is not difficult to verify that the remainder is given by
$$
e^{z}-P_{n-1}(z)=\frac{z^{n}}{\left(n-1\right)!}\int_{0}^{\infty}\frac{\exp\left(\frac{zt}{1+t}\right)}{(1+t)^{n+1}}dt
$$
for $n\geq1$. For example, one can verify by integrating by parts that
$$
\int_{0}^{\infty}\frac{\exp\left(\frac{zt}{1+t}\right)}{(1+t)^{n+1}}dt=\frac{1}{n}+\frac{z}{n}\int_{0}^{\infty}\frac{\exp\left(\frac{zt}{1+t}\right)}{(1+t)^{n+2}}dt,
$$
where $n\geq1$. Then, by successively integrating by parts starting with
$$
\frac{z^{n}}{\left(n-1\right)!}\int_{0}^{\infty}\frac{\exp\left(\frac{zt}{1+t}\right)}{(1+t)^{n+1}}dt,
$$
one recovers the complete power series expansion for $e^{z}-P_{n-1}(z)$.
Now, for $\mathrm{Re}(z)<0$ can write
$$
\left|\frac{z^{n}}{\left(n-1\right)!}\int_{0}^{\infty}\frac{\exp\left(\frac{zt}{1+t}\right)}{(1+t)^{n+1}}dt\right|\leq\left|\frac{z^{n}}{\left(n-1\right)!}\right|\int_{0}^{\infty}\frac{\left|\exp\left(\frac{zt}{1+t}\right)\right|}{(1+t)^{n+1}}dt\leq\left|\frac{z^{n}}{\left(n-1\right)!}\right|\int_{0}^{\infty}\frac{dt}{(1+t)^{n+1}}=\frac{\left|z^{n}\right|}{n!}.
$$
Of course, the expression for the remainder holds for any $z$, but for $\mathrm{Re}(z)>0$ the exponential $\left|\exp\left(\frac{zt}{1+t}\right)\right|$ cannot be bounded by $1$ for all $0\leq t<\infty$.
By the way, I am sorry I had to present it like this - "It is not difficult...", "It is obvious..." - but I really just engineered this solution by hand, taking the Stieltjes asymptotic series as an example - see the Chapter on asymptotic series in "Advanced Mathematical Methods for Scientists and Engineers" by C.M. Bender and S.A. Orszag, McGraw-Hill (1978) for details.
Cheers!
P.S. For context, this result is quoted in books on quantum field theory when discussing how perturbative expansions (in field theory, but also in quantum mechanics in general) generate asymptotic series that most often do not converge, rather than convergent Taylor series with finite radii of convergence. For example, see "Condensed Matter Field Theory" by Altland and Simons or "Quantum Many-Particle Systems" by Negele and Orland. The following toy example is quoted as a typical ersatz-field-theoretic average:
$$
I(g)=\int_{-\infty}^{\infty}\frac{dx}{\sqrt{2\pi}}\,\exp\left(-\frac{1}{2}x^{2}-gx^{4}\right).
$$
Now, the perturbation theory approach would be to expand the integrand in powers of $g$, but such an expansion obviously has zero radius of convergence w.r.t. $g=0$ since the integral diverges for $g<0$. It is then stated then that we are dealing with an asymptotic series, so we should also be mindful of how many terms we wish to include in the expansion, since the error, given by the integral of the remainder of the Taylor series expansion of $\exp(-gx^{4})$), first decreases, and then increases by including more and more powers in the expansion. This is demonstrated by using the result quoted in this discussion, however none of the physics books justify it! - They simply pull it out of their magicians' hats!