I think, one has to be a bit more careful with which quantity you are considering. In addition the time dependence of the state vectors and the operators representing observables is dependent on the picture of the time dependence, but the physical outcome of quantum mechanics is independent of the choice of the picture. For simplicity's sake, I'll work in the Schrödinger picture, where the state kets (or equivalently and more general the statistical operator) evolve with the time-evolution operator generated from the full Hamiltonian of the system, which I assume to be not explicitly time dependent. The equation of motion for that operator reads
\frac{\partial}{\partial t} \hat{C}(t,t_0)=-\mathrm{i} \hat{H} \hat{C}(t,t_0), \quad \hat{C}(t,t_0)=\hat{1}.
I've set \hbar=1. The solution of this equation reads
\hat{C}(t,t_0)=\exp[-\mathrm{i} H (t-t_0)].
This is obviously a unitary operator.
The operators representing observables are constant in time (if not explicitly time dependent, which case I exclude here).
Now let |o_1,\ldots,o_n \rangle denote the simultaneous (generalized) eigenvector of a complete set of operators. Further let |\psi_0 \rangle denote the state ket at time, t_0. The state ket at any later time is then given by
|\psi,t \rangle=\hat{C}(t,t_0) |\psi_0 \rangle.
The wave function wrt. to the complete set of (generalized) eigenvectors is given by
\psi(o_1,\ldots,o_n,t)=\langle o_1,\ldots,o_n|\psi,t \rangle
and generally time dependent. It's modulus squared gives the probability (densitity) for simultaneously getting the values o_j when measuring the complete set of compatible observables. This probability is of course also time dependent.
The wave function obeys the equation of motion,
\mathrm{i} \partial_t \psi(o_1,\ldots,o_n,t) = \hat{H}(o_1,\ldots o_n) \psi(o_1,\ldots,o_n,t),
where \hat{H}(o_1,\ldots,o_n) denotes the Hamiltonian in the representation of the complete set of operators. In general that's a differential operator (as, e.g., in the position representation) or an integral kernel (as, e.g., in the momentum representation). This is of course nothing else than the Schrödinger equation for the wave function and, btw., independent on the picture of time evolution as it should be since the wave function (or better said its modulus squared) is a physically observable probability.
Now, if you have a complete set of energy eigenvectors, |E,\alpha_j, where the \alpha_j count a possible degeneracy of the energy eigenvectors, then the time evolution is particularly simple since
\psi(E,\alpha_j,t)=\langle E,\alpha_j|\exp(-\mathrm{i} \hat{H} t)|\psi_0 \rangle<br />
=\langle \exp(+\mathrm{i} \hat{H} t E,\alpha_j|\psi_0 \rangle = \langle \exp(+\mathrm{i} E t E,\alpha_j|\psi_0 \rangle = \exp(-\mathrm{i} E t) \psi(E,\alpha_j,t_0).
The probability to measure an energy, E (eventually simultaneously with some other observables if there exist observables which are compatible with the Hamiltonian, i.e., that are conserved quantities that are not pure functions of the Hamiltonian itself) is therefore time independent.