Using above formula, I could calculate the given commutator.
$$
[\epsilon^{\mu\nu\rho\sigma} M_{\mu \nu}M_{\rho\sigma},M_{\alpha\beta}]=i\epsilon^{\mu\nu\rho\sigma}(M_{\mu \nu}[M_{\rho\sigma},M_{\alpha\beta}]+[M_{\rho\sigma},M_{\alpha\beta}]M_{\mu \nu})
$$
(because...