It seems the condition on $d$ needs to be revised somewhat. In the following we use the coefficient of operator $[q^k]$ to denote the coefficient of $q^k$ of a series. Given the Laurent polynomial
\begin{align*}
\sum_{s=N}^{M} \gamma_s q^s
\end{align*}
we have $M\geq N$. We write the condition $d\geq M-N$ as
\begin{align*}
\color{blue}{d=M-N+e}\tag{1}
\end{align*}
with $e$ a non-negative integer.
We obtain
\begin{align*}
\color{blue}{[q^d]}&\color{blue}{\sum_{s=N}^M\gamma_sq^s\sum_{i=0}^{\infty}q^i}\\
&=\sum_{s=N}^M\gamma_s[q^{d-s}]\sum_{i=0}^{\infty}q^i\tag{2.1}\\
&=\sum_{s=0}^{M-N}\gamma_{s+N}[q^{d-s-N}]\sum_{i=0}^{\infty}q^i\tag{2.2}\\
&=\sum_{s=0}^{M-N}\gamma_{M-s}[q^{d-M+s}]\sum_{i=0}^{\infty}q^i\tag{2.3}\\
&\,\,\color{blue}{=\sum_{s=0}^{M-N}\gamma_{M-s}[q^{e-N+s}]\sum_{i=0}^{\infty}q^i}\tag{2.4}\\
&\stackrel{?}=\sum_{s=0}^{M-N}\gamma_{M-s}\tag{2.5}\\
\end{align*}
Comment:
In (2.1) we use the linearity of the coefficient of operator and apply the rule $[q^{k-l}]A(q)=[q^k]q^lA(q)$.
In (2.2) we shift the index to start with $s=0$.
In (2.3) we change the order of summation $s\to M-N-s$.
In (2.4) we apply (1).
Let's have a closer look at $[q^{e-N+s}]$ in (2.4). Since for each $s$ we want to select the summand $\gamma_{M-s}$ in order to obtain (2.5), we need a non-negative exponent of $q$. We therefore have the condition
\begin{align*}
\color{blue}{e-N+s\geq 0\qquad\qquad 0\leq s\leq M-N}\tag{3}
\end{align*}
$N\leq 0$: Since $e\geq 0$, (3) is fulfilled whenever $N\leq 0$.
$N>0$: The situation is different when $N$ is positive. Setting >$s=0$ we get from (1) and (3)
\begin{align*}
e&\geq N\\
d=M-N+e&\geq (M-N)+N\\
\color{blue}{d}&\color{blue}{\geq M} \tag{4}
\end{align*}
Conclusion: In case $N>0$ we need $d\geq M$ instead of $d\geq M-N$. The condition on $d$ can therefore be stated as
\begin{align*}
\color{blue}{d\geq \max\{M,M-N\}}
\end{align*}
Plausibility check $N=2, M=3, d=1$:
Here we have $N>0$ and violate (4) since $d=M-N=3-2=1$ is less than $M=3$. We obtain
\begin{align*}
[q^1]&\left(\gamma_2q^2+\gamma_3q^3\right)\sum_{i=0}^{\infty}q^i\\
&=\left(\gamma_2[q^{-1}]+\gamma_3[q^{-2}]\right)\sum_{i=0}^{\infty}q^i\\
&=0
\end{align*}
We see although we have $1=d\geq M-N=1$, we get $0$ instead of $\gamma_2+\gamma_3$.