Obtenemos para el LHS que es
$$[z^k] \frac{1}{1-z} \left(1+2z+2z^2+\cdots\right)^n = [z^k] \frac{1}{1-z} \left(1+2\frac{z}{1-z}\right)^n \\ = [z^k] \frac{(1+z)^n}{(1-z)^{n+1}}.$$
Extrayendo los coeficientes obtenemos
$$\sum_{q=0}^k {n\choose q} {k-q+n\choose n} = \sum_{q=0}^k {n\choose q} {k-q+n\choose k-q}.$$
Introducir la integral $${k-q+n\choose k-q} = \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{k-q+n}}{z^{k-q+1}} \; dz.$$
Esto refuerza el alcance (el polo se desvanece cuando $q\gt k$ ) por lo que podemos extender la suma a $n$ , obteniendo
$$\frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{k+n}}{z^{k+1}} \sum_{q=0}^n {n\choose q} \frac{z^q}{(1+z)^q} \; dz \\ = \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{k+n}}{z^{k+1}} \left(1+\frac{z}{1+z}\right)^n \; dz \\ = \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{k}}{z^{k+1}} (1+2z)^n \; dz.$$
Lo evaluamos utilizando el negativo del residuo en el infinito para encontrar (los residuos suman cero)
$$\mathrm{Res}_{z=\infty} \frac{(1+z)^{k}}{z^{k+1}} (1+2z)^n = -\mathrm{Res}_{z=0} \frac{1}{z^2} (1+1/z)^k z^{k+1} (1+2/z)^n \\ = -\mathrm{Res}_{z=0} \frac{1}{z} (1+z)^k \frac{(2+z)^n}{z^n} = -\mathrm{Res}_{z=0} \frac{1}{z^{n+1}} (1+z)^k (2+z)^n.$$
Extrayendo el residuo y volteando el signo obtenemos
$$\sum_{q=0}^n {n\choose q} 2^{n-q} {k\choose n-q}.$$
Este es el RHS como se reclama.
Observación. Esto también funciona cuando $k \gt n$ porque cuando bajamos el límite superior de la suma a $n$ estamos omitiendo términos que de todos modos habrían sido nulos debido a la ${n\choose q}$ plazo.