In this contribution, we introduce higher spin Dirac operators, i.e. a specific class of differential operators in Clifford analysis of several vector variables, motivated by equations from theoretical physics. In particular, the higher spin Dirac operator in three vector variables will be explicitly constructed, starting from a description of the so-called twisted Rarita-Schwinger operator.