Equivalently, in the Lie Algebra $M(2,q)$, how can I show that there are precisely $q^2$ solutions M to $[M,X]=0,$ where $X$ is a non central element of $\mathop{GL}(2,q)$, where $q$ is an odd prime?
It is easy to show that there are at least $q^2$ solutions and less than $q^4$, as the set of solutions forms a Lie subalgebra of $M(2,q)$ with the identity matrix and $X$ itself in this set of solutions, so the dimension of this Lie subalgebra is at least $2$. Assuming the solution set is all of $M(2,q)$ leads to a contradiction that $X$ is non-central, as substituting the matrix of 1s for $M$ and then the matrix with 1s in the top row and 0s in the bottom row leads to $X$ being central in $M(2,q)$. So all that remains is to show that dimension of the subalgebra of solutions cannot be $3$.
Write $M(2,q) = \langle I, J, D, E \rangle_{\mathbb F_q}$ where $$ I = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix}, \quad J = \begin{bmatrix} 1 & 0 \\ 0 & -1 \end{bmatrix}, \quad D = \begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix}, \quad E = \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix}. $$ Since this problem is invariant under a change of basis (I mean the basis of $\mathbb F_q^2$, not of $M(2,q)$), we can assume without loss of generality that $X = J$ (note that $J$ is not central and that $[M,X] = 0$ if and only if $[PMP^{-1},PXP^{-1}] = 0$ for $P \in \mathop{GL}(2,q)$, so the "without loss of generality" makes sense). You can explicitly compute the Lie brackets $[J,I]=0$, $[J,J] = 0$, $[J,D] = \begin{bmatrix} 0 & 2 \\ -2 & 0 \end{bmatrix}$ and $[J,E] = \begin{bmatrix} 0 & 2 \\ 2 & 0 \end{bmatrix}$ to notice that the subspace of matrices $Y$ satisfying $[J,Y]=0$ has dimension $2$.
Not exactly a smart argument, but it works.
Hope that helps,