Starting with your third to last line, we begin by rewriting \begin{equation} \begin{split} (\psi^\dagger)^*(\gamma^0)^* (\gamma^\mu)^* \psi^* &= \psi^T \big[(\gamma^0)^\dagger\big]^T \big[(\gamma^\mu)^\dagger\big]^T (\psi^\dagger)^T \\ &= \big[\psi^\dagger (\gamma^\mu)^\dagger (\gamma^0)^\dagger \psi \big]^T\\ &= \psi^\dagger (\gamma^\mu)^\dagger (\gamma^0)^\dagger \psi \end{split} \end{equation}
where in going from penultimate to last line we have used that the components of the current are complex numbers and thus not matrix valued, such that we may drop the transpose. We may then proceed in a way similar to my answer to this question, using the following properties of the gamma matrices
\begin{align} (\gamma^0)^\dagger &= \gamma^0, \\ (\gamma^\mu)^\dagger &= \gamma^0 \gamma^\mu \gamma^0, \\ (\gamma^0)^2 &= \mathbb{I}_{4}, \end{align}
where $\mathbb{I}_{4}$ is the identity to write
\begin{equation} \begin{split} \psi^\dagger(\gamma^\mu)^\dagger(\gamma^0)^\dagger \psi &= \bar{\psi}\gamma^\mu(\gamma^0)^2\psi\\ &= \bar{\psi} \gamma^\mu \psi. \end{split} \end{equation}
This is then the result $j^\mu_c = j^\mu$. This is a consequence of the charge conjugation symmetry of quantum electrodynamics.