I edited the question as a result of the discussion in the comments. Originally my question was how to interpret the four discarded solutions. Now I'm making a step back and hope that someone can clarify in what sense it is sensible to discard four of the eight original solutions of the Dirac equation.
From making the ansatz ${\mathrm{e}}^{+ipx}$ and ${\mathrm{e}}^{-ipx}$, with $E=\pm \sqrt{ (\vec p)^2 +m^2} $ we get eight solutions of the Dirac equation. $u_1, u_2, u_3 , u_4$ and $v_1,v_2,v_3,v_4$.
Conventionally the four solutions ($u_3 , u_4,v_3,v_4$.) following from $E=- \sqrt{ (\vec p)^2 +m^2}$ are said to be linearly dependent of the remaining four solutions with $E=+\sqrt{ (\vec p)^2 +m^2}$ two ($u_1,u_2$) are commonly interpreted as particle and two ($v_1,v_2$) as antiparticle solutions.
Nevertheless, in order to be able to construct chirality eigenstates we need the other four solutions and I'm unsure in how far we can then say that four of the eight solutions are really linearly dependent. A chiral eigenstate must always be of the form $ \psi_L= \begin{pmatrix} f \\ -f \end{pmatrix} $ for some two component object $f$. In order to construct such an object we need all eight solutions. For example $\psi_L= u_1 - u_3$, as can be seen from the explicit form of the solutions recited below.
In addition, I'm unable to see that the eight solutions are really linearly dependent, because for me this means that we can find numbers $a,b,c,d,e,f,g,h \neq 0$, such that $a u_1 + b u_2 + c u_3 +d u_4 + e v_1 + f v_2 + g v_3 + h v_4=0$. As pointed out in the comments, this can be done, but only for one point in time. Is this really enough? In what sense is then for example the basis used in the Fourier expansion $\sum_n (a_n e^{in x} + b_n e^{-in x}) $ linearly independent? With the same reasoning we could find numbers for one $x$ to show that all these $e^{in x}$ and $e^{-in x}$ are linearly dependent...
The explicit solutions
This is derived for example here
Two solutions follow from the ansatz ${\mathrm{e}}^{-ipx}$ with $E=+ \sqrt{ (\vec p)^2 +m^2}$ and two with $E=- \sqrt{ (\vec p)^2 +m^2}$ .
In the rest frame the solutions are
$$ E=+ \sqrt{ (\vec p)^2 +m^2} \rightarrow u_1 = \begin{pmatrix} 1 \\ 0 \\0 \\ 0 \end{pmatrix} {\mathrm{e}}^{-imt} \qquad u_2 = \begin{pmatrix} 0 \\ 1 \\0 \\ 0 \end{pmatrix} {\mathrm{e}}^{-imt} $$
$$ E=- \sqrt{ (\vec p)^2 +m^2} \rightarrow u_3 = \begin{pmatrix} 0 \\ 0 \\1 \\ 0 \end{pmatrix} {\mathrm{e}}^{-imt} \qquad u_4 = \begin{pmatrix} 0 \\ 0 \\0 \\ 1 \end{pmatrix} {\mathrm{e}}^{-imt} $$
Analogous four solutions from the ansatz ${\mathrm{e}}^{+ipx}$, we get four solutions.
$$ E=+ \sqrt{ (\vec p)^2 +m^2} \rightarrow v_1 = \begin{pmatrix} 0 \\ 0 \\1 \\ 0 \end{pmatrix} {\mathrm{e}}^{imt} \qquad v_2 = \begin{pmatrix} 0 \\0 \\0 \\ 1 \end{pmatrix} {\mathrm{e}}^{imt} $$
$$ E=- \sqrt{ (\vec p)^2 +m^2} \rightarrow v_3 = \begin{pmatrix} 1 \\ 0 \\0 \\ 0 \end{pmatrix} {\mathrm{e}}^{imt} \qquad v_4 = \begin{pmatrix} 0 \\ 1 \\0 \\ 0 \end{pmatrix} {\mathrm{e}}^{imt} $$
Examples for chiral eigenstate are, with some two component object $f$
$$\psi_L = \begin{pmatrix} f \\ -f \end{pmatrix} \hat = u_1 - u_3 = \begin{pmatrix} 1 \\ 0 \\0 \\ 0 \end{pmatrix} {\mathrm{e}}^{-imt} - \begin{pmatrix} 0 \\ 0 \\1 \\ 0 \end{pmatrix} {\mathrm{e}}^{-imt} \qquad \text{ or } \qquad \psi_L = \begin{pmatrix} f \\ -f \end{pmatrix} \hat = u_2 -u_4 $$
$$\psi_L = \begin{pmatrix} f \\ -f \end{pmatrix} \hat = v_1 - v_3 \qquad \text{ or } \qquad \psi_L = \begin{pmatrix} f \\ -f \end{pmatrix} \hat = u_2 - u_4 $$
And similar for $\Psi_R = \begin{pmatrix} f \\ f \end{pmatrix}$.
Are four of the eight solutions really dependent? If yes, how can this be shown explicitly ? Any source, book, pdf would be awesome. Is it possible to interpret the solutions $(u_3,u_4,v_3,v_4)$ that can be discarded for many applications, but that are needed in order to create chirality eigenstates?