The first-order fermionic kinetic term produces second-class momentum constraints. After eliminating them, the canonical graded Dirac brackets are
with other mixed brackets zero. The bracket of two odd variables is symmetric, consistently with the graded Poisson bracket. The constraint algebra is
It is the worldline supersymmetry algebra: the square of the fermionic constraint generates the Hamiltonian constraint.
With , quantization converts these to
Represent the fermions by Gamma matrices, , where , and . The wavefunction is a Dirac spinor. Its physical-state constraint is
which is exactly the massless Dirac equation. Moreover , so this equation implies the remaining mass-shell constraint. Conversely, a solution of the massless Dirac equation satisfies both constraints. The Klein-Gordon constraint alone would not impose the spinorial first-order equation.