[guided]Let $X \in \mathfrak{X}(P)$ be horizontal and let $\xi \in \mathfrak{g}$. We want to evaluate the proposed identity on $(X,\xi_P)$, where $\xi_P$ is the fundamental vertical vector field defined by
\begin{align*}
\xi_P(q)=\frac{d}{dt}\bigg|_{t=0}q\exp(t\xi).
\end{align*}
The curvature form is horizontal by definition: it only sees the horizontal components of its inputs. Since the vertical vector field $\xi_P$ has horizontal component $0$,
\begin{align*}
\Omega(X,\xi_P)=0.
\end{align*}
Now compute the bracket term. Because $X$ is horizontal, $\omega(X)=0$. Because $\omega$ is a principal connection form, $\omega(\xi_P)=\xi$. Therefore
\begin{align*}
\frac{1}{2}[\omega,\omega](X,\xi_P) = [\omega(X),\omega(\xi_P)]_{\mathfrak{g}} = [0,\xi]_{\mathfrak{g}} = 0.
\end{align*}
So the only remaining point is the exterior derivative term. The formula for the exterior derivative of a $\mathfrak{g}$-valued one-form gives
\begin{align*}
d\omega(X,\xi_P) = X(\omega(\xi_P))-\xi_P(\omega(X))-\omega([X,\xi_P]).
\end{align*}
The first term is zero because $\omega(\xi_P)$ is the constant $\mathfrak{g}$-valued function with value $\xi$. The second term is zero because $\omega(X)=0$.
Thus we must prove $\omega([X,\xi_P])=0$. This is exactly where principal equivariance enters, but we need its infinitesimal form rather than only the statement that right translations preserve horizontal subspaces. For $g \in G$, define $R_g: P \to P$ by $R_g(q)=qg$, and let $(dR_g)_q: T_qP \to T_{qg}P$ be its differential. The flow of $\xi_P$ is given by
\begin{align*}
\Phi_t^{\xi_P}=R_{\exp(t\xi)},
\end{align*}
where $\exp: \mathfrak{g} \to G$ is the Lie group exponential map. The connection form satisfies the principal equivariance identity
\begin{align*}
R_g^*\omega=\operatorname{Ad}_{g^{-1}}\omega.
\end{align*}
Here $\operatorname{Ad}: G \to GL(\mathfrak{g})$ is the adjoint representation. Therefore
\begin{align*}
(\Phi_t^{\xi_P})^*\omega = \operatorname{Ad}_{\exp(-t\xi)}\omega.
\end{align*}
Differentiating at $t=0$ gives
\begin{align*}
\mathcal{L}_{\xi_P}\omega = -\operatorname{ad}_{\xi}\circ \omega,
\end{align*}
where $\operatorname{ad}_{\xi}: \mathfrak{g} \to \mathfrak{g}$ is the linear map $\eta \mapsto [\xi,\eta]_{\mathfrak{g}}$.
Now evaluate this one-form identity on the horizontal vector field $X$. The Lie derivative of a one-form satisfies
\begin{align*}
(\mathcal{L}_{\xi_P}\omega)(X)=\xi_P(\omega(X))-\omega([\xi_P,X]).
\end{align*}
Since $\omega(X)=0$, the left side also equals $-\omega([\xi_P,X])$. The infinitesimal equivariance identity gives
\begin{align*}
-\omega([\xi_P,X]) = -[\xi,\omega(X)]_{\mathfrak{g}} = 0.
\end{align*}
Hence $\omega([\xi_P,X])=0$. Since $[X,\xi_P]=-[\xi_P,X]$, we also have
\begin{align*}
\omega([X,\xi_P])=0.
\end{align*}
Substituting into the exterior derivative formula gives
\begin{align*}
d\omega(X,\xi_P)=0.
\end{align*}
Combining the three computations,
\begin{align*}
\left(d\omega+\frac{1}{2}[\omega,\omega]\right)(X,\xi_P)=0=\Omega(X,\xi_P).
\end{align*}[/guided]