It seems to have escaped attention that this sum may be evaluated
using harmonic summation techniques.
Put $$S(x) = \zeta(3) +
\sum_{n\ge 1} \frac{-1+\coth(n\pi x)}{n^3}$$
and introduce the sum
$$T(x) = \sum_{n\ge 1} \frac{-1+\coth(n\pi x)}{n^3}.$$
The sum term is harmonic and may be evaluated by inverting its Mellin
transform. We will construct a functional equation for $T(x).$
Recall the harmonic sum identity
$$\mathfrak{M}\left(\sum_{k\ge 1} \lambda_k g(\mu_k x);s\right) =
\left(\sum_{k\ge 1} \frac{\lambda_k}{\mu_k^s} \right) g^*(s)$$
where $g^*(s)$ is the Mellin transform of $g(x).$
In the present case we have
$$\lambda_k = \frac{1}{k^3}, \quad \mu_k = k
\quad \text{and} \quad
g(x) = 2\frac{e^{-2\pi x}}{1-e^{-2\pi x}}.$$
We need the Mellin transform $g^*(s)$ of $g(x)$ which is
$$2 \int_0^\infty \frac{e^{-2\pi x}}{1-e^{-2\pi x}} x^{s-1} dx
\\ = 2 \int_0^\infty \sum_{q\ge 1} e^{-2q\pi x} x^{s-1} dx
= 2 \sum_{q\ge 1} \int_0^\infty e^{-2q\pi x} x^{s-1} dx
\\= 2 \Gamma(s) \sum_{q\ge 1} \frac{1}{(2\pi q)^s}
= \frac{2}{2^s} \frac{1}{\pi^s} \Gamma(s) \zeta(s).$$
It follows that the Mellin transform $Q(s)$ of the harmonic sum
$T(x)$ is given by
$$Q(s) =
\frac{2}{2^s} \frac{1}{\pi^s} \Gamma(s) \zeta(s) \zeta(s+3)
\quad\text{because}\quad
\sum_{k\ge 1} \frac{\lambda_k}{\mu_k^s} =
\sum_{k\ge 1} \frac{1}{k^3} \frac{1}{k^s}
= \zeta(s+3)$$
for $\Re(s) > -2.$
The Mellin inversion integral here is
$$\frac{1}{2\pi i} \int_{3/2-i\infty}^{3/2+i\infty} Q(s)/x^s ds$$
which we evaluate by shifting it to the left for an expansion about
zero.
Fortunately the trivial zeros of the two zeta function terms cancel
the poles of the gamma function term. Shifting to $\Re(s) = -3 -1/2$
we get
$$T(x)
= \frac{\pi^3 x^3}{90}
+ 4\zeta'(-2)\pi^2 x^2
+ \frac{\pi^3 x}{18}
- \zeta(3)
+ \frac{\pi^3}{90x}
+ \frac{1}{2\pi i} \int_{-7/2-i\infty}^{-7/2+i\infty} Q(s)/x^s ds.$$
Substitute $s = -2 - t$ in the remainder integral to get
$$- \frac{1}{2\pi i}
\int_{3/2+i\infty}^{3/2-i\infty}
\frac{2}{2^{-2-t}}
\frac{1}{\pi^{-2-t}} \Gamma(-2-t) \zeta(-2-t) \zeta(1-t)
x^{t+2} dt$$
which is
$$\frac{x^2}{2\pi i}
\int_{3/2-i\infty}^{3/2+i\infty}
2^{3+t} \pi^{2+t} \Gamma(-2-t) \zeta(-2-t) \zeta(1-t)
x^t dt.$$
In view of the desired functional equation we now use the functional
equation of the Riemann zeta function on $Q(s)$ to prove that the
integrand of the last integral is in fact $-Q(t).$
Start with the functional equation
$$\zeta(1-s) = \frac{2}{2^s\pi^s}
\cos\left(\frac{\pi s}{2}\right) \Gamma(s) \zeta(s)$$
and substitute this into $Q(s)$ to obtain
$$Q(s) =
\frac{2}{2^s} \frac{1}{\pi^s}
\frac{\zeta(1-s) 2^s \pi^s}{2\cos\left(\frac{\pi s}{2}\right)}
\zeta(s+3)
= \frac{\zeta(3+s)}{\cos\left(\frac{\pi s}{2}\right)}
\zeta(1-s).$$
Apply the functional equation again (this time to $\zeta(s+3)$) to get
$$Q(s) = \frac{1}{\cos\left(\frac{\pi s}{2}\right)}
\frac{2}{2^{-2-s} \pi^{-2-s}}
\cos\left(\frac{\pi (-2-s)}{2}\right)
\Gamma(-2-s) \zeta(-2-s) \zeta(1-s)$$
Observe that
$$\frac{\cos\left(-\pi-\frac{\pi s}{2}\right)}
{\cos\left(\frac{\pi s}{2}\right)}
= - \frac{\cos\left(-\frac{\pi s}{2}\right)}
{\cos\left(\frac{\pi s}{2}\right)} = -1$$ so we finally get
$$Q(s) =
- 2^{3+s} \pi^{2+s} \Gamma(-2-s) \zeta(-2-s) \zeta(1-s),$$
thus proving the claim.
We have established the functional equation
$$T(x)
= \frac{\pi^3 x^3}{90}
+ 4\zeta'(-2)\pi^2 x^2
+ \frac{\pi^3 x}{18}
- \zeta(3)
+ \frac{\pi^3}{90x}
- x^2 T(1/x).$$
Finally returning to the sum that was the initial goal we see that it
has the value $$\zeta(3) + T(x) + x^2 (\zeta(3) + T(1/x))$$
or
$$\zeta(3) + T(x) + x^2 \zeta(3) + x^2 T(1/x).$$
Using the functional equation for $T(x)$ this becomes
$$\zeta(3) + T(x) + x^2 \zeta(3) +
\frac{\pi^3 x^3}{90}
+ 4\zeta'(-2)\pi^2 x^2
+ \frac{\pi^3 x}{18}
- \zeta(3)
+ \frac{\pi^3}{90x}
- T(x)$$ which is
$$x^2 \zeta(3)
+ \frac{\pi^3 x^3}{90}
+ 4\zeta'(-2)\pi^2 x^2
+ \frac{\pi^3 x}{18}
+ \frac{\pi^3}{90x}.$$
The inspiration for this calculation is from the paper "Mellin Transform and its Applications" by Szpankowski.
Addendum.
In view of the fact that $$\zeta(3) + 4\zeta'(-2)\pi^2 =0 $$
(consult e.g. MathWorld)
this finally becomes
$$\frac{\pi^3 x^3}{90}
+ \frac{\pi^3 x}{18}
+ \frac{\pi^3}{90x}
= \frac{\pi^3}{90x}
\left(x^4 + 5x^2 + 1\right).$$
Addendum II. There is another functional equation of a harmonic sum at this MSE link, this one somewhat more advanced.
Best Answer
The first identity is (from my point of view) a really hard one. It was already many times discussed on MSE. Particular interest for you can present this discussion, where also a reference concerning the method probably used by Ramanujan was given.
The others two sums can be evaluated by elementary means of residue calculus. The general method is described here. The same method can be applied whenever a polynomial is present in the denominator of the expression. Let us demonstrate this on the example of the series $$ \sum_{n=0}^{\infty} \frac{(-1)^n}{(2n+1)^{4K+1} \cosh \left(\frac{(2n+1)\pi}{2} \right)} $$
For this we consider the integral: $$ \oint_{\Gamma_\nu}\frac{dz}{z^k\cos z\cosh z} $$ where $\Gamma_\nu$ is a square contour $$ (-\nu,-\nu)\to(\nu,-\nu)\to(\nu,\nu)\to(-\nu,\nu)\to(-\nu,-\nu) $$ avoiding the poles of the denominator situated at $(n+\frac12)\pi$ and $(n+\frac12)\pi i$ ($n\in\mathbb Z$). Exactly the same reasoning as in the cited answer will result in the identity: $$ 4\sum_{n=0}^{\infty} \frac{(-1)^n}{\left(\frac{(2n+1)\pi}{2} \right)^{4K+1} \cosh \left(\frac{(2n+1)\pi}{2} \right)}=A_{4K} $$ with $A_k$ defined by the series expansion $$ \frac1{\cos z\cosh z}=\sum_{k=0}^\infty A_kz^k. $$
As, particularly, $A_4=\frac16$ one obtains for $K=1$ the value of the sum cited in your question.