Let's set $c=1$ for simplicity.

Using your observations, it suffices  to show that (just combine the second and third equations you write down)
$$
  \dot \gamma = \vec u \cdot \frac{d}{dt}(\gamma \vec u).
$$
To prove this, the following facts are useful:
$$
  \dot \gamma = \gamma^3\vec u \cdot\dot{\vec u}, \qquad \gamma^2\vec u^2 +1 = \gamma^2. 
$$
Now just compute
\begin{align}
\vec u \cdot \frac{d}{dt}(\gamma \vec u)
  &=\vec u \cdot (\dot \gamma \vec u + \gamma \dot{\vec u}) \\
&=\vec u \cdot (\gamma^3(\vec u \cdot \dot{\vec u})\vec u + \gamma \dot{\vec u}) \\
&= \gamma \vec u \cdot \dot{\vec u}(\gamma^2 \vec u^2 + 1) \\
&= \gamma^3\vec u \cdot\dot {\vec u} \\
&= \dot \gamma
\end{align}