Taking the $Z-$transform for the equation that you proved, i.e.
\begin{equation}
u_{n+1}=-\frac{3}{4}+\frac{3}{4}(n+1)u_n
\end{equation}
\begin{equation}
zU(z)-zu(0)=-\frac{3}{4}\frac{z}{z-1}-\frac{3z}{4}\frac{\text{dU(z)}}{\text{dt}}+\frac{3}{4}U(z)
\end{equation}
Assuming $u(0)=1$, and after some manipulation of the above equation and taking inverse $Z-$transform, we get,
\begin{equation}
u_n=\frac{1}{4}\Bigg((4c_1+3)\Big(\frac{3}{4}\Big)^n\Gamma(n+1)-4e^\frac{4}{3}E_{-n}\Big(\frac{4}{3}\Big)\Bigg)
\end{equation}
To find $c_1$, we substitute $u_0=1$ and use the asymptotic expansion of $E_n(x)$ viz.,
\begin{equation}
E_n(x)= \frac{e^{-x}}{x}\Bigg[1-\frac{n}{x}+\frac{n(n+1)}{x^2}-...\Bigg]
\end{equation}
\begin{equation}
E_0\Big(\frac{4}{3}\Big)=\frac{3}{4}e^{-\frac{4}{3}}
\end{equation}
Thus, $c_1=1$
Thus,
\begin{equation}
u_n=\frac{1}{4}\Bigg(7\cdot\Big(\frac{3}{4}\Big)^n\Gamma(n+1)-4e^{\frac{4}{3}}E_{-n}\Big(\frac{4}{3}\Big)\Bigg).
\end{equation}
Now, replacing $n$ with $-n$ and $x$ by $\frac{4}{3}$ in the Exponential Integral function $E_n(x)$, we get,
\begin{equation}
u_n=\frac{1}{4}\Bigg(7\cdot\Big(\frac{3}{4}\Big)^nn!-3\displaystyle \sum_{k=0}^{n} \Bigg(\frac{3}{4}\Bigg)^k{}^nP_k\Bigg)
\end{equation}
This can be represented in other form as:
\begin{equation}
y_n=\frac{4^{-n}\Big(7\cdot 3^{n+1}\cdot (n+1)!-4^{n+1}e^{4/3}E_{-n}\Big(\frac{4}{3}\Big)(n+3)\Big)}{3(n+1)}
\end{equation}
where, $ y_n = 4u_n-\frac{4}{n+1}$.
Since this is a strictly increasing function, it is proved that $u_n\ge\frac{1}{n+1}$

Wolframalpha plot for $y_n$
\frac
in exponents or limits of integrals. It looks bad and confusing, and it rarely appears in professional mathematics typesetting. $\endgroup$