A closed form for a triple integral with sines and cosines $$\small\int^\infty_0 \int^\infty_0 \int^\infty_0 \frac{\sin(x)\sin(y)\sin(z)}{xyz(x+y+z)}(\sin(x)\cos(y)\cos(z) + \sin(y)\cos(z)\cos(x) + \sin(z)\cos(x)\cos(y))\,dx\,dy\,dz$$
I saw this integral $I$ posted on a page on Facebook . The author claims that there is a closed form for it. 

My Attempt
This can be rewritten as 
$$3\small\int^\infty_0 \int^\infty_0 \int^\infty_0 \frac{\sin^2(x)\sin(y)\cos(y)\sin(z)\cos(z)}{xyz(x+y+z)}\,dx\,dy\,dz$$
Now consider 
$$F(a) = 3\int^\infty_0 \int^\infty_0 \int^\infty_0\frac{\sin^2(x)\sin(y)\cos(y)\sin(z)\cos(z) e^{-a(x+y+z)}}{xyz(x+y+z)}\,dx\,dy\,dz$$
Taking the derivative 
$$F'(a) = -3\int^\infty_0 \int^\infty_0 \int^\infty_0\frac{\sin^2(x)\sin(y)\cos(y)\sin(z)\cos(z) e^{-a(x+y+z)}}{xyz}\,dx\,dy\,dz$$
By symmetry we have
$$F'(a) = -3\left(\int^\infty_0 \frac{\sin^2(x)e^{-ax}}{x}\,dx \right)\left( \int^\infty_0 \frac{\sin(x)\cos(x)e^{-ax}}{x}\,dx\right)^2$$
Using W|A I got 
$$F'(a) = -\frac{3}{16} \log\left(\frac{4}{a^2}+1 \right)\arctan^2\left(\frac{2}{a}\right)$$
By integeration we have 
$$F(0) = \frac{3}{16} \int^\infty_0\log\left(\frac{4}{a^2}+1 \right)\arctan^2\left(\frac{2}{a}\right)\,da$$
Let $x = 2/a$
$$\tag{1}I = \frac{3}{8} \int^\infty_0\frac{\log\left(x^2+1 \right)\arctan^2\left(x\right)}{x^2}\,dx$$
Question
I seem not be able to verify (1) is correct nor find a closed form for it, any ideas ?

$\newcommand{\bbx}[1]{\,\bbox[15px,border:1px groove navy]{\displaystyle{#1}}\,}
 \newcommand{\braces}[1]{\left\lbrace\,{#1}\,\right\rbrace}
 \newcommand{\bracks}[1]{\left\lbrack\,{#1}\,\right\rbrack}
 \newcommand{\dd}{\mathrm{d}}
 \newcommand{\ds}[1]{\displaystyle{#1}}
 \newcommand{\expo}[1]{\,\mathrm{e}^{#1}\,}
 \newcommand{\ic}{\mathrm{i}}
 \newcommand{\mc}[1]{\mathcal{#1}}
 \newcommand{\mrm}[1]{\mathrm{#1}}
 \newcommand{\pars}[1]{\left(\,{#1}\,\right)}
 \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}}
 \newcommand{\root}[2][]{\,\sqrt[#1]{\,{#2}\,}\,}
 \newcommand{\totald}[3][]{\frac{\mathrm{d}^{#1} #2}{\mathrm{d} #3^{#1}}}
 \newcommand{\verts}[1]{\left\vert\,{#1}\,\right\vert}$
Given the 'cyclic symmetry' of your integral $\ds{\color{#f00}{\mc{J}}}$, it's equivalent to
\begin{align}
\color{#f00}{\mc{J}} & \equiv 3\int_{0}^{\infty}\!\!\int_{0}^{\infty}\!\!\int_{0}^{\infty}\!\!
\mrm{sinc}\pars{x}\mrm{sinc}\pars{y}\mrm{sinc}\pars{z}\sin\pars{x}\cos\pars{y}
\cos\pars{z}\ \times
\\[3mm] & \phantom{\equiv 3}
\underbrace{\bracks{\int_{0}^{\infty}\expo{-\pars{x + y + z}t}\,\dd t}}
_{\ds{1 \over x + y + z}}\dd x\,\dd y\,\dd z =
3\int_{0}^{\infty}\mrm{f}\pars{t}\mrm{g}^{2}\pars{t}\,\dd t
\\[5mm] & \mbox{where}\qquad
\left\{\begin{array}{rcl}
\ds{\mrm{f}\pars{t}} & \ds{=} &
\ds{\Im\mc{I}\pars{t}}
\\[2mm]
\ds{\mrm{g}\pars{t}} & \ds{=} &
\ds{\Re\mc{I}\pars{t}}
\\[2mm]
\ds{\mc{I}\pars{t}} & \ds{\equiv} &
\ds{\int_{0}^{\infty}\mrm{sinc}\pars{x}\expo{-\pars{t - \ic}x}\dd x}
\end{array}\right.\label{1}\tag{1}
\end{align}

Then
\begin{align}
\left.\mc{I}\pars{t}\right\vert_{\large\ \color{#f00}{t\ >\ 0}} & \equiv
\int_{0}^{\infty}\mrm{sinc}\pars{x}\expo{-\pars{t - \ic}x}\dd x =
\int_{0}^{\infty}
\overbrace{\pars{{1 \over 2}\int_{-1}^{1}\expo{\ic kx}\,\dd k}}
^{\ds{\mrm{sinc}\pars{x}}}\ \expo{-\pars{t - \ic}x}\dd x
\\[5mm] & =
{1 \over 2}\int_{-1}^{1}\int_{0}^{\infty}\expo{-\pars{t - \ic - \ic k}x}
\dd x\,\dd k =
{1 \over 2}\int_{-1}^{1}{\dd k \over t - \ic - \ic k} =
{1 \over 2}\int_{-1}^{1}{t + \pars{k + 1}\ic \over
\pars{k + 1}^{2} + t^{2}}\,\dd k
\\[5mm] & =
{1 \over 2}\int_{0}^{2}{t + k\ic \over k^{2} + t^{2}}\,\dd k =
{1 \over 2}\int_{0}^{2/t}{1 + k\ic \over k^{2} + 1}\,\dd k =
{1 \over 2}\arctan\pars{2 \over t} +
{1 \over 4}\ln\pars{{4 \over t^{2}} + 1 }\ic
\end{align}

$\ds{\color{#f00}{\mc{J}}}$ becomes ( see \eqref{1} ):
\begin{align}
\color{#f00}{\mc{J}} & =
3\int_{0}^{\infty}\bracks{{1 \over 4}\ln\pars{{4 \over t^{2}} + 1}}
\bracks{{1 \over 2}\arctan\pars{2 \over t}}^{2}\dd t
\\[5mm] & \stackrel{2/t\ \mapsto\ t}{=}\,\,\,
{3 \over 16}\int_{\infty}^{0}\ln\pars{t^{2} + 1}
\arctan^{2}\pars{t}\pars{-\,{2\,\dd t \over t^{2}}} =
{3 \over 8}\ \underbrace{\int_{0}^{\infty}
{\ln\pars{t^{2} + 1}\arctan^{2}\pars{t} \over t^{2}}\,\dd t}
_{\ds{{\Large\color{#f00}{\S}}: {\pi^{3} \over 12} + \pi\ln^{2}\pars{2}}}
\end{align}


$\ds{{\Large\color{#f00}{\S}}}$: The integral was already evaluated in the
  $\texttt{@Zaid Alyafeai}$ fine answer.


Finally, the answer to the proposed OP integral is given by
$$
\bbox[15px,#ffe,border:1px dotted navy]{\ds{{\color{#f00}{\mc{J}} =
{\pi^{3} \over 32} + {3 \over 8}\,\pi\ln^{2}\pars{2}}}} \approx 1.5350
$$
