I am not good at the properties of matrix commutators, thus I basically brute force it by picking entries from $A=\begin{pmatrix}a & b \\ c & d\end{pmatrix},B=\begin{pmatrix}e & f \\ g & h\end{pmatrix}$. Note how the diagonals give the strongest restriction to the elements as they have to add up to 1, i.e. after expanding $[A,B]=AB-BA$, we get $bg-fc=cf-hc=1$. Now setting $h=0$ to narrow down cases, we get $bg=2cf$.
Now because in the original equations they have to add up to 1, it means $b,g,c,f\neq 0$. Looking at the equations for the off diagonals, we have $af-eb-fd=-(d-a)f-eb=0=ce+dg-ga…