The 1st equation:
$$|\psi\rangle=c_+|+\rangle+c_-|-\rangle$$
is the generic form for any pure state in a two-level system. The space of all normed pairs of ($c_+$, $c_-$) are called the Bloch Sphere.
Being that it applies to any two level system, it does not elucidate the physical nature of a spinor particularly well. In fact, in many teaching scenarios, it may lead to confusion. This is because of the terms "spin-up" and "spin-down". This leads the uninitiated into thinking the spin is either aligned or anti-aligned with magnitude $\frac 1 2$ (in units of $\hbar$), this is not the case.
A spin 1/2 particle has an angular momentum of
$$\sqrt{s(s+1)} = \frac{\sqrt 3} 2$$
and the eigenstates allows to know only the projection on a single axis, which is $\frac 1 2$. That means, that in an $S_z$ eigenstate, there is still:
$$ S_{x \& y} = \sqrt{s(s+1)-s^2} = \frac 1 {\sqrt 2}$$
units of angular momentum that is equally distributed in the $x-y$ plane. That fact, along with the commutations relations:
$$ [S_i, S_j] = 2i\epsilon_{ijk}S_k$$
and the ability to quantize the spin on any axis are critical to understanding spinors.
In short, you can write the coefficients as:
$$ c_+ = \cos{\theta/2}e^{-i\phi/2} $$
$$ c_- = \sin{\theta/2}e^{+i\phi/2} $$
and, while any pair of ($\theta,\phi$) represents a superposition in the original basis, it is a pure spin up state when quantized along an axis pointing in the ($\theta, \phi$) direction.
So that's the spin part of a wave function. If you write the spin wave function as a column vector (your 2nd equations) and combine it with a factored spatial wave function, you get:
$$\Psi(\vec{x})=\psi( \vec x)\begin{pmatrix}c_+\\ c_-\end{pmatrix}$$
That represents a particle, with a spatial wave function, and an unentangled spin. In the most general case, the space and spin components can be entangled, and that requires your 3rd equations:
$$\psi(\vec{x})=\begin{pmatrix}\psi_{1/2}( \vec x) \\ \psi _{-1/2}(\vec{x})\end{pmatrix}$$
This is required for systems like the Stern-Gerlach experiment.
If we look at just one spatial component of the wave function, e.g., $\psi_+(\vec x)$, it's written here as a function of a vector, $\vec x$. That can be done in cartesian coordinates:
$$ \psi(\vec x) = \psi(x,y,z)$$
other coordinates, or, what is generally the most useful when angular momentum is conserved: spherical coordinates. In that case, the angular part is captured in spherical harmonics, per that standard projection in your 5th equation:
$$\sum_{l,m}\sum_{s_z=-s}^s\langle \vec{x} | l \ m \ s_z\rangle \langle l \ m \ s_z | \psi\rangle=\sum _{l,m,s_z}Y_{l,m}(\theta , \phi)u_{s_z}c^{s_z}_{l,m}(r) \ \ \ \ \ \ (1)$$
what we see from this equation is that the orbital angular momentum ($(l,m)$) couple to each spin state (per the usual rules of angular momentum addition). Note that eigenstates of total angular momentum quantum numbers ($J^2, J_z$) mix different $(l, m)$ with spin up and spin down, as captured by the $c^{\pm}_{lm}$. This is standard Clebsch-Gordan stuff, e.g.:
$$|\frac 3 2 \frac 1 2\rangle = \begin{pmatrix} \sqrt{\frac 2 3}Y_{1,0}(\theta, \phi) \\ \sqrt{\frac 1 3}Y_{1,1}(\theta, \phi)\end{pmatrix}$$
So that represents and entanglement of the spin and angular coordinates such that total angular momentum quantum numbers are good.
When one sees this in, say, the Hydrogen atom, the radial wave function is entirely factored out. Your final equation, however, allows radial dependence in the $c^{\pm}_{lm}(r)$. This introduces spatial entanglement between spin and the space coordinate in the radial direction. (I cannot think of a system off hand where that comes up...since, if you're using spherical harmonics, the problem should be separable in spherical coordinates). Nevertheless, it is the most general spinor wave function in spherical coordinates.
Addendum:
Note that the $\theta$ and $\phi$ dependence on the RHS of (1) is directly related to $\vec x$ on the LHS. Here I am going to $\vec r$ for a 3D position vector, as $\vec x$ will cause confusion:
In cartesian coordinates $(x, y, z)$ and the cartesian basis:
$${\bf r} = x{\bf\hat x}+y{\bf\hat y}+z{\bf\hat z} $$
which is, in spherical coordinates $(r, \theta, \phi)$ and the cartesian basis:
$${\bf r} = r\sin\theta\cos\phi{\bf\hat x}+r\sin\theta\sin\phi{\bf\hat y}+r\cos\theta{\bf\hat z} $$
Now we can switch from the cartesian basis vectors to spherical basis vectors (spherical basis vectors are the 3 eigenvectors of rotations about the z-axis, they are not the spherical-coordinate basis vectors ($\hat r, \hat\theta, \hat\phi)$, and understanding them is very helpful in understanding spin, and tensor operators, etc):
$${\bf r} = r\sin\theta\cos\phi\frac 1 {\sqrt 2}({\bf\hat e^+}+{\bf\hat e^-})+r\sin\theta\sin\phi \frac i {\sqrt 2}({\bf\hat e^+}-{\bf\hat e^-})+r\cos\theta{\bf\hat e^0} $$
which can be rearranged:
$${\bf r} = \frac 1 {\sqrt 2}r\sin\theta e^{+i\phi}{\bf\hat e^+}+\frac i {\sqrt 2}r\sin\theta e^{-i\phi} {\bf\hat e^-}+r\cos\theta{\bf\hat e^0} $$
note that the angular coefficients are the $l=1$ spherical harmonics:
$${\bf r} = r\sqrt{\frac{4\pi} 3}\big[Y_1^1(\theta,\phi){\bf\hat e^+} + Y_1^0(\theta,\phi){\bf\hat e^0} +Y_1^{-1}(\theta,\phi){\bf\hat e^-}
\big ]$$
That formula may take some getting used to. In the spherical basis, the spherical harmonics are the coordinates. It's strange, how can a function be a coordinate? Well, look at $Y_1^0(\theta,\phi)$. On the unit sphere, it is $z$. Well $z$ is a function (it's a polynomial), and so are $x$ and $y$. So we use polynomial coefficients in the straight cartesian rep, we just don't usually think of them as such, we think of them as labeled numbers.