Como se indicó en los comentarios, por los siguientes Feynman del enfoque integral tenemos:
$$\begin{eqnarray*} I'(s) &=& \int_{0}^{1}\left(1+\frac{\log x}{1-x}-\frac{\log x}{2}\right)\frac{x^s}{1-x}dx\\ &=& \frac{1}{2}\zeta(2,s+1)+\int_{0}^{1}\left(1+\frac{\log x}{1-x}\right)\frac{x^s}{1-x}dx\\&=&\frac{1}{2}\zeta(2,s+1)-\frac{1}{2s+2}\phantom{}_3 F_2\left(1,1,2;3,s+2;1\right),\end{eqnarray*}$$
pero probablemente no hay una mejor forma para el segundo plazo.
Como cuestión de hecho, al tratar el segundo término de esta manera:
$$\int_{0}^{1}\left(1+\frac{\log x}{1-x}\right)\frac{x^s}{1-x}dx=\lim_{a\to 0^+}\int_{0}^{1}\left(1+\frac{\log x}{1-x}\right)\frac{x^s}{(1-x)^{1-a}}$$
obtenemos:
$$ I'(s)=\frac{1}{2}\psi'(s+1)+s\,\psi'(s)+\psi(s)-\psi(s+1)-1\tag{1}$$
donde $\psi$ es la función digamma $\psi(s)=\frac{\Gamma'(s)}{\Gamma(s)}$. Mediante la integración de $(1)$ tenemos:
$$ I(s) =
I(1)+\frac{3\gamma+1}{2}-s+\frac{\psi(s+1)}{2}-\log\Gamma(s+1)+s\,\psi(s).\la etiqueta{2}$$
Todos necesitamos encontrar un cerrado expresión para $I(s)$ es ahora a evaluar $I(1)$ o $I(0)$, dado por $(2)$:
$$ I(0)=I(1)+\gamma-\frac{1}{2}$$
sostiene. En virtud de $(2)$ también tenemos:
$$\lim_{s\to +\infty} I(s) = I(1) + \frac{3\gamma-\log(2\pi)}{2},$$
pero ya
$$ f(x)=\left(\frac{1}{\log x}+\frac{1}{1-x}-\frac{1}{2}\right)\frac{x}{1-x}$$
es un continuo, de manera positiva, creciente y acotada de la función ( $f(x)\leq\frac{1}{12}$ )$(0,1)$, mientras $s\to+\infty$ tenemos $I(s)\leq\frac{1}{12s}$, lo $\lim_{s\to +\infty}I(s)=0$. Por lo tanto tenemos:
$$ I(1) = \frac{\log(2\pi)-3\gamma}{2},\qquad I(0) = \frac{\log(2\pi)-\gamma-1}{2} \tag{3}$$
y:
$$ I(s) =
\frac{1+\log(2\pi}{2}-s+\frac{\psi(s+1)}{2}-\log\Gamma(s+1)+s\,\psi(s).\la etiqueta{4}$$
Como un bono, por la expansión de $(4)$ en un barrio de $+\infty$ y la explotación de la convexidad, también se puede mostrar que: $$\frac{1}{12s+8}\leq I(s)\leq\frac{1}{12s+6}$$ holds for any $s\in\mathbb{R}^+$.