<s>[INST] Prove $\sum_{k=1}^{\infty}\frac{(-1)^k}{(2k+1)^2}(2H_{2k}+H_k)=\frac{\pi^3}{32}-2G\ln2$ How to prove

$$\sum_{k=1}^{\infty}\frac{(-1)^k}{(2k+1)^2}(2H_{2k}+H_k)\stackrel ?=\frac{\pi^3}{32}-2G\ln2,$$
  where $G$ is the Catalan's constant.

Attempt
For the first sum,
$$\sum_{k=1}^{\infty}\frac{(-1)^k}{(2k+1)^2}H_{2k}=\Re\left\{\sum_{k=1}^{\infty}\frac{i^k}{(k+1)^2}H_{k}\right\},$$
which can be evaluated by using the formula in this post:
$$\sum_{n=1}^\infty \frac{H_n}{n^2}\, x^n=\zeta(3)+\frac{\ln(1-x)^2\ln(x)}{2}+\ln(1-x)\operatorname{Li}_2(1-x)+\operatorname{Li}_3(x)-\operatorname{Li}_3(1-x),$$
but we cannot apply the similar approach to the second sum 
$$\sum_{k=1}^{\infty}\frac{(-1)^k}{(2k+1)^2}H_k.$$
Then, I tried to write the sum as
$$\sum_{k=1}^{\infty}\frac{(-1)^k}{(2k+1)^2}\int_0^1\frac{2x^{2k}+x^k-3}{x-1}~\mathrm dx$$
and it become more complicated.
Edit: 
Are we able to evaluate the sum directly (avoid calculating integrals and  polylogs as much as possible)? The integral given by @Jack D'Aurizio is a bit complicated (see this post).
 [/INST] this summation also has some relationship with integral identity, which asked almost the same time, see this, rewrite as
$$\frac3{2} \int_{0}^{1} {\frac{\arctan x \ln(1+x^2)}{x} \mathrm{d}x} = \int_{0}^{1} {\frac{\arctan x \ln(1+x)}{x} \mathrm{d}x}$$
given some famous series
$$\arctan x = \sum_{n=0}^{\infty} {\frac{(-1)^n x^{2n+1}}{2n+1}}$$
$$\frac1{2} \arctan x \ln (1+x^2) = \sum_{n=0}^{\infty} {\frac{(-1)^{n+1}H_{2n}}{2n+1} x^{2n+1}}$$
on the left
$$\frac3{2} \int_{0}^{1} {\frac{\arctan x \ln(1+x^2)}{x} \mathrm{d}x} = 3\int_{0}^{1} {\sum_{n=0}^{\infty} {\frac{(-1)^{n+1}H_{2n}}{2n+1} x^{2n}} \>\mathrm{d}x} = 3\sum_{n=0}^{\infty} {\frac{(-1)^{n+1}H_{2n}}{(2n+1)^2}}$$
on the right
$$\begin{aligned}
\int_{0}^{1} {\frac{\arctan x \ln(1+x)}{x} \mathrm{d}x}
& = \int_{0}^{1} {\sum_{n=0}^{\infty} {\frac{(-1)^n}{2n+1}} x^{2n}\ln(1+x) \>\mathrm{d}x}\\
& = 2\ln2\sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^2}} - \sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^2} \int_{0}^{1} {\frac{1+x^{2n+1}}{1+x}}\mathrm{d}x}\\
& = 2G\ln2 - \sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^2} \int_{0}^{1} {\sum_{k=0}^{2n} {(-x)^k}}\mathrm{d}x}\\
& = 2G\ln2 + \sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^2} \sum_{k=1}^{2n+1} {\frac{(-1)^k}{k}}}\\
& = 2G\ln2 + \sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^2}(H_{n}-H_{2n+1})}
\end{aligned}$$
unscramble this identity
$$\sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^2}(H_{n}+3H_{2n}-H_{2n+1})} = -2G\ln2$$
or
$$\begin{aligned}
\sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^2}(H_{n}+2H_{2n})} 
& = \sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^2}(H_{2n+1}-H_{2n})} - 2G\ln2\\
& = \sum_{n=0}^{\infty} {\frac{(-1)^n}{(2n+1)^3}} - 2G\ln2 = \frac{\pi^3}{32} - 2G\ln2
\end{aligned}$$
where you may need some particular values of Dirichlet Beta Function.
</s>