Consider a system with Hamiltonian $\hat{H}=\hat{H_0}+\hat{V}$. We define the interaction picture kets $|\psi(t)\rangle _I$ by $$\tag{1} |\psi(t)\rangle _I=\exp\left(\frac{i}{\hbar}\hat{H}_0(t-t_0)\right)|\psi(t)\rangle_{S}$$ and interaction picture operators $\hat{O}_I$ by $$\tag{2} \hat{O}_I= \exp\left(\frac{i}{\hbar}\hat{H}_0(t-t_0)\right)\hat{O}\exp\left(-\frac{i}{\hbar}\hat{H}_0(t-t_0)\right).$$
From the ordinary Schrödinger equation, we obtain the equation of motion for the state: $$\tag{3} i\hbar\frac{d}{dt}|\psi(t)\rangle _I = \hat{V}_I|\psi(t)\rangle _I.$$ As is done e.g. in Sakurai's book, we also define the time-evolution operator $\hat{U}(t,t_0)_I$ by $$\tag{4} |\psi(t)\rangle _I = \hat{U}(t,t_0)_I |\psi(t_0)\rangle _I.$$ Differentiating (4) w.r.t $t$ yields $$\tag{5} \begin{align} i\hbar\frac{d}{dt}|\psi(t)\rangle _I &= i\hbar\frac{d}{dt}\left(\hat{U}(t,t_0)_I\right) |\psi(t_0)\rangle_{I} \\ &=i\hbar\frac{d}{dt}\left(\hat{U}(t,t_0)_I\right)\hat{U}(t,t_0)_I^{-1} |\psi(t)\rangle_{I}. \end{align}$$ Comparison of (3) and (5) shows that we must have the following equation of motion for $\hat{U}(t,t_0)_I$: $$\tag{6} i\hbar\frac{d}{dt}\hat{U}(t,t_0)_I = \hat{V}_I\hat{U}(t,t_0)_I.$$ However, if we take $\hat{O}=\hat{U}(t,t_0)$ in equation (2), then equation (2) would yield $$\tag{7} \begin{align} &i\hbar\frac{d}{dt}\hat{U}(t,t_0)_I\\ &= \exp\left(\frac{i}{\hbar}\hat{H}_0(t-t_0)\right)(-\hat{H}_0\hat{U}(t,t_0)+\hat{U}(t,t_0)\hat{H}_0)\exp\left(-\frac{i}{\hbar}\hat{H}_0(t-t_0)\right)\phantom{abc}\\&+\exp\left(\frac{i}{\hbar}\hat{H}_0(t-t_0)\right)\hat{H}\hat{U}(t,t_0)\exp\left(-\frac{i}{\hbar}\hat{H}_0(t-t_0)\right)\\ &= \hat{V}_I\hat{U}(t,t_0)_I+\hat{U}(t,t_0)_I\hat{H}_0. \end{align}$$
Evidently, equations (6) and (7) do not agree. Am I misunderstanding something here?