Hay muchos buenos identidades jugando aquí:
$$\sum_{j+k=u}\frac{1}{jk}=\frac{2H_{u-1}}{u}\etiqueta{1},$$
$$\sum_{h=1}^{n}\frac{H_h}{h}=\frac{1}{2}\left(H_n^2+H_n^{(2)}\right)\etiqueta{2},$$
$$\sum_{h=1}^{n}\frac{H_{h-1}}{h}=\frac{1}{2}\left(H_n^2-H_n^{(2)}\right)\etiqueta{3},$$
$$\sum_{j=1}^{n}\frac{H_j}{n+1-j}=H_{n+1}^2-H_{n+1}^{(2)}\la etiqueta{4}.$$
Pruebas:
$$(1)\quad \sum_{j+k=u}\frac{1}{jk}=\sum_{j=1}^{u-1}\frac{1}{j(u-j)}=\frac{1}{u}\sum_{j=1}^{u-1}\left(\frac{1}{j}+\frac{1}{u-j}\right)=\frac{2H_{u-1}}{u}.$$
$$(2,3)\quad H_n^2=H_n^{(2)}+2\sum_{j=1}^{n}\frac{1}{j}\sum_{k<j}\frac{1}{k}=H_n^{(2)}+2\sum_{j=1}^{n}\frac{H_{j-1}}{j}=-H_n^{(2)}+2\sum_{j=1}^{n}\frac{H_{j}}{j}.$$
Mediante la definición de $H_0=0$, tenemos:
$$(4)\quad\sum_{j=1}^{n}\frac{H_j}{n+1-j}=[x^{n+1}]\frac{\log^2(1-x)}{1-x},$$
pero sabemos que la serie de coeficientes de $\log^2(1-x)$, entonces:
$$(4)\quad\sum_{j=1}^{n}\frac{H_j}{n+1-j}=\sum_{k=1}^{n+1}\frac{2H_{k-1}}{k}$$
y $(4)$ de la siguiente manera a partir de $(3)$. Ahora vamos a jugar un poco con la última identidad que Greg Martin demostró:
$$I_n = \frac{1}{2}\sum_{i+j+k\leq n+3}\frac{1}{ijk}=\frac{1}{2}\sum_{t=3}^{n+3}\sum_{i+j+k=t}\frac{1}{ijk}=\frac{1}{2}\sum_{h=3}^{n+3}\sum_{i=1}^{t-2}\frac{1}{i}\sum_{j+k=t-i}\frac{1}{jk};$$
el uso de $(1)$ tenemos:
$$I_n=\sum_{t=3}^{n+3}\sum_{i=1}^{t-2}\frac{H_{t-i-1}}{i(t-i)}=\sum_{t=3}^{n+3}\frac{1}{t}\sum_{i=1}^{t-2}H_{t-i-1}\left(\frac{1}{i}+\frac{1}{t-i}\right),$$
a continuación, volver a indexar el interior de la suma:
$$I_n=\sum_{t=3}^{n+3}\frac{1}{t}\sum_{i=1}^{t-2}H_{i}\left(\frac{1}{t-i-1}+\frac{1}{i+1}\right),$$
donde:
$$\sum_{i=1}^{t-2}\frac{H_i}{i+1}=\sum_{i=1}^{t-1}\frac{H_{i-1}}{i}=\frac{1}{2}\left(H_{t-1}^2-H_{t-1}^{(2)}\right)$$
por $(3)$ y
$$\sum_{i=1}^{t-2}\frac{H_i}{t-1-i}=H_{t-1}^2-H_{t-1}^{(2)}$$
por $(4)$. Poniendo todo junto, se obtiene:
$$I_n=\frac{3}{2}\sum_{t=3}^{n+3}\frac{H_{t-1}^2-H_{t-1}^{(2)}}{t}.\la etiqueta{5}$$
Esto le da a ese $I_n$ se comporta como $\frac{1}{2}\log^3$n.
La parte de debajo de esta línea es superado por mi próxima respuesta.
Por otra parte, la sumación por partes, se obtiene:
$$\sum_{k=1}^{m+1}\frac{H_{k-1}^{(2)}}{k}=H_{m}^{(2)}H_{m+1}-\sum_{j=1}^m \frac{H_j}{j^2},$$
$$\sum_{k=1}^{m+1}\frac{H_{k-1}^{2}}{k}=H_{m}^{2}H_{m+1}-\sum_{j=1}^m \frac{(H_j+H_{j-1})H_j}{j}=H_{m}^{2}H_{m+1}-\sum_{j=1}^m \frac{(2H_{j-1}+1/j)(H_{j-1}+1/j)}{j},$$
así:
$$\sum_{k=1}^{m+1}\frac{H_{k-1}^{2}}{k}=H_{m}^{2}H_{m+1}-2\sum_{j=1}^{m}\frac{H_{j-1}^2}{j}-3\sum_{j=1}^m\frac{H_{j-1}}{j^2}-H_{m}^{(3)},$$
y reordenando, se obtiene:
$$\sum_{k=1}^{m+1}\frac{H_{k-1}^{2}}{k}=\frac{1}{3}H_{m}^{2}H_{m+1}+\frac{2H_m^2}{3(m+1)}-\sum_{j=1}^m\frac{H_{j-1}}{j^2}-\frac{1}{3}H_m^{(3)}.$$
Sumando todo, tenemos que $(5)$ puede ser escrita en la forma siguiente:
$$I_{n}=\frac{1}{2}H_{n+2}^2 H_{n+3}+\frac{H_{n+2}^2}{n+3}-\frac{3}{2}H_{n+2}^{(2)}H_{n+3}+H_{n+2}^{(3)}+\frac{3}{2}\sum_{j=1}^{n+2}\frac{H_j}{j^2}.\tag{6}$$
La última suma es ahora claramente delimitado por una absoluta constante ($\zeta(3)$, para istance), y creo firmemente que no simplificar aún más (me equivoqué).
