If we consider the the relativistic Lorentz force law:

$$\frac{d}{dt} (m\gamma \vec{u})=e(\vec{E}+\vec{u} \times \vec{B})$$

How can we deduce:

$$\frac{d}{dt} (m\gamma c^2)=e \vec{E} \cdot \vec{u}$$

Clearly dotting with $\vec{u}$ will give us the RHS.
Which leaves us:

$$\vec{u} \cdot \frac{d}{dt} (m\gamma \vec{u})=e \vec{u} \cdot \vec{E}$$

Could anyone help explain how to proceed and if this is the correct method?