Aunque he visto muy pocas pruebas en este campo para poder comparar, este enfoque podría ser muy interesante.
Podemos transformar la suma a una cuádruple integral que Mathematica puede resolver de inmediato. Espero que debería ser posible para resolver la integral "matemáticamente" así, que luego de completar la prueba.
Tenemos que calcular
$$s=\sum _{k=2}^{\infty } \frac{(-1)^k}{k^2} H(k) H(k-1) $$
Escrito
$$\frac{1}{n^2}=\int_0^1 \frac{1}{x}\,dx \int_0^x y^{n-1} \, dy $$
$$\frac{1}{n}=\int_0^1 r^{n-1} \, dr$$
y
$$H(k)=\sum _{n=1}^k \frac{1}{n}=\int_0^1 \left(\sum _{n=1}^k r^{n-1}\right) \, dr=\int_0^1 \frac{1-r^k}{1-r} \, dr$$
la suma de $s$ por debajo de la integral se convierte en
$$si=\frac{1}{x(1-r)(1-s)}\sum _{k=2}^{\infty } (-1)^k \left(1-r^k\right) \left(1-s^{k-1}\right) y^{k-1}$$
En el que se evalúa a
$$si = \frac{y \left(r^2 s^2 y+r^2 s-r^2 y-r^2-r s^2 y+r y-s+1\right)}{(1-r) (1-s) x (y+1) (r y+1) (s y+1) (r s y+1)}$$
Ahora la integral que se evalúa es
$$s4 = \int _0^1 dx\int _0^x dy\int _0^1 dr\int _0^1 ds \; si$$
Mathematica encuentra inmediatamente
$$s4 = \frac{\pi ^4}{480} $$
Desde
$$\zeta (4)=\frac{\pi ^4}{90}$$
y
$$\frac{90}{480} = \frac{3}{16} $$
por fin hemos
$$s = \frac{3}{16} \zeta(4)$$