The set of all matrix-valued first-order differential operators that commute with the Dirac equation in n-dimensional complex Euclidean space is computed. In four dimensions it is shown that all matrix-valued second-order differential operators that commute with the Dirac operator in four dimensions are obtained as products of first-order operators that commute with the Dirac operator. Finally some additional coordinate systems for which the Dirac equation in Minkowski space can be solved by separation of variables are presented. These new systems are comparable to the separation in oblate spheroidal coordinates discussed by Chandrasekhar [S. Chandrasekhar, The Mathematical Theory of Black Holes (Oxford U.P., Oxford, 1983)].