Clebsh-Gordan coefficients don't play a part in the definition of LQG. The notion of intertwining tensors is enough to build spin network states.
They are, however, a very convenient tool for practical calculations.
This has to do with representation theory of $\mathfrak{su}_2$. A spin-$j$ irrep can be thought of as the totally symmetric part of a tensor product of $2j$ fundamental (spin-$1/2$) irreps.
You can also implement intertwiners between $V_{\{1,2,3\}}$ as tensors acting on the tensor product $V_{1/2}^{\otimes 2 j_1} \otimes V_{1/2}^{\otimes 2 j_2} \otimes V_{1/2}^{\otimes 2 j_3}$, and later taking only the projection of these tensors on the symmetric subspace. This is convenient, because there exists a very simple graphical calculation model for these – see Penrose's binor calculus.
Clebsh-Gordan coefficients are essentially the ratios between the normalization of intertwiners in Penrose's binor calculus, and the renormalization of intertwiners in LQG (given by the spin network inner product).
They can be used, for example, to easily calculate the values of the $6j$ symbol that enters the definition of the Ponzano-Regge model (quantum gravity in $3d$ Euclidean space-time). You replace normalized intertwiners with symmetrized binor calculus diagrams which makes their evaluation trivial, and multiply by the Clebsh-Gordan coefficients to get the correct answer.