Consider the construction of the Borel Functional Calculus for a self adjoint operator $A$ as descibred here:
Continuity of the functional calculus form of the Spectral Theorem
or better yet, here:
http://www.math.mcgill.ca/jakobson/courses/ma667/mendelsontomberg-spectral.pdf on page 4.
Assume the continuous functional calculus was established.
I want to show: $AB = BA \implies B\Phi(f) = \Phi(f)B$ where $f \in \mathbb{B}(\mathbb{R})$ the bounded Borel functions on $\mathbb{R}$.
See my suggestion for a simpler proof below please.
The selfadjoint operators enjoy monotone convergence properties in the strong operator topology. Specifically, suppose $\{ A_n \}_{n=1}^{\infty}$ is a sequence of bounded selfadjoint operators on a Hilbert space $\mathcal{H}$ such that $$ A_1 \le A_2 \le A_3 \le \cdots \le A_{\infty} $$ for some selfadjoint operator $A_{\infty}$. Then $\lim_n A_n x=Ax$ exists for all $x$ and some selfadjoint operator $A$. Using this, you can extend the functional calculus for continuous functions $f$ to include characteristic functions $f=\chi_{E}$ of open, closed, and half-open intervals through strong limits. Hence, you can extend the commutativity $B\Phi(f)=\Phi(f)B$ to $f=\chi_{E}$, where $E$ is in the sigma algebra generated by such sets, which includes the Borel sigma algebra; indeed, you may extended to countable unions of sets because of the above monotone property, while intersections can be handled starting with $\Phi(f)\Phi(g)=\Phi(fg)$ for continuous functions. Thus $B\Phi(\chi_{E})=\Phi(\chi_{E})B$ must hold for all Borel sets $E$. From there, using monotone convergence properties, it follows that $B\Phi(f)=\Phi(f)B$ for all bounded Borel measurable functions $f$.