Can we claim that for an orthogonal matrix $A$ (satysfying $AA^T=I$), its symmetric and skew-symetric parts are always commuting?
Symmetric part is calculated as $S=\frac{1}{2}(A+A^T)$ and skew symmetric part as $K=\frac{1}{2}(A-A^T)$.
It's easy to show the claim for dimension $2$.
Also for dimension $3$ we can see commutation for orthogonal matrices with determinant equal to 1.
For a rotation matrix $R$ it's evident from the Rodrigues formula,
when $R=\color{blue}{I+(1-\cos\theta) V^2 }+ \color{green}{\sin\theta V }$ for some skew-symmetric matrix $V $.
Here commutation follows because symmetric part and skew-symmetric part are expressed as polynomials of the same matrix $V $.
Is it possible to find similar explanation for higher dimensions of orthogonal matrices?
One can verify it using a direct computation:
$$ SK = \frac{1}{2} \left( A + A^T \right) \cdot \frac{1}{2} \left(A - A^T \right) = \frac{1}{4} \left( A^2 - AA^T + A^TA - \left( A^T \right)^2 \right) = \frac{1}{4} \left( A^2 - \left( A^T \right)^2 \right)$$ while $$ KS = \frac{1}{2} \left( A - A^T \right) \cdot \frac{1}{2} \left( A + A^T \right) = \frac{1}{4} \left( A^2 + AA^T - A^TA - \left( A^T \right)^2 \right) = \frac{1}{4} \left( A^2 - \left( A^T \right)^2 \right)$$
where we used the fact that $AA^T = A^TA$ (as they both are equal to $I_n$). More generally, the symmetric and and skew-symmetric part of a matrix commute iff the matrix is normal and an orthogonal matrix is normal.