Here is another way to approach this problem.
The function $$f(z) = 1 - z/2 + z^2/3 + \ldots + (-1)^{k+1} z^k/(k+1) + \ldots $$ is analytic on the unit disc $ \{ z : |z| < 1\}$, which implies $ g(z) = \exp f(z)$ is also analytic on $ \{ z : |z| < 1\}$ and hence can be expanded as a power series $$g(z) = a_0 + a_1 z + a_2 z^2 + \dots + $$ in $ \{ z : |z| < 1\}$. We can easily compute the the first few coefficients as $a_0 = \exp f(0) = e$, $a_1 = \exp f(0) f^{'}(0) = -e/2$, and similarly $a_3 = 11e/24.$
However, $f(x) = \log(1+x)/x$ for for all real $x$ with $ 0 < x < 1$, so $g(x) = (1+x)^{1/x}$ for $ 0 < x < 1$ and the above series for real $x$ is an analytic extension of $(1+x)^{1/x}$ to $-1 < x < 1$.
Writing $$(1+x)^{1/x} = e - ex/2 + 11x^2/24 + \dots $$ from which we get
$(1+x)^{1/x} + cx = e + (c - e/2) x + 11ex^2/24 + \dots + \dots $.
The derivative of the above function at 0 is c -e/2, which is < 0, if c < e/2, by the continuity of the derivative, there is an interval $[0,\epsilon]$ on which the derivative of the function above is strictly negative and hence it decreases. Since 1/n decreases and lies in $[0,\epsilon]$ for all large $n$ this means $(1+1/n)^{n} + c/n$ increases eventually for any $ c < e/2$. This holds for any $x_n$ that strictly decreases to 0 not only $1/n,$, $(1+x_n)^{1/x_n} + c x_n$ eventually increases if $ c < e/2$. We can similarly argue that $(1+x_n)^{1/x_n} + c x_n$ eventually increases if $ c > e/2$ if $x_n$ strictly decreases to 0. For $c = e/2$, the positivity of the coefficient of $x^2$ implies $(1+x_n)^{1/x_n} + e x_n / 2$ eventually starts decreasing.