Nirenberg's paper On elliptic PDEs claims that if a function $f$ on $\mathbb{R}^n$ tends to zero at infinity or is in $L^q$ for any $q < \infty$ then the "interpolation" inequality $$ \lVert∇ f \rVert_{L^{2p}} ≤ C \left(\lVert f\rVert_{L^\infty} \lVert ∇^2 f \rVert_{L^p}\right)^{1/2} $$ holds (theorem on p. 125). This is proved by establishing the inequality $$ \lVert∇ f \rVert_{L^{q}} ≤ C \left(\lVert f\rVert_{L^r} \lVert ∇^2 f \rVert_{L^p}\right)^{1/2} \label{1}\tag{1} $$ (equation (2.5) in the proof of that theorem) for any $1 < r,p < \infty$ and $\frac{2}{q} = \frac{1}{r} + \frac{1}{p}$ with a constant $C$ independent of $r, p$, and the taking the limit $r \to \infty$. I haven't really looked at Nirenberg's proof, but as far as I can see it is "calculus-based".
On the other hand inequality \eqref{1} looks like the interpolation inequality $$ \lVert f \rVert_{W^{1,q}} ≤ C_{p, r} \left(\lVert f\rVert_{L^r} \lVert f \rVert_{W^{p, 2}}\right)^{1/2} \label{2}\tag{2} $$ for $1 < r, p < \infty$ and $q$ as above, that you get by using harmonic analysis (the Mihlin multiplier theorem and the Littlewood--Paley characterization of Sobolev spaces, to be precise) to establish that $W^{1,q}$ is the $\frac{1}{2}$-interpolation space between $W^{2,p}$ and $W^{0,r} = L^r$ (e.g. Bergh--Löfström, An introduction to interpolation spaces, Thm 6.4.5).
Is there some way to get \eqref{1} using the same methods as for \eqref{2} (maybe only with $\lVert f \rVert_{W^{p,2}}$ on the RHS instead of $\lVert ∇^2 f \rVert_{L^p}$), but with the constant independent of $p, r$, so that you can conclude the case $r = \infty$?
Or even better, is there some way to modify the interpolation arguments to give the $r = \infty$ case directly?