Considering the algebraic identity
\begin{align*}
&(a-b)^3b = a^3b - 3a^2b^2 + 3ab^3 - b^4 = -2a^3b +3(a^3b+ab^3) -3a^2b^2 -b^4\\
&\Longrightarrow \ \ \ 2a^3b = -{b^4 \over 2} -{b^4 + 6a^2b^2\over 2} + 3(a^3b+ab^3) - (a-b)^3b
\end{align*} with $a = \ln(1-x)$ and $b= \ln (1+x)$ it follows that
\begin{align*}
2\int_0^1 {\ln^3(1-x)\ln(1+x)\over x}dx =& - \frac 1 2\int_0^1 {\ln^4(1+x)\over x}d x \\
&-\frac 12 \int_0^1 \frac{\ln^4(1+x) + 6\ln^2(1-x)\ln^2(1+x)}{x}dx\\
&+3\int_0^1 \frac{\ln^3(1-x)\ln(1+x) + \ln(1-x)\ln^3(1+x)}{x}dx\\
&- \int_0^1 \frac{\ln^3\left(\frac{1-x}{1+x}\right)\ln(1+x)}{x}dx\\
=:& -I_1 - I_2 + I_3 -I_4.
\end{align*}
For $I_1$, make substitution $y = \frac x {1+x}$ to get:
\begin{align*}
I_1 =& \frac 1 2 \int_0^{\frac 12} \frac{\ln^4(1-y)}{y(1-y)} dy \\
=& \frac 1 2\underbrace{ \int_0^{\frac 12} \frac{\ln^4(1-y)}{y} dy}_{z=1-y}+ \frac 1 2 \int_0^{\frac 12} \frac{\ln^4(1-y)}{1-y} dy\\
=& \frac 1 2 \int_{\frac 1 2 }^1 \frac{\ln^4 z} {1-z} dz + \frac {\ln^5 2}{10}\\
=& \frac 12 \sum_{n=1}^\infty \int_{\frac 1 2}^1 z^{n-1}\ln^4 z\ dz + \frac {\ln^5 2}{10}\\
=& \frac 12 \sum_{n=1}^\infty \frac{\partial^4}{\partial n^4}\left[\frac 1 n - \frac 1 {n2^n}\right] + \frac {\ln^5 2}{10}\\
=& \frac 12 \sum_{n=1}^\infty \left[\frac{24}{n^5} - \frac {24}{n^52^n} - \frac{24 \ln 2}{n^42^n}-\frac{12\ln^2 2}{n^3 2^n}-\frac{4\ln^3 2}{n^2 2^n} - \frac{\ln^4 2}{n2^n}\right] + \frac {\ln^5 2}{10}\\
=&12\zeta(5) - 12\text{Li}_5(1/2) - 12\ln 2 \text{Li}_4(1/2) -6\ln^2 2 \text{Li}_3(1/2) -2\ln^3 2\text{Li}_2(1/2)-\frac {2}{5}\ln^5 2\\
=&\boxed{-12\Big(\text{Li}_5(1/2) + \ln 2\text{Li}_4(1/2)-\zeta(5)\Big)-{21 \over 4}\zeta(3)\ln^2 2 +{1\over 3} \pi^2 \ln^3 2-{2 \over 5} \ln^5 2}
\end{align*} where the well-known values
\begin{align*}\text{Li}_2(1/2) = {\pi^2 \over 12}-{\ln^2 2\over 2} , \qquad \text{Li}_3(1/2) ={7\zeta(3) \over 8} -{\pi^2 \ln 2\over 12} + {\ln^3 2 \over 6}
\end{align*} are used.
Actually, $I_2$ was already evaluated by the OP here using the algebraic identity $$b^4 + 6a^2b^2 = \frac {(a-b)^4} 2+\frac{(a+b)^4}{2} -a^4.$$
It holds that
$$
\boxed{I_2 = \frac {21}{8} \zeta(5).}
$$
In fact, the value of $I_3$ can also be found in the previous answer of @Przemo's. For $I_3$, one can use the algebraic relation $3(a^3b + ab^3) =\frac 3 8 \left[ (a+b)^4 - (a-b)^4\right]$.
This gives
\begin{align*}
I_3=& \underbrace{\frac 3 8 \int_0^1 \frac{\ln^4(1-x^2)}{x} dx}_{x^2 = y} - \underbrace{\frac 3 8 \int_0^1 \frac{\ln^4\left(\frac{1-x}{1+x}\right)}{x} dx}_{\frac{1-x}{1+x} = y}\\
=&\frac 3 {16}\underbrace{\int_0^1 \frac{\ln^4(1-y)}{y} dy }_{1-y\mapsto y}- \frac 3 4 \int_0^1 \frac{\ln^4 y}{1-y^2} dy\\
=&\frac 3 {16}\int_0^1 \frac{\ln^4 y}{1-y} dy - \frac 3 4 \sum_{n=0}^\infty \int_0^1 y^{2n} \ln^4 y \ dy\\
=&\frac 3 {16}\sum_{n=1}^\infty \int_0^1 y^{n-1}\ln^4 y \ dy - \frac 3 4 \sum_{n=0}^\infty \frac {24}{(2n+1)^5}\\
=&\frac 3 {16}\sum_{n=1}^\infty \frac{24}{n^5} - 18 \sum_{n=0}^\infty \frac {1}{(2n+1)^5}\\
=&\frac {9}{2} \zeta(5)- 18\cdot \frac {31}{32}\zeta(5)\\
=&\boxed{-\frac{207}{16}\zeta(5)}
\end{align*} as can be found in @Przemo's answer.
For $I_4$, make substitution $ \frac{1-x}{1+x}\mapsto x$ to get
\begin{align*} I_4 = &2\int_0^1 \frac{\ln^3 x \ln\left(\frac 2 {1+x}\right)}{1-x^2} dx \\
=&2\ln 2 \int_0^1 \frac{\ln^3 x}{1-x^2} dx - \underbrace{2\int_0^1\frac{\ln^3 x \ln(1+x)}{1-x^2} dx }_{=:J}\\
=& 2\ln 2\sum_{n=0}^\infty \int_0^1 x^{2n} \ln^3 x\ dx - J\\
=& - 12\ln 2 \underbrace{\sum_{n=0}^\infty \frac 1 {(2n+1)^4}}_{\frac{15}{16}\zeta(4) = \frac{\pi^4}{96}} - J \\
=& -\frac{\pi^4 \ln 2}{8} - J.
\end{align*}
\begin{align*}
J = &\int_0^1\frac{2\ln^3 x \ln(1+x)}{1-x^2} dx \\
=& \underbrace{\int_0^1 \frac{\ln^3 x \ln(1+x)}{1+x}dx}_{=:A} + \int_0^1 \frac{\ln^3 x \ln(1+x)}{1-x}dx\\
=& A + \int_0^1 \frac{\ln^3 x \ln(1-x^2)}{1-x}dx -\int_0^1 \frac{\ln^3 x \ln(1-x)}{1-x}dx\\
=&A + \int_0^1 \frac{(1+x)\ln^3 x \ln(1-x^2)}{1-x^2}dx -\int_0^1 \frac{\ln^3 x \ln(1-x)}{1-x}dx\\
=&A + \underbrace{\int_0^1 \frac{\ln^3 x \ln(1-x^2)}{1-x^2}dx }_{=:B}+\underbrace{\int_0^1 \frac{x\ln^3 x \ln(1-x^2)}{1-x^2}dx}_{x^2 \mapsto x}-\int_0^1 \frac{\ln^3 x \ln(1-x)}{1-x}dx\\
=&A + B - \underbrace{\frac {15}{16} \int_0^1 \frac{\ln^3 x \ln(1-x)}{1-x}dx}_{=:C}\\
=&A + B - C.
\end{align*}
For $A$, we can use the McLaurin series of
$$
\frac{\ln (1+x)}{1+x} = \sum_{n=0}^\infty (-1)^{n-1}H_n x^n
$$ ($H_0= 0$) to get
\begin{align*}
A = & \sum_{n=0}^\infty (-1)^{n-1}H_n \int_0^1 x^n\ln^3 x \ dx \\
=&6 \sum_{n=0}^\infty \frac{(-1)^{n}H_n}{(n+1)^4}\\
=&6 \sum_{n=0}^\infty \frac{(-1)^{n}H_{n+1}}{(n+1)^4} - 6\sum_{n=0}^\infty \frac{(-1)^{n}}{(n+1)^5}\\
=&6 \sum_{n=1}^\infty \frac{(-1)^{n-1}H_{n}}{n^4} - 6\sum_{n=1}^\infty \frac{(-1)^{n-1}}{n^5}\\
=& 6\left(\frac{59}{32}\zeta(5) - \frac{\pi^2\zeta(3)}{12}\right)-6\cdot \frac{15}{16}\zeta(5)\\
=& \frac{87}{16}\zeta(5) - \frac{\pi^2 \zeta(3)}{2}.
\end{align*}
Here, the known value of $ \sum_{n=1}^\infty (-1)^{n-1}{H_n \over n^4}$ is used.
For $B$, make substitution $u = x^2$ to get
\begin{align*}
B =& \frac 1 {16} \int_0^1 \frac{\ln^3 u \ln(1-u)}{\sqrt u (1-u)} du \\
=& \frac 1 {16} \left[\frac{\partial^4}{\partial x^3\partial y} \text{B}(x,y)\right]_{x=\frac 1 2, y = 0^+}
\end{align*} where $\text{B}(\cdot,\cdot)$ is Euler's Beta function. We can use the fact that
\begin{align*}
\lim_{y\to 0^+}\frac{\partial^2}{\partial x\partial y} \text{B}(x,y) = -\frac 1 2 \psi''(x) + \psi'(x) \big[\psi(x) + \gamma\big]
\end{align*} to get
\begin{align*}
B =& \frac 1 {16}\frac{d^2}{dx^2}\left[-\frac 1 2 \psi''(x) + \psi'(x) \big[\psi(x) + \gamma\big]\right]_{x=\frac 1 2}\\
=&\frac 1 {16} \left[-\frac 1 2 \psi''''(1/2) + \psi'''(1/2)\big[\psi(1/2) + \gamma\big] + 3\psi'(1/2)\psi''(1/2)\right]\\
=& \frac 1 {16}\left[-21\pi^2 \zeta(3) + 372\zeta(5) - 2\pi^4 \ln 2\right]
\end{align*} which can be evaluated using the series representations of polygamma functions $$\psi(x) +\gamma = - \frac 1 x +\sum_{n=1}^\infty \frac 1 n - \frac 1 { n+x},\\
\psi'(x) = \sum_{n=0}^\infty \frac 1 {(n+x)^2}$$ and the derived fact that $\psi(\tfrac 1 2 )+\gamma = -2\ln 2$ and $\psi^{(k)}(\tfrac 1 2)=(-1)^{k+1}k!(2^{k+1}-1)\zeta(k+1)$ for $k\ge 1$.
For $C$, we can use the same method as used in the evaluation of $B$. It holds that
\begin{align*}
C =& \frac {15}{16} \left[\frac{\partial^4}{\partial x^3\partial y} \text{B}(x,y)\right]_{x=1, y = 0^+}\\
=&\frac {15} {16}\left[-\frac 1 2 \psi''''(1) + \psi'''(1)\big[\psi(1) + \gamma\big] + 3\psi'(1)\psi''(1)\right]\\
=&\frac{15}{16}\left[12\zeta(5) -6\zeta(2)\zeta(3)\right]\\
=&\frac {45}{4}\zeta(5) -\frac {15\pi^2 \zeta(3)}{16}
\end{align*} where $\psi(1) +\gamma = 0$, $\psi'(1) = \zeta(2)$, $\psi''(1) = -2\zeta(3)$ and $\psi''''(1) = -24\zeta(5)$ are used.
Combining $A,B,C$, we have that $$J =A+B-C= \frac{279}{16}\zeta(5) -\frac{7\pi^2\zeta(3)}{8} - \frac{\pi^4 \ln 2}{8}$$ and
$$
\boxed{I_4 = -\frac{\pi^4 \ln 2}{8} - J = -\frac{279}{16}\zeta(5)+\frac{7\pi^2\zeta(3)}{8}}
$$
Finally, these evaluate $\int_0^1 {\ln^3(1-x)\ln(1+x)\over x}dx =\frac 1 2\big[-I_1-I_2+I_3-I_4\big]$ as follows.
\begin{align*}
\int_0^1 {\ln^3(1-x)\ln(1+x)\over x}dx =&\ 6\text{Li}_5(1/2) + 6\ln 2\ \text{Li}_4(1/2)-\frac{81}{16}\zeta(5)-{7\pi^2 \over 16}\zeta(3)\\
&+\frac{21\ln^2 2}{8}\zeta(3)- \frac{1}{6}\pi^2\ln^3 2+\frac{1}{5}\ln^5 2.
\end{align*}
Using the identity given in the OP, we get the desired integral $I$
\begin{align*}
\int_0^{\frac 1 2}\frac{\text{Li}_2^2(x)}{x} dx = &-2\text{Li}_5(1/2) -2\ln 2\ \text{Li}_4(1/2)+\frac{27}{32}\zeta(5) +\frac{7\pi^2}{48}\zeta(3)-\frac{7\ln^2 2}{8}\zeta(3) \\
&-\frac{\pi^4\ln 2}{144} +\frac{\pi^2\ln^3 2}{12} - \frac{7\ln^5 2}{60}.
\end{align*}
We will look into the integral
\begin{align*}
I = &\int_{0}^{1}\arcsin^4 x\frac{ \ln x}{\sqrt{1-x^2}}\ \mathrm dx
\end{align*} taking the @nospoon's novel approach presented here. Using the MacLaurin series of $\arcsin^4 x$
$$
\arcsin^4 x =\frac 3 2 \sum_{n=1}^\infty \frac{4^{n}H_{n-1}^{(2)}}{n^2{2n \choose n}}x^{2n}
$$
and the fact that
$$
\small\operatorname{B}(n+\tfrac 1 2,\tfrac 1 2) = \int_0^1 x^{n-1/2}(1-x)^{-1/2}\ \mathrm dx = 2\int_0^{\frac\pi 2} \sin^{2n}\theta\ \mathrm d\theta = \frac{\pi}{4^n}{2n \choose n},\tag{$\small x\mapsto \sin^2\theta$}
$$
\begin{align*}\small
\psi(n+\tfrac 12 ) -\psi(n+1) =&\small \sum_{k=1}^\infty \frac 1{\scriptsize k+n} - \frac 1{\scriptsize k+n-\tfrac 1 2} \\
=&\small\sum_{k=1}^\infty \left(\frac 1{\scriptsize k} - \frac 1{\scriptsize k-\tfrac 1 2}\right)-\sum_{k=1}^n\frac 1 {\scriptsize k} + \sum_{k=1}^n\frac 1{\scriptsize k-\tfrac 1 2}\\
=&\small-2\ln 2 -H_n +2(H_{2n}-\tfrac 1 2H_n)\\
=&\small 2(H_{2n}-H_n-\ln 2),
\end{align*}
\begin{align*}
\Longrightarrow \ {\int_{ 0}^{1 }x^{2n}\frac{ \ln x}{\sqrt{1-x^2}}\ \mathrm dx} = & \frac 1 4\int_{0 }^{1 }x^{n-1/2} { \ln x \over \sqrt{1-x}}\ \mathrm dx\tag{$\small x^2\mapsto x$}\\
=& \frac 1 4 \left[\frac{\partial }{\partial x}\operatorname{B}(x,y) \right]_{x=n+1/2,y=1/2}\\
=&\frac 1 4\Big[ \operatorname{B}(x,y)\big[\psi(x) -\psi(x+y) \big]\Big]_{x=n+1/2,y=1/2}\\
=& \frac 1 4 \operatorname{B}(n+\tfrac 1 2,\tfrac 1 2)\big[\psi(n+\tfrac 12 ) -\psi(n+1) \big]\\
=& \frac{\pi}2\frac{{2n \choose n}}{4^{n}} \left(H_{2n} - H_n -\ln 2\right),
\end{align*} where $\operatorname{B}(x,y)$ and $\psi(x)$ are the Beta and digamma function, respectively, we have
\begin{align*}
I = &\frac 3 2\sum_{n=1}^\infty \frac{4^{n}H_{n-1}^{(2)}}{n^2{2n \choose n}}\int_{0}^{1}x^{2n}\frac{ \ln x}{\sqrt{1-x^2}}\ \mathrm dx \\
=&\frac {3\pi}4 \sum_{n=1}^\infty \frac{H^{(2)}_{n-1}}{n^2}\left(H_{2n} - H_n -\ln 2\right) \\
=&\frac {3\pi}4\sum_{n=1}^\infty \frac{H^{(2)}_{n-1}H_{2n}}{n^2}-\frac {3\pi}4\underbrace{\sum_{n=1}^\infty \frac{H^{(2)}_{n-1}H_{n}}{n^2}}_{=-2\zeta(5) +2\zeta(2)\zeta(3)}-\frac {3\pi\ln 2}4\underbrace{\sum_{n=1}^\infty \frac{H^{(2)}_{n-1}}{n^2}}_{=\frac{3}4 \zeta(4)}\\
=&\frac{3\pi}{4} \sum_{n=1}^\infty \frac{H^{(2)}_{n}H_{2n}}{n^2} -\frac{3\pi}4\underbrace{\sum_{n=1}^\infty \frac{H_{2n}}{n^4}}_{=\frac{37}{4}\zeta(5)-4\zeta(2)\zeta(3)} +\frac{3\pi}2 \zeta(5) -\frac{\pi^3}4\zeta(3) -\frac{\pi^5\ln 2}{160}\\
=&\boxed{3\pi S -\frac{87\pi}{16} \zeta(5) +\frac{\pi^3}{4}\zeta(3) -\frac{\pi^5\ln 2}{160}}
\end{align*} where $S = \sum_{n=1}^\infty \frac{H_{2n}H^{(2)}_{n}}{4n^2}$ is the sum in question, and the known values of several Euler sums
$$
\sum_{n=1}^\infty \frac{H^{(2)}_{n-1}H_{n}}{n^2}=-2\zeta(5) +2\zeta(2)\zeta(3),\tag{1}
$$
$$\sum_{n=1}^\infty \frac{H^{(2)}_{n}}{n^2}=\frac{7}4 \zeta(4),\tag{2}
$$
\begin{align*}\sum_{n=1}^\infty \frac{H_{2n}}{n^4} =& 8\sum_{n=1}^\infty \frac{H_{n}}{n^4}-8\sum_{n=1}^\infty \frac{(-1)^{n-1} H_{n}}{n^4}\\
=&8\big(3\zeta(5)-\zeta(2)\zeta(3)\big)-8\left(\frac{59}{32}\zeta(5)-\frac 1 2\zeta(2)\zeta(3)\right)\\
=&\frac{37}4\zeta(5) - 4\zeta(2)\zeta(3)\tag{3}
\end{align*} are used.
Note: $(1)$ is in @nospoon's answer here, $(2)$ can be found here, and for $(3)$ you can see Euler's formula and here.
Evaluation of $I$: By making substitution $x = \sin \theta$ and using the Fourier series of
$$
\ln (\sin\theta) = -\ln 2 -\sum_{k=1}^\infty \frac{ \cos(2k \theta)}{k},
$$ we get
\begin{align*}
I =& \int_{0}^{\frac\pi 2} \theta^4 \ln(\sin\theta)\ \mathrm d\theta\\
=&\int_{0}^{\frac\pi 2} \theta^4\left(-\ln 2 -\sum_{k=1}^\infty \frac{ \cos(2k \theta)}{k}\right)\ \mathrm d\theta\\
=& -\ln 2\int_0^{\frac \pi 2}\theta^4\ \mathrm d\theta-\sum_{k=1}^\infty \frac{1}{k}\underbrace{\int_{0}^{\frac\pi 2}\theta^4 \cos(2k \theta) \ \mathrm d\theta}_{\text{IBP}\times 4}\\
=& -\frac{\pi^5\ln 2}{160}-\sum_{k=1}^\infty \frac{1}{k}\cdot\left(-\frac{\pi^3}{8}\frac{(-1)^{k-1}}{k^2} +\frac{3\pi}{4}\frac{(-1)^{k-1}}{k^4}\right)\\
=&-\frac{\pi^5\ln 2}{160}+\frac{\pi^3}8\underbrace{\sum_{k=1}^\infty \frac{(-1)^{k-1}}{k^3}}_{=\frac 3 4 \zeta(3)} - \frac{3\pi}4\underbrace{\sum_{k=1}^\infty \frac{(-1)^{k-1}}{k^5}}_{=\frac{15}{16}\zeta(5)}\\
=&\boxed{-\frac{\pi^5\ln 2}{160}+\frac{3\pi^3}{32}\zeta(3) -\frac{45\pi}{64}\zeta(5).}
\end{align*}
Combining these, we get the equation
$$
3\pi S-\frac{87\pi}{16} \zeta(5) +\frac{\pi^3}{4}\zeta(3) -\frac{\pi^5\ln 2}{160}=-\frac{\pi^5\ln 2}{160} +\frac{3\pi^3}{32}\zeta(3)-\frac{45\pi}{64}\zeta(5),
$$hence it follows
$$
\boxed{S = \frac{101}{64}\zeta(5) -\frac{5\pi^2}{96}\zeta(3).}
$$
Addendum: By considering MacLaurin series of
\begin{align*}
\ln(1-x)\ln(1+x)
=&-\sum_{k=1}^\infty \left(\frac{ H_{2k}}k-\frac{H_k}{k} + \frac1{2k^2}\right)x^{2k}
\end{align*} and
\begin{align*}
\frac{H_k}{k^2} + \frac{H_k^{(2)}}{k} -\frac{\zeta(2)}{k}
=& \frac{\partial }{\partial k}\left[-\frac{H_k}{k}\right]\\
=& \int_0^1 x^{k-1}\ln x\ln(1-x)\ \mathrm dx\\
=&4\int_0^1 x^{2k-1}\ln x \ln(1-x^2)\ \mathrm dx
\end{align*} we have that
\begin{align*}
&\int_{0}^{1}\ln(1-x)\ln(1+x) \frac{\ln x\ln(1-x^2)}x \ \mathrm dx \\&=-\sum_{k=1}^\infty \left(\frac{ H_{2k}}k-\frac{H_k}{k} + \frac1{2k^2}\right)\int_{0}^{1}x^{2k-1} \ln x \ln(1-x^2)\ \mathrm dx \\
&=-\frac 1 4\sum_{k=1}^\infty \left(\frac{ H_{2k}}k-\frac{H_k}{k} + \frac1{2k^2}\right)\left(\frac{H_k}{k^2} + \frac{H_k^{(2)}}{k} -\frac{\zeta(2)}{k}\right).
\end{align*} The integral can be attacked by considering algebraic identity
$$
ab(a+b) = \frac 1 3 (a+b)^3 - \frac {a^3}3 -\frac{b^3}3
$$ with $a=\ln(1-x)$ and $b=\ln(1+x)$, and extant results.
For the sum, after expanding the summand, the only tricky part is
$$
\sum_{k=1}^\infty\frac{H_{2k}H_k}{k^3},
$$ which can be found here. Then, the sum $\sum_{k=1}^\infty \frac{H_{2k}H_k^{(2)}}{4k^2}$ can be evaluated by solving the equation obtained.
Best Answer
Here is the Cornel's reduction of the $\operatorname{arctanh}$ integral to manageable integrals by using clever integrations by parts and exploiting Dilogarithm reflection formula in $(5)$ here http://mathworld.wolfram.com/Dilogarithm.html.
$$I=\frac{1}{16} \underbrace{\int_0^1 \frac{\text{Li}_2(x) \log ^2(x)}{x} \textrm{d}x}_{\text{known}}-\frac{1}{16} \underbrace{\int_0^1 \frac{\text{Li}_2(x) \log ^2(2-x)}{x} \textrm{d}x}_{\text{may use series with powers of 2 in denominator}}\\-\frac{1}{4} \underbrace{\int_0^1 \frac{\text{Li}_2(x) \log ^2(x)}{2-x} \textrm{d}x}_{\text{may use series with powers of 2 in denominator}}-\frac{3}{8} \underbrace{\int_0^1 \frac{\text{Li}_2(x) \log (x) \log (2-x)}{2-x} \textrm{d}x}_{\text{may use series with powers of 2 in denominator}}\\+\frac{5}{64} \underbrace{\int_0^1 \frac{\log ^4(1-x)}{x} \textrm{d}x}_{\text{known}}-\frac{1}{12} \underbrace{\int_0^1 \frac{\log (1-x) \log ^3(1+x)}{x} \textrm{d}x}_{\text{known}}\\-\frac{7}{16} \underbrace{\int_0^1 \frac{\log^2(1-x) \log (x) \log (1+x)}{1+x} \textrm{d}x}_{\text{Calculated in (Almost) Impossible Integrals, Sums, and Series}}-\frac{1}{4}\underbrace{\int_0^1\frac{\log^3(x)\log(1-x)}{2-x}\textrm{d}x}_{\text{may use series with powers of 2 in denominator}}\\-\frac{7}{32} \underbrace{\int_0^1 \frac{\log ^2(1-x) \log ^2(1+x)}{x} \textrm{d}x}_{\text{known}}.$$
Some of these integrals are already known. The harder ones can be easily reduced to advanced alternating harmonic series of weight $5$ or to advanced harmonic series of weight $5$ with powers of $2$ in denominator (and these ones are all already known). Also, playing with these integrals it's possible to get even more shortcuts. With this form of the main integral to calculate you're done.
A NOTE: The integral $\displaystyle \int_0^1 \frac{\log^2(1-x) \log (x) \log (1+x)}{1+x} \textrm{d}x$ may be found calculated in the book, (Almost) Impossible Integrals, Sums, and Series, pages $527$-$528$. In general, the integrals above where it is suggested a reduction to the series with powers of $2$ in denominator may also be rearranged and reduced to alternating harmonic series of weight $5$. End of story.