Un lugar agradable expresión para la suma se puede encontrar utilizando la serie de Fourier
$$ \sum \limits_{n=1}^\infty \frac{\sin(n y)}{n} = \frac{\pi}{2} \left[1 - 2 \left\{\frac{y}{2\pi}\right\} \right] \, , \, y \in \mathbb{R} \, , $$
donde $\{z\} = z - \lfloor z \rfloor$ es la parte fraccionaria de $z \in \mathbb{R}$ . Con esta definición satisface $\{-z\} = 1- \{z\}$ para $z \in \mathbb{R} \setminus \mathbb{Z}$ .
Podemos escribir la suma como
\begin{align}
S &\equiv - \sum \limits_{n=1}^\infty \frac{I_n}{n} = \frac{\pi}{2} \int \limits_0^\infty \mathrm{e}^{-x} \left[2 \left\{\frac{\ln(x)}{2\pi}\right\} - 1\right] \, \mathrm{d} x = \frac{\pi}{2} \int \limits_{-\infty}^\infty \mathrm{e}^{-\mathrm{e}^{t}+t} \left[2 \left\{\frac{t}{2\pi}\right\} - 1\right] \, \mathrm{d} t\\
&= \frac{\pi}{2} \int \limits_0^\infty \left(\mathrm{e}^{-\mathrm{e}^{t}+t} \left[2 \left\{\frac{t}{2\pi}\right\} - 1\right] + \mathrm{e}^{-\mathrm{e}^{-t}-t} \left[2 \left(1-\left\{\frac{t}{2\pi}\right\}\right) - 1\right] \right)\, \mathrm{d} t \\
&= \frac{\pi}{2} \int \limits_0^\infty \left(\mathrm{e}^{-\mathrm{e}^{-t}-t} - \mathrm{e}^{-\mathrm{e}^{t}+t} \right) \left(1 - 2\left\{\frac{t}{2\pi}\right\} \right) \, \mathrm{d}t \\
&= \frac{\pi}{2} \left[1 - \frac{1}{\mathrm{e}} - \frac{1}{\mathrm{e}} - 2 \sum \limits_{n=0}^\infty \int \limits_{2 \pi n}^{2 \pi (n+1)} \left(\mathrm{e}^{-\mathrm{e}^{-t}-t} - \mathrm{e}^{-\mathrm{e}^{t}+t} \right) \left(\frac{t}{2\pi} - n\right) \, \mathrm{d} t \right] \\
&= \frac{\pi}{2} \left[1 - \frac{2}{\mathrm{e}} - \frac{1}{\pi} \int \limits_0^\infty \left[-\ln(x) \mathrm{e}^{-x}\right] \, \mathrm{d} x + 2 \sum \limits_{n=0}^\infty n \left(\mathrm{e}^{-\mathrm{e}^{-2\pi(n+1)}}-\mathrm{e}^{-\mathrm{e}^{-2\pi n}} - \mathrm{e}^{-\mathrm{e}^{2\pi n}} + \mathrm{e}^{-\mathrm{e}^{2\pi(n+1)}}\right) \right] \\
&\equiv \frac{\pi}{2}-\frac{\pi}{\mathrm{e}}- \frac{\gamma}{2} + \pi (S_1 - S_2) \, .
\end{align}
La segunda suma es esencialmente igual a cero:
$$ S_2 = \sum \limits_{n=1}^\infty n \left( \mathrm{e}^{-\exp[2\pi n]} - \mathrm{e}^{-\exp[2\pi(n+1)]}\right) = \sum \limits_{n=1}^\infty \mathrm{e}^{- \exp(2 \pi n)} \simeq 3 \cdot 10^{-233} \, . $$
La primera suma puede ser simplificado con el uso de la exponencial de la serie:
$$ S_1 = \sum \limits_{n=1}^\infty n \left( \mathrm{e}^{-\exp[-2\pi (n+1)]} - \mathrm{e}^{-\exp[-2\pi n]}\right) = \sum \limits_{k=1}^\infty \frac{(-1)^{k-1}}{k!} \sum \limits_{n=1}^\infty \mathrm{e}^{-2 \pi k n} = \sum \limits_{k=1}^\infty \frac{(-1)^{k-1}}{k! (\mathrm{e}^{2 \pi k} - 1)} \, .$$
Mantener sólo los tres primeros términos de $S_1$ e ignorando $S_2$ , obtenemos
$$ S \simeq \pi\left[\frac{1}{2} - \frac{1}{\mathrm{e}} - \frac{\gamma}{2 \pi} + \frac{1}{\mathrm{e}^{2\pi} - 1} - \frac{1}{2(\mathrm{e}^{4\pi} - 1)} + \frac{1}{6(\mathrm{e}^{6\pi} - 1)}\right] \simeq 0.13233339071\color{red}{3} \, .$$
El rojo tres es la primera desviación de la 'exacta' valor. No creo que hay expresiones cerradas para las $S_1$ e $S_2$ , pero al menos este resultado no es demasiado lejos.