After knowing that $$ I_2=\int_0^{+\infty} \frac{\ln ^3 x}{(1+x)^{2}} d x \stackrel{x\mapsto\frac{1}{x}}{=} -\int_0^{+\infty} \frac{\ln ^3 x}{(1+x)^{2}} d x \Rightarrow \int_0^{+\infty} \frac{\ln ^3 x}{(1+x)^{2}} d x=0 $$ I start to investigate the integrals with higher powers $$I_n=\int_0^{+\infty} \frac{\ln ^3 x}{(1+x)^{n}} d x $$ where $n$ is a natural number greater than $2$.
We first split the interval into two as $$ \begin{aligned} \int_0^{+\infty} \frac{\ln ^3 x}{(1+x)^3} d x & =\int_0^1 \frac{\ln ^3 x}{(1+x)^3} d x+\int_1^{+\infty} \frac{\ln ^3 x}{(1+x)^3} d x \\ & =\int_0^1 \frac{\ln ^3 x}{(1+x)^3} d x+\int_1^0 \frac{\ln ^3\left(\frac{1}{x}\right)}{\left(1+\frac{1}{x}\right)^3} \frac{d x}{-x^2} \\ & =\int_0^1 \frac{\ln ^3 x}{(1+x)^3} d x-\int_0^1 \frac{x \ln ^3 x}{(1+x)^3} d x \end{aligned} $$ For any $|x|<1$, we have $$ \frac{1}{1+x}=\sum_{k \rightarrow 0}^{\infty}(-1)^k x^k $$ Differentiating both sides w.r.t. $x$ twice yields $$ \frac{1}{(1+x)^3}=\frac{1}{2} \sum_{k=0}^{\infty}(-1)^k(k+2)(k+1) x^k $$
Plugging into the integrand, we have $$ \int_0^{\infty} \frac{\ln ^3 x}{(1+x)^3}dx=\frac{1}{2}\left[\sum_{k=0}^{\infty}(-1)^k(k+2)(k+1)\left( \int_0^1 x^k \ln ^3 x d x-\int_0^1 x^{k+1} \ln ^3 x d x \right)\right] $$ Noting that $$ \begin{aligned} \int_0^1 x^n \ln ^3 x d x & =\left.\frac{\partial^3}{\partial a^3} \int_0^1 x^a d x\right|_{x=w} \\ & =\left.\frac{\partial^3}{\partial a^3}\left(\frac{1}{a+1}\right)\right|_{x=n} \\ & =-\frac{6}{(n+1)^4} \end{aligned} $$ $$ \begin{aligned}I_3&=\frac{1}{2} \sum_{k=0}^{\infty}(-1)^k(k+2)(k+1)\left[-\frac{6}{(k+1)^4}+\frac{6}{(k+2)^4}\right]\\&= 3 \left[-\sum_{k=0}^{\infty} \frac{(-1)^k(k+2)}{(k+1)^3}+\sum_{k=0}^{\infty} \frac{(-1)^k(k+1)}{(k+2)^3}\right]\\&= -3\left[\sum_{k=0}^{\infty} \frac{(-1)^k(k+2)}{(k+1)^3}+\sum_{k=1}^{\infty} \frac{(-1)^kk}{(k+1)^3}\right]\\&=-3 \sum_{k=0}^{\infty} \frac{(-1)^k(2 k+2)}{(k+1)^3}\\&=-6 \sum_{k=0}^{\infty} \frac{(-1)^k}{(k+1)^2}\\&=-\frac{\pi^2}{2}\end{aligned} $$
Let’s continue with $I_4$ by considering their difference $$ \begin{aligned} D & =\int_0^{+\infty} \frac{\ln ^3 x}{(1+x)^3} d x -\int_0^{+\infty} \frac{\ln ^3 x}{(1+x)^4} d x \\ & =\int_0^{+\infty} \frac{x \ln ^3 x}{(1+x)^4} d x \\ & =\int_0^{+\infty} \frac{\frac{1}{x} \ln ^3 \frac{1}{x}}{\left(1+\frac{1}{x}\right)^4} \frac{d x}{x^2} \\ & =-\int_0^{+\infty} \frac{x \ln ^3 x}{(1+x)^4} d x \\ & =-D\\\Rightarrow D&=0 \end{aligned} $$ Hence $$\boxed{I_4=I_3=-\frac{\pi^2}{2} }$$
By Wolfram-alpha, we get $$ \begin{aligned} & I_5=-\frac{1}{24}\left(6+11 \pi^2\right) \\ & I_6=-\frac{1}{12}\left(6+5 \pi^2\right) \\ & I_7=-\frac{1}{360}(255+137 \pi^2) \end{aligned} $$
I guess that in general, $$ I_n=A+B \pi^2 \textrm{ for some rational numbers } A \textrm{ and }B. $$
Can we prove it further? Your help or suggestions are highly appreciated.
By differentiation
We may start our journey with the integral $$ I(a)=\int_0^{+\infty} \frac{\ln ^3 x}{(a+x)^2} d x $$ Using substitution $ax\mapsto x$ yields $$ \begin{aligned} I(a) & =\frac{1}{a} \int_0^{+\infty} \frac{\ln ^3(a x)}{(1+x)^2} d x \\ & =\frac{1}{a} \int_0^{\infty} \frac{\ln ^3 a+3 \ln ^2 a\ln x+3 \ln a \ln ^2 x+\ln ^3 x}{(1+x)^2} dx \end{aligned} $$ Noting that $\int_0^{+\infty} \frac{\ln ^{2n+1} x}{(1+x)^2} d x=0$, we have $$ I(a)=\frac{\ln ^3 a}{a} \underbrace{ \int_0^{\infty} \frac{1}{\left(1+x\right)^2}d x}_{=1} + \frac{3 \ln a}{a} \underbrace{ \int_0^{\infty} \frac{\ln ^2 x}{(1+x)^2} d x}_{\frac{\pi^2}{3} } =\frac{\ln ^3 a}{a}+\frac{\pi^2 \ln a}{a} $$ By differentiating $I(a)$ w.r.t. $a$ by $n-2$ times, we get
$$\boxed{I_n=\left.\frac{(-1)^n}{(n-1) !} \frac{\partial^{n-2}}{\partial a^{n-2}} \left(\frac{\ln ^3 a}{a}+\frac{\pi^2 \ln a}{a} \right)\right|_{a=1}=A_n+B_n\pi^2 } $$ where $ \left. \displaystyle A_n= \frac{(-1)^n}{4(n-1) !} \frac{\partial^{n-1}}{\partial a^{n-1}} \left(\ln ^4 a\right) \right|_{a=1} $ and $\displaystyle \left. B_n= \frac{(-1)^n}{2(n-1) !} \frac{\partial^{n-1}}{\partial a^{n-1}} \left(\ln^2 a\right) \right|_{a=1} $ are rational.
For examples, $$ \left. I_5=-\frac{1}{24} \frac{\partial^4}{\partial a^4}\left(\frac{\ln ^4 a}{4}+\frac{ \pi^2 \ln^2 a}{2} \right)\right|_{a=1} =-\frac{1}{24}\left(6+11 \pi^2\right) $$ $$I_6=\left.\frac{1}{120} \frac{\partial^5}{\partial a^5}\left(\frac{\ln ^4 a}{4}+\frac{\pi^2 \ln^2 a}{2}\right)\right|_{a=1}= -\frac{1}{12}\left(6+5 \pi^2\right)$$