Let's work in the Schrödinger picture. Then the state vectors carry the entire time evolution, while the essentially self-adjoint (Hermitean is not sufficient!) operators that represent the observables are time-independent (we leave out the somwhat more complicated case of explicitly time-dependent observables). Then the state vector obeys
$$\mathrm{i} \hbar \partial_t |\psi(t) \rangle = \hat{H} |\psi(t) \rangle.$$
Obviously the formal solution solution of this equation is
$$|\psi(t) \rangle = \exp \left (-\frac{\mathrm{i}}{\hbar} \hat{H} t \right ) |\psi(0) \rangle.$$
Now suppose at the initial time ##t=0## the system has been prepared in an energy eigenstate ##|\psi(0) \rangle=|u_E \rangle##. Then you have [corrected in view of #6]
$$|\psi(t) \rangle = \exp \left (-\frac{\mathrm{i}}{\hbar} \hat{H} t \right ) |u_E \rangle =
\exp \left (-\frac{\mathrm{i}}{\hbar} E t \right ) |u_E \rangle,$$
which means that the time evolution of the state is just given by multiplying ##|u_E \rangle## with a phase factor.
Now since the physical content of the state ket is just given by Born's rule, i.e., if you have a complete set of compatible observables ##A_i## (##i \in \{1,\ldots,n \}##) and ##|a_1,\ldots,a_n \rangle## a complete set of common eigenvectors of the corresponding essentially self-adjoint operators ##\hat{A}_i##, then the probability (distribution) to measure at time ##t## the values ##(a_i)## is given by
$$P_{\psi}(t,a_1,\ldots,a_n) = |\langle a_1,\ldots a_n |\psi(t) \rangle|^2,$$
a phase factor doesn't play any role, i.e., for the above case that an energy eigenstate is the initial state of the system, the probability distribution is time-independent
$$P_{\psi}(t,a_1,\ldots a_n)=|\langle a_1,\ldots a_n |u_E \rangle|^2.$$
Thus, the eigenvectors of the Hamiltonian represent stationary states of the quantum system.