I am trying to find the closed form expression for $${_4F_3}\!\left(\begin{array}c\tfrac12,\tfrac12,\tfrac12,\tfrac12\\\tfrac32,\tfrac32,\tfrac32\end{array}\middle|\tfrac12\right).$$ I encountered this expression when evaluating the integral $$\int_0^1 \frac{\ln^2(1-x)}{\sqrt{1-x^2}}\,dx.$$ I tried utilizing the series representation of the generalized hypergeometric function Pochhammer symbol and applying Legendre duplication, which left me with the infinite series $$\sum_{n=0}^\infty \frac{\binom{2n}{n}}{8^n(2n+1)^3}.$$ Using Mathematica, it can be approximated to be $${_4F_3}\!\left(\begin{array}c\tfrac12,\tfrac12,\tfrac12,\tfrac12\\\tfrac32,\tfrac32,\tfrac32\end{array}\middle|\tfrac12\right) \approx 1.01015469252627.$$
I have reason to believe that the closed form may contain the inverse tangent integral function, as the hypergeometric function is similar to @Cleo's answer. Are there any neat properties of the $_4F_3$ hypergeometric function that can be exploited in order to obtain the closed form result?