This one can also be done using complex variables.
Starting from the generating function
$$\sum_{n\ge 0} P_n(x) t^n
= \frac{1}{\sqrt{1-2xt+t^2}}$$
Call $P_n(0) = Q_n$ so that
$$\sum_{n\ge 0} Q_n t^n
= \frac{1}{\sqrt{1+t^2}}.$$
Now using Lagrange Inversion we get
$$Q_{2n} =
\frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{1}{z^{2n+1}}
\frac{1}{\sqrt{1+z^2}} \; dz
\\ = \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{1}{z^{2n+2}}
\frac{1}{\sqrt{1+z^2}} \; z \; dz.$$
Put $1+z^2 = w^2$ so that $z\; dz = w\; dw$ to obtain
$$2\times \frac{1}{2\pi i}
\int_{|w-1|=\epsilon} \frac{1}{(w^2-1)^{n+1}}
\frac{1}{w} \; w\; dw.$$
The scalar two in front is because when $z$ makes one turn around the
origin, $w$ makes two, but the contour only counts one turn.
Continuing we have
$$\frac{2}{2\pi i}
\int_{|w-1|=\epsilon} \frac{1}{(w-1)^{n+1}}
\frac{1}{(w+1)^{n+1}} \; dw
\\ = \frac{2}{2\pi i}
\int_{|w-1|=\epsilon} \frac{1}{(w-1)^{n+1}}
\frac{1}{(2+w-1)^{n+1}} \; dw
\\ = \frac{1}{2^{n}} \frac{1}{2\pi i}
\int_{|w-1|=\epsilon} \frac{1}{(w-1)^{n+1}}
\frac{1}{(1+(w-1)/2)^{n+1}} \; dw
\\ = \frac{1}{2^{n}} \frac{1}{2\pi i}
\int_{|w-1|=\epsilon} \frac{1}{(w-1)^{n+1}}
\sum_{q\ge 0} (-1)^q {q+n\choose n}
\frac{(w-1)^q}{2^q} \; dw
.$$
Extracting coefficients this gives
$$\frac{1}{2^{n}} (-1)^n {2n\choose n} \frac{1}{2^n}.$$
This is
$$\frac{(-1)^n}{4^n} \frac{(2n)!}{(n!)^2}$$
as claimed.