In my answer to this question, I come across the following case of the Meijer G-function:
$$F(b)=G^{2~2}_{3~3}\left(1\middle|\begin{array}c1,1;b+1\\b,b;0\end{array}\right), b>0$$
and based on my experiments, $F(b)$ have the following closed form:
$$F(b) \stackrel?=\frac{\Gamma(b)}{b}\left(-\gamma-\psi(b)+\frac{2^{1-b}}{b}{_2F_1}\left(\begin{array}c1,1\\b+1\end{array}\middle|-1\right)+b{_3F_2}\left(\begin{array}c1,1,b+1\\2,2\end{array}\middle|-1\right)\right)$$
Is there any chance of proving this?
Edit: Using the definition of the Meijer G-Function, $F(b)$ have the integral representation: $$F(b)=\frac{1}{2\pi}\int^{+\infty}_{-\infty}\frac{\Gamma(\tfrac{b}{2}+ix)\Gamma(\tfrac{b}{2}-ix)}{\tfrac{b^2}{4}+x^2}dx.$$
Edit 2: I've found a further generalization: $$F(b,z)=G^{2~2}_{3~3}\left(z\middle|\begin{array}c1,1;b+1\\b,b;0\end{array}\right)=\frac{1}{2\pi}\int^{+\infty}_{-\infty}\frac{\Gamma(\tfrac{b}{2}+ix)\Gamma(\tfrac{b}{2}-ix)}{\tfrac{b^2}{4}+x^2}e^{(b/2+ix)\log z}dx\\ \stackrel?=\frac{\Gamma(b)z^b}{b}\left(-\log z-\gamma-\psi(b)+\frac{(z+1)^{1-b}}{b}{_2F_1}\left(\begin{array}c1,1\\b+1\end{array}\middle|-z\right)+bz~{_3F_2}\left(\begin{array}c1,1,b+1\\2,2\end{array}\middle|-z\right)\right), b\not\in\{0,-1,-2,\dots\},z\neq0,-1.$$ Edit 3: Further simplified the ${_2F_1}$ part.
Case $F(1,1)$ is this claim: $$ 2 \int_{-\infty}^\infty \frac{dx}{(4 x^2+1)\cosh(\pi x)} = 2\log 2 . $$ Maybe start with a proof of that.
added March 24
We will attempt this $F(1,1)$ case using residues. I changed variables slightly to make the problem $$ \int_{-\infty}^{+\infty} \frac{dx}{(x^2+1)\cosh(\pi x/2)} = 2\log 2 . $$ The function $$ F(z) := \frac{1}{(z^2+1)\cosh(\pi z/2)} $$ is analytic in the complex plane except for poles at $(2n+1)i$, $n \in \mathbb Z$. Let $M$ be a large integer, and consider the contour integral of $F$ around the square with vertices $$ -4 M ,\quad 4 M ,\quad 4M + 8 M i,\quad -4 M\pi + 8 M i . $$ The integral along the bottom side is $$ \int_{-4M}^{4M} \frac{dx}{(x^2+1)\cosh(\pi x/2)} $$ so the limit as $M \to \infty$ is the integral to be computed. We can do estimates for the other three sides to show their integrals go to zero as $M \to \infty$. (Postponed...)
The poles for $F(z)$ inside the square are at $z=(2n+1) i$, for $n=0,1,\dots, 4M-1$. Now $z=i$ is a double pole, the residue there is $$ \mathrm{Res}_{z=i} F(z) = \frac{-i}{2\pi} . $$ For $n>0$, the residue is $$ \mathrm{Res}_{z=(2n+1) i} F(z) = \frac{(-1)^n i}{2 n (n+1)\pi} $$ Now as $M \to \infty$, the square encloses more and more of these poles, so in the limit the sum is the sum of all the residues in the upper half-plane: $$ \frac{-i}{2\pi}+\sum_{n=1}^\infty \frac{(-1)^n i}{2n(n+1)\pi} =\frac{-i\log 2}{\pi} $$ so finally the integral around the contour (which converges to the integral we want) also converges to $2\pi i$ times this sum of residues, $$ \int_{-\infty}^{+\infty} F(z)\,dz = 2 \log 2 . $$