Let $H = {\left( {{L^2}\left( {{\mathbb{R}^3},{\mathbb{C}^4}} \right)} \right)^4}$ be the Hilbert space in question and for $u = \left( {{u_1},{u_2},{u_3},{u_4}} \right),v = \left( {{v_1},{v_2},{v_3},{v_4}} \right) \in H$ let the inner product on $H$ be defined by $\left\langle {u,v} \right\rangle = \sum\limits_{i = 1}^4 {\int_{ - \infty }^{ + \infty } {{u_i}\left( {{x_i}} \right)\overline {{v_i}\left( {{x_i}} \right)} d{x_i}} } $.
${\left[ {\nabla u} \right]_{i,j}} = {\partial _j}{u_i}$ so
${\left[ {Du} \right]_{i,j}} = - i\sum\limits_{k = 1}^4 {{\alpha _{i,k}}{{\left[ {\nabla u} \right]}_{k,j}}} + {\beta _{i,j}} = - i\sum\limits_{k = 1}^4 {{\alpha _{i,k}}{\partial _j}{u_k}} + {\beta _{i,j}}$
$D:H \to {\left( {{M_{4 \times 4}}\left( \mathbb{C} \right)} \right)^H}$, where ${\left( {{M_{4 \times 4}}\left( \mathbb{C} \right)} \right)^H}$ is a vector space of functions from $H$ to $4 \times 4$ matrices with complex-valued entries.
This means that $\left\langle {Du,v} \right\rangle $ is ill-defined, since $\left\langle { \cdot \,,\, \cdot } \right\rangle :H \times H \to \mathbb{C}$.
Suppose that we somehow redefine $D$ or the inner product so that $\left\langle {Du,v} \right\rangle $ becomes well-defined.
In that case, it would suffice to show that
1) $\left\langle {Du,v} \right\rangle - \left\langle {u,Dv} \right\rangle = 0$ (integration by parts combined with boundary conditions $u\left( { \pm \infty } \right) = v\left( { \pm \infty } \right) = 0$ - which must hold since $u,v \in H$ - should provide the needed equality). We conclude that $D$ is symmetric.
2) $D$ is well-defined on its domain. There should be no problems verifying this
It's easy to show that a symmetric operator is linear, so Hellinger-Toeplitz theorem then implies that $D$ is bounded and hence hermitian.
Note: In the text above, the definitions in the accepted answer here are used, meaning that hermitian implies self-adjoined implies symmetric and converse implications don't necessarily hold.