Closed-form of $\sum_{n=1}^\infty\frac{(-1)^{n+1}}{n}\Psi_3(n+1)=-\int_0^1\frac{\ln(1+x)\ln^3 x}{1-x}\,dx$ Does the following series or integral have a closed-form

\begin{equation}
\sum_{n=1}^\infty\frac{(-1)^{n+1}}{n}\Psi_3(n+1)=-\int_0^1\frac{\ln(1+x)\ln^3 x}{1-x}\,dx
\end{equation}

where $\Psi_3(x)$ is the polygamma function of order $3$.

Here is my attempt. Using equation (11) from Mathworld Wolfram:
\begin{equation}
\Psi_n(z)=(-1)^{n+1} n!\left(\zeta(n+1)-H_{z-1}^{(n+1)}\right)
\end{equation}
I got
\begin{equation}
\Psi_3(n+1)=6\left(\zeta(4)-H_{n}^{(4)}\right)
\end{equation}
then
\begin{align}
\sum_{n=1}^\infty\frac{(-1)^{n+1}}{n}\Psi_3(n+1)&=6\sum_{n=1}^\infty\frac{(-1)^{n+1}}{n}\left(\zeta(4)-H_{n}^{(4)}\right)\\
&=6\zeta(4)\sum_{n=1}^\infty\frac{(-1)^{n+1}}{n}-6\sum_{n=1}^\infty\frac{(-1)^{n+1}H_{n}^{(4)}}{n}\\
&=\frac{\pi^4}{15}\ln2-6\sum_{n=1}^\infty\frac{(-1)^{n+1}H_{n}^{(4)}}{n}\\
\end{align}
From the answers of this OP, the integral representation of the latter Euler sum is
\begin{align}
\sum_{n=1}^\infty\frac{(-1)^{n+1}H_{n}^{(4)}}{n}&=\int_0^1\int_0^1\int_0^1\int_0^1\int_0^1\frac{dx_1\,dx_2\,dx_3\,dx_4\,dx_5}{(1-x_1)(1+x_1x_2x_3x_4x_5)}
\end{align}
or another simpler form
\begin{align}
\sum_{n=1}^\infty\frac{(-1)^{n+1}H_{n}^{(4)}}{n}&=-\int_0^1\frac{\text{Li}_4(-x)}{x(1+x)}dx\\
&=-\int_0^1\frac{\text{Li}_4(-x)}{x}dx+\int_0^1\frac{\text{Li}_4(-x)}{1+x}dx\\
&=\text{Li}_5(-1)-\int_0^{-1}\frac{\text{Li}_4(x)}{1-x}dx\\
\end{align}
I don't know how to continue it, I am stuck. Could anyone here please help me to find the closed-form of the series preferably with elementary ways? Any help would be greatly appreciated. Thank you.

Edit :
Using the integral representation of polygamma function
\begin{equation}
\Psi_m(z)=(-1)^m\int_0^1\frac{x^{z-1}}{1-x}\ln^m x\,dx
\end{equation}
then we have
\begin{align}
\sum_{n=1}^\infty\frac{(-1)^{n+1}}{n}\Psi_3(n+1)&=-\sum_{n=1}^\infty\frac{(-1)^{n+1}}{n}\int_0^1\frac{x^{n}}{1-x}\ln^3 x\,dx\\
&=-\int_0^1\sum_{n=1}^\infty\frac{(-1)^{n+1}x^{n}}{n}\cdot\frac{\ln^3 x}{1-x}\,dx\\
&=-\int_0^1\frac{\ln(1+x)\ln^3 x}{1-x}\,dx\\
\end{align}
I am looking for an approach to evaluate the above integral without using residue method or double summation.

Computation of $\displaystyle U=\int_0^1 \frac{\ln(1+x)\ln^3 x}{1-x}\,dx$
\begin{align*}
U&\overset{\text{IBP}}=\left[\left(\int_0^x \frac{\ln^3 t}{1-t}\,dt\right)\ln(1+x)\right]_0^1-\int_0^1 \frac{1}{1+x}\left(\int_0^x\frac{\ln^3 t}{1-t}\,dt\right)\,dx\\
&=-6\zeta(4)\ln 2+\int_0^1\int_0^1 \left(\frac{\ln^3(tx)}{(1+t)(1+x)}-\frac{\ln^3(tx)}{(1+t)(1-tx)}\right)\,dt\,dx\\
&=-6\zeta(4)\ln 2+6\left(\int_0^1\frac{\ln^2 t}{1+t}\,dt\right)\left(\int_0^1\frac{\ln x}{1+x}\,dx\right)+\\
&2\left(\int_0^1\frac{\ln^3 t}{1+t}\,dt\right)\left(\int_0^1\frac{1}{1+x}\,dx\right)-\int_0^1 \frac{1}{t(1+t)}\left(\int_0^t \frac{\ln^3 u}{1-u}\,du\right)\,dt\\
&=-\frac{33}{2}\zeta(4)\ln 2-\frac{9}{2}\zeta(2)\zeta(3)-\int_0^1 \frac{1}{t(1+t)}\left(\int_0^t \frac{\ln^3 u}{1-u}\,du\right)\,dt\\
&\overset{\text{IBP}}=-\frac{33}{2}\zeta(4)\ln 2-\frac{9}{2}\zeta(2)\zeta(3)-\left[\ln\left(\frac{t}{1+t}\right)\left(\int_0^t \frac{\ln^3 u}{1-u}\,du\right)\right]_0^1+\int_0^1 \frac{\ln\left(\frac{t}{1+t}\right)\ln^3 t}{1-t}\,dt\\
&=-\frac{45}{2}\zeta(4)\ln 2-\frac{9}{2}\zeta(2)\zeta(3)+\int_0^1 \frac{\ln\left(\frac{t}{1+t}\right)\ln^3 t}{1-t}\,dt\\
&=-\frac{45}{2}\zeta(4)\ln 2-\frac{9}{2}\zeta(2)\zeta(3)+24\zeta(5)-U\\
U&=\boxed{-\frac{45}{4}\zeta(4)\ln 2-\frac{9}{4}\zeta(2)\zeta(3)+12\zeta(5)}
\end{align*}
