I consider the relation
$$E = \sqrt{ \vec{p}^2 c^2 + (m c^2)^2}$$
by far not "banal" because it reflects the invariance of the 4-vector length in the Minkowski-metric. I will be more explicit:
Energy and momentum can be grouped together into a what is called in special relativity 4-vector, here called 4-momentum vector $p^{\mu} = (p^0,p^1,p^2,p^3)= (E/c,\vec{p})$.
As the squares of 3-vectors in euclidian space don't change -- i.e. are invariant under 3-dimensional rotations like in the Newton's mechanics
$\vec{p}^2= p^2_x+ p^2_y + p^2_z$.
squares $X^2$ of 4-vectors $X^{\mu}$ in Minkowski space are invariant under "rotations" of the latter -- i.e. invariant under Lorentz-transformations:
$$p^2:=\sum\limits_{\mu=0}^3 p^{\mu}p_{\mu} = (\frac{E}{c})^2 - \vec{p}^2 = m^2c^2$$
the "banal" term $ (mc)^2 $ does not change under Lorentz-transformations which actually shows that the definition $p^{\mu} = (p^0,p^x,p^y,p^z)= (E/c,\vec{p})$ makes sense.
In Euclidian space the 4-momentum $p^{\mu}p_{\mu}$ would only contain positive terms, however, in the Minkowski-space, in which Lorentz-transformations operate
however the length of 4-vector contains a minus sign due to its special metric and makes the introduction of 2 types of 4-vectors reasonable:
$p^{\mu} =(E/c,\vec{p})$ and $p_{\mu} =(E/c,-\vec{p})$
The evaluation of
$\sum\limits_{\mu=0}^3 p^{\mu}p_{\mu}$
according to the given definitions exactly yields the equation above.
It could be argued that it is just a stronger formalisation of the same relation. However, the importance of invariant quantities in relativistic physics cannot be stressed enough. I give another example
The same concept applies to $x^{\mu}=(t,\vec{x})$ (and $x_{\mu}=(t,-\vec{x})$ ) which can also be recognized as 4-vector whose 4-length $x^{\mu} x_{\mu}$ respectively $\Delta x^{\mu}\Delta x_{\mu}$ is invariant under Lorentz-transformations. Einstein discovered in particular this invariance and constructed from this recognition the Lorentz-transformations.
$(\Delta s)^2\equiv \Delta x^{\mu}\Delta x_{\mu}= c^2 (t_2 - t_1)^2 - (x_2-x_1)^2$
What photons concerns, so it was found that the well-known dispersion relation of light $\vec{k}^2 =(\frac{\omega}{c})^2$ in vacuum looks like
$$\sum\limits_{\mu=0}^3 k^{\mu}k_{\mu} = (\frac{\omega}{c})^2 - \vec{k}^2 = 0$$
or together with Plancks relationship $E =\hbar \omega$ and $\vec{p} = \hbar \vec{k}$ like
$$\sum\limits_{\mu=0}^3 p^{\mu}p_{\mu} = (\frac{E}{c})^2 - \vec{p}^2 = 0$$
so we recognize that $(\omega/c, \vec{k})$ can also be written as a 4-vector $k^{\mu}=(\omega/c, \vec{k})$ just with the little distinction that the 4-length of this vector $k^2: = \sum\limits_{m=0}^3 k^{\mu}k_{\mu}=0$,
which is of course also invariant under Lorentz-transformations.
The dispersion relation for light was of course known long before the discovery of special relativity, now that we see how nicely it fits into the theory with light quanta of rest mass zero, it is so natural distinguish between massless
particles with $m=0$ and $m\neq 0$. In particular the relation
$$(\frac{E}{c})^2 - \vec{p}^2 = m^2c^2$$
serves as dispersion relation for almost all elementary particles with rest mass in the standard model as does $$(\frac{\omega}{c})^2 - \vec{k}^2 = 0$$ for massless particles like photons.
EDIT:
Apart from the importance to construct 4-vectors in a way that their 4-length (squared) is Lorentz invariant, the definition of the energy is chosen to fit in a 4-vector so that it transforms under Lorentz transformation as the $0^{th}$ component of a 4-vector. If the another definition was preferred, for instance $E_{kin}= (\gamma-1)mc^2$ (which has to property to reproduce in the Newton's limit $E_{kin}\approx\frac{1}{2}mv^2$ ) we don't know easily how to transform it under Lorentz transformations. However, this feature is important if in scattering experiments one wants to switch quickly from the CMSystem to the labor -system or vice-versa.