[guided]The goal is to turn $F$, which is an abstract bounded linear functional on $M_*$, into evaluation at an element of $M$. The vector functionals $\omega_{\eta,\zeta}$ provide the bridge, because they are normal functionals and therefore lie in the domain of $F$. We use the convention that Hilbert-space inner products are linear in the first variable and conjugate-linear in the second variable.
Define
\begin{align*}
b_F:H_u\times H_u\to\mathbb C
\end{align*}
by
\begin{align*}
b_F(\eta,\zeta)=F(\omega_{\eta,\zeta}).
\end{align*}
This is well-defined because $\omega_{\eta,\zeta}\in M_*$ for every $\eta,\zeta\in H_u$. The usual vector functional estimate says
\begin{align*}
\|\omega_{\eta,\zeta}\|_{M^*}\le \|\eta\|_{H_u}\|\zeta\|_{H_u}.
\end{align*}
Indeed, for every $y\in M$ with $\|y\|_{\mathrm{op}}\le 1$,
\begin{align*}
|\omega_{\eta,\zeta}(y)|=|(\pi_u(y)\eta,\zeta)_{H_u}|\le \|\pi_u(y)\eta\|_{H_u}\|\zeta\|_{H_u}\le \|\eta\|_{H_u}\|\zeta\|_{H_u}.
\end{align*}
Taking the supremum over such $y$ proves the estimate. Since $F$ is bounded,
\begin{align*}
|b_F(\eta,\zeta)|\le \|F\|_{(M_*)^*}\|\omega_{\eta,\zeta}\|_{M^*}\le \|F\|_{(M_*)^*}\|\eta\|_{H_u}\|\zeta\|_{H_u}.
\end{align*}
Thus $b_F$ is a bounded sesquilinear form. The [Riesz representation theorem](/theorems/221) for bounded sesquilinear forms gives a unique bounded operator $T_F\in\mathcal{L}(H_u)$ satisfying
\begin{align*}
(T_F\eta,\zeta)_{H_u}=b_F(\eta,\zeta)
\end{align*}
for all $\eta,\zeta\in H_u$.
We now prove that $T_F$ comes from the von Neumann algebra represented by $\pi_u$. To do this, it is enough to show that $T_F$ commutes with every operator in the commutant $\pi_u(M)'$. Let $R\in\pi_u(M)'$. Since $R$ commutes with every $\pi_u(x)$, for every $x\in M$ we compute
\begin{align*}
\omega_{R\eta,\zeta}(x)=(\pi_u(x)R\eta,\zeta)_{H_u}=(R\pi_u(x)\eta,\zeta)_{H_u}.
\end{align*}
Using the definition of the adjoint of $R$,
\begin{align*}
(R\pi_u(x)\eta,\zeta)_{H_u}=(\pi_u(x)\eta,R^*\zeta)_{H_u}.
\end{align*}
Therefore
\begin{align*}
\omega_{R\eta,\zeta}=\omega_{\eta,R^*\zeta}.
\end{align*}
Apply $F$ to this identity:
\begin{align*}
(T_FR\eta,\zeta)_{H_u}=F(\omega_{R\eta,\zeta})=F(\omega_{\eta,R^*\zeta}).
\end{align*}
By the definition of $T_F$ again,
\begin{align*}
F(\omega_{\eta,R^*\zeta})=(T_F\eta,R^*\zeta)_{H_u}=(RT_F\eta,\zeta)_{H_u}.
\end{align*}
Thus
\begin{align*}
(T_FR\eta,\zeta)_{H_u}=(RT_F\eta,\zeta)_{H_u}
\end{align*}
for all $\eta,\zeta\in H_u$, so $T_FR=RT_F$. Since $R\in\pi_u(M)'$ was arbitrary,
\begin{align*}
T_F\in \pi_u(M)''.
\end{align*}
Since $\pi_u$ is a faithful normal representation of the von Neumann algebra $M$, the range $\pi_u(M)$ is strongly closed in $\mathcal{L}(H_u)$. Thus $\pi_u(M)$ is a von Neumann algebra acting on $H_u$, and the bicommutant theorem identifies $\pi_u(M)''$ with $\pi_u(M)$. Hence $T_F=\pi_u(x_F)$ for a unique $x_F\in M$, where uniqueness follows from faithfulness of $\pi_u$.[/guided]