We have two equivalent expressions for the four-momentum,
$$\textbf{P} = \gamma m (c,\vec{u}) \tag{1}\label{1}$$
$$\textbf{P} = (E/c,\vec{p}) \tag{2}\label{2}$$
where $\gamma = 1/\sqrt{1-(\frac{u}{c})^2}$. The four-force and the three-force are given by,
$$\textbf{F} = \frac{d \textbf{P}}{d\tau} \tag{3}\label{3}$$
$$\vec{f} = \frac{d \vec{p}}{dt} \tag{4}\label{4}$$
so that the four-force can be written as,
$$\textbf{F} = \frac{d}{d\tau}(E/c,\vec{p}) = \gamma \Big(\frac{dE}{cdt},\frac{d\vec{p}}{dt}\Big) = \gamma \Big(\frac{\vec{f} \cdot \vec{u}}{c},\vec{f}\Big)\tag{5}\label{5}.$$
There are several ways to derive,
$$\frac{dE}{dt} = \vec{f} \cdot \vec{u} \tag{6}\label{6}$$
to which I won't discuss the other ways on how to do it. I'm going to derive it the direct way, using the relativistic energy $E = \gamma mc^2$. So, $\frac{dE}{dt} = mc^2 \frac{d\gamma}{dt}$ (in this post I'm assuming a constant mass).
$$\frac{d\gamma}{dt} = \frac{\gamma^3}{c^2} u \frac{du}{dt} = \frac{\gamma^3}{c^2} \vec{u} \cdot \frac{d\vec{u}}{dt} \tag{7}\label{7}$$
$$\frac{dE}{dt} = \gamma^3 m \vec{u} \cdot \frac{d\vec{u}}{dt} \tag{8}\label{8}$$
The question now is how can we establish $\gamma^3 m \frac{d\vec{u}}{dt} = \vec{f}$? We need to have a form $\gamma^3 m \frac{d\vec{u}}{dt} = \frac{d(\gamma m\vec{u})}{dt} = \vec{f}$. So we study this form,
$$\frac{d(\gamma m\vec{u})}{dt} = \gamma m \frac{d\vec{u}}{dt} + m\vec{u} \frac{d\gamma}{dt} = \gamma m \frac{d\vec{u}}{dt} + \frac{\gamma^3 m}{c^2} \vec{u} \Big(\vec{u} \cdot \frac{d\vec{u}}{dt}\Big). \tag{9}\label{9}$$
Using vector algebra identity $A \times (B \times C) = (A \cdot C)B - (A \cdot B)C$, where $A=B=\vec{u}$ and $C=\frac{d\vec{u}}{dt}$,
$$\frac{d(\gamma m\vec{u})}{dt} = \gamma m \frac{d\vec{u}}{dt} + \frac{\gamma^3 m}{c^2} \vec{u} \Big(\vec{u} \cdot \frac{d\vec{u}}{dt}\Big) = \gamma m \frac{d\vec{u}}{dt} + \frac{\gamma^3 m}{c^2} \Bigg[(\vec{u} \cdot \vec{u}) \frac{d\vec{u}}{dt} + \vec{u} \times \Big(\vec{u} \times \frac{d\vec{u}}{dt}\Big)\Bigg]. \tag{10}\label{10}$$
If I assume $\vec{u} \times \frac{d\vec{u}}{dt} = 0$, then
$$\frac{d(\gamma m\vec{u})}{dt} = \gamma m \frac{d\vec{u}}{dt} + \frac{\gamma^3 m}{c^2} (\vec{u} \cdot \vec{u}) \frac{d\vec{u}}{dt} = \gamma m \frac{d\vec{u}}{dt} + \frac{\gamma^3 m}{c^2} u^2 \frac{d\vec{u}}{dt} = \gamma^3 m \frac{d\vec{u}}{dt}, \tag{11}\label{11}$$
the calculation of the equality in $\eqref{11}$ is plain algebra. Thus we have confirmed the form we needed in the statement above $\eqref{9}$. However, this can only be true if the three-force and velocity are completely aligned, but $\eqref{6}$ is a general form which does not say anything about the orientation of $\vec{f}$ and $\vec{u}$. Am I missing something here or is this really the case?