Of course, it must read
\left (\frac{\partial \hat{\rho}}{\partial t} \right )_{\text{expl}}+\frac{1}{\mathrm{i} \hbar} [\hat{\rho},\hat{H}]=0.
[I corrected the previous posting for a sign mistake and the missing =0].
This von Neuamnn equation is valid in any picture of time evolution.
To make ist more clear one must remember that there is a set of operators generating the algebra of observables that are by definition not explicitly time dependent. In QFT these are the fundamental field operators, in non-relativistic QM (in the first-quantized formulation) it's the position, moomentum, and spin operators, etc.
The mathematical time dependence of observables and states is only determined up to a time-dependent unitary transformation and this leads to the different choices of the picture of time evolution. The most general choice is the Dirac picture. It turns out that the time evolution is always local in time, and you have two equations of motion for the unitary time-evolution operators for the observable operators and the statistical operator (using \hbar=1 for simplification),
\begin{split}<br />
\partial_t \hat{A}(t,t_0)&=\mathrm{i} \hat{H}_1(t) \hat{A}(t,t_0), \\<br />
\partial_t \hat{C}(t,t_0)&=\mathrm{i} \hat{H}_2(t) \hat{C}(t,t_0). \\<br />
\end{split}<br />
The Hermitean operators \hat{H}_1 and \hat{H}_2 are arbitrary but constrained by the demand that
\hat{H}=\hat{H}_1+\hat{H}_2.
The time evolution of the observable statistical operators are then given by
<br />
\begin{split}<br />
\hat{O}(t) &=\hat{A}(t,t_0) \hat{O}(t_0) \hat{A}^{\dagger}(t,t_0),\\<br />
\hat{\rho}(t) &=\hat{C}(t,t_0) \hat{\rho}(t_0) \hat{C}^{\dagger}(t,t_0).<br />
\end{split}
From this you get the equations of motion in the so chosen picture as
\begin{split}<br />
\frac{\mathrm{d}}{\mathrm{d} t} \hat{O}&= -\mathrm{i}[\hat{O},\hat{H}_1] + \left (\frac{\partial \hat{O}}{\partial t} \right )_{\text{expl}},\\<br />
\frac{\mathrm{d}}{\mathrm{d} t} \hat{\rho}&= \mathrm{i}[\hat{\rho},\hat{H}_2].<br />
\end{split}<br />
The explicit time dependence refers to the time dependence which does not come by implicit time dependence of the fundamental operators that build up \hat{O}.
Of course, as a function of observables and perhaps explicitly of time the Statistical operator also must fulfill the equation of motion
\frac{\mathrm{d}}{\mathrm{d} t} \hat{\rho}= -\mathrm{i}[\hat{O},\hat{H}_1] + \left (\frac{\partial \hat{\rho}}{\partial t} \right )_{\text{expl}}.
Subtracting both expressions for the time derivative of the Statistical Operator one obtains von Neumann's equation:
-\mathrm{i}[\hat{\rho},\hat{H}]+\left (\frac{\partial \hat{\rho}}{\partial t} \right )_{\text{expl}}=0,
because \hat{H}_1+\hat{H}_2=\hat{H}.
